ORIGINAL RESEARCH article

Front. Genet., 18 June 2025

Sec. Genetics of Common and Rare Diseases

Volume 16 - 2025 | https://doi.org/10.3389/fgene.2025.1516872

Identification and functional characterization of ABCA4 gene variants in three patients with Stargardt disease or retinitis pigmentosa

  • QL

    Qi Luo 1

  • JH

    Juan Huang 2

  • LS

    Lu Shi 3,4

  • GZ

    Guanghong Zhang 1

  • LX

    Linping Xue 2

  • KW

    Kehu Wu 1

  • XL

    Xiaoyu Li 1

  • LY

    Lei Yang 1

  • DL

    Dujun Li 1

  • LM

    Liangwei Mao 4

  • JL

    Jihong Luo 1*

  • 1. Department of Ophthalmology, Hubei Provincial Hospital of Traditional Chinese Medicine, Affiliated Hospital of Hubei University of Chinese Medicine, Hubei Province Academy of Traditional Chinese Medicine, Wuhan, China

  • 2. The First Clinical Medical School of Hubei University of Chinese Medicine, Wuhan, China

  • 3. Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China

  • 4. Department of Medicine, Wuhan Primbio Medical Laboratory, Wuhan, China

Abstract

Introduction:

The diversity of phenotypes, ranging from inherited retinal dystrophies (such as Stargardt disease 1, cone–rod dystrophy 3, and retinitis pigmentosa 19) to late-onset age-related macular degeneration 2, has been attributed to loss-of-function variants in the ABCA4 gene. In this study, we aimed to identify and analyze potential pathogenic ABCA4 variants in patients with Stargardt disease or retinitis pigmentosa and to explore the impact of an intronic variant (NM_000350.3:c.6386 + 4A>G) on mRNA splicing.

Methods:

We enrolled three patients from unrelated families with Stargardt disease or retinitis pigmentosa after comprehensive ophthalmological evaluations were performed. Whole-exome sequencing and Sanger sequencing were applied for mutation screening, focusing on inherited retinal dystrophy-related genes. Additionally, the splicing alteration caused by c.6386 + 4A>G was functionally characterized by a minigene splicing assay.

Results:

Five ABCA4 germline variants were detected in three patients: one frameshift, one nonsense, one splicing, and two missense variants. Furthermore, two pathogenic and two likely pathogenic variants and one variant of uncertain significance were determined according to ACMG/AMP and ClinGen sequence variant interpretation (SVI) guidelines. The minigene splicing assay result proved that c.6386 + 4A>G affected the wild-type donor splice-site recognition of intron 46 and yielded a truncated transcript with a 47-bp deletion in exon 46.

Discussion:

Our study identified two novel ABCA4 variants, expanding the mutational spectrum of the ABCA4 gene in Stargardt disease and retinitis pigmentosa while providing new insights into the molecular pathology of ABCA4 splicing defects.

Introduction

Inherited retinal dystrophies (IRDs) are a group of hereditary diseases that cause severe vision loss due to photoreceptor degeneration or dysfunction (; ; ). They affect approximately 1 in 2,000–3,000 individuals worldwide (). IRDs are extremely heterogeneous in terms of genotype and phenotype. One phenotype could result from variants in multiple genes, and diverse variants in the same gene may lead to different phenotypic characteristics (; ). A total of 306 identified genes have been documented in the Retinal Information Network (RetNet, https://web.sph.uth.edu/RetNet/; accessed October 2024). ABCA4 is one of the most prevalent IRD-causing genes across populations, with its variants accounting for 2%–5% of non-syndromic retinitis pigmentosa (RP) cases (; ). For instance, Haer-Wigman et al. reported that approximately 7% of the causative variants in 266 Dutch patients with visual impairment were detected in the ABCA4 gene (). Bernardis et al. found that the most recurrent pathogenic mutations in 109 IRD-affected patients were in the ABCA4 and USH2A genes (). Perea-Romero et al. identified ABCA4 as the most prevalent gene in autosomal recessive nonsyndromic RP families among 6089 IRD-affected cases ().

The human ABCA4 gene (OMIM: 601691) is located on chromosome 1p22.1, spans more than 128 kb, and contains 50 exons (). The protein encoded by this gene, ATP-binding cassette superfamily transmembrane protein, is primarily expressed in retinal photoreceptors. The ABCA4 protein transports N-retinylidene-phosphatidylethanolamine (NrPE) from the lumen to the cytoplasmic leaflet of rod and cone disc membranes by utilizing the energy derived from adenosine triphosphate (ATP) binding and hydrolysis (). If the ABCA4 protein is dysfunctional, NrPE and all-trans-retinal (ATR) accumulate and condense into N-retinylidene-N-retinylphosphatidylethanolamine (A2PE) in the photoreceptor outer segments. Subsequently, the shed photoreceptor outer segments are phagocytosed by retinal pigment epithelial (RPE) cells, which hydrolyze A2PE via lysosomal enzymes to form N-retinylidene-N-retinylethanolamine (A2E) (; ). Consequently, A2E that cannot be further metabolized accumulates within the RPE, forming the main component of lipofuscin, a substance toxic to the RPE, leading to degeneration of RPE and loss of photoreceptor cells (). Clinically, this process manifests as a reduction or loss of visual acuity.

Functional impairment of the ABCA4 protein is responsible for a wide range of retinal degeneration phenotypes, including Stargardt disease 1 (STGD1), cone–rod dystrophy 3 (CRD3), retinitis pigmentosa 19 (RP19), and other conditions collectively referred to as ABCA4-associated retinopathies (Sheremet et al., 2018; ). Among these, RP19 is characterized by the progressive degeneration and death of rod photoreceptors, resulting in night blindness, decay of visual acuity, and concentric constriction of the visual field (; ). In RP19 patients, fundus examinations reveal bone spicule-like pigmentation, peripheral scattered pigmentation, or severe atrophy of the RPE, while the electroretinogram (ERG) shows abolished rod responses and markedly diminished cone responses (). The main clinical manifestation of STGD1 is progressive bilateral centrifugal vision loss caused by macular atrophy and subretinal deposition of lipofuscin-like substances, typically occurring in childhood or early adolescence (; ). It is crucial to report and functionally characterize novel potentially pathogenic ABCA4 sequence variants to enable precise diagnosis and therapy of ABCA4-related IRDs.

We hereby describe three unrelated patients who were clinically diagnosed with RP or STGD1. We analyzed the mutational spectrum of ABCA4 in the patients using whole-exome sequencing (WES) and explored the pathogenesis of the novel splicing variant by an in vitro minigene assay.

Materials and methods

Clinical examination

This study was conducted in accordance with guidelines of the Declaration of Helsinki. Ethical approval was obtained from the Human Research Ethics Committee of Hubei Provincial Hospital of Traditional Chinese Medicine (Wuhan, China). Each proband underwent a detailed clinical history review and comprehensive ophthalmological evaluation by a retina specialist. Age at onset and disease progression was determined by patients’ self-reported history of visual loss. The following examinations were performed: fundus photography (Canon CF-1, Tokyo, Japan), visual field examination (AP-7000 automatic perimeter, Kowa, Japan), spectral domain optical coherence tomography (SD-OCT; Heidelberg Engineering GmbH 69121, Heidelberg, Germany), and fluorescein fundus angiography (Canon, Tokyo, Japan). The patients’ parents denied any family history of ocular disease.

DNA isolation and next-generation sequencing

Three milliliters of peripheral blood was collected from the three patients and their available family members. Genomic DNA (gDNA) was isolated from peripheral blood leukocytes using the magnetic bead method blood gDNA extraction kit (Tiangen Biotech, Beijing, China). The gDNA concentration was quantified using a Qubit 3.0 fluorometer (Qubit double-stranded DNA detection kit, Invitrogen, USA), while integrity was assessed by 1% agarose gel electrophoresis. Qualified gDNA samples were fragmented using the VAHTS® Universal Plus DNA Library Prep Kit (Vazyme, Nanjing, China) and processed with the xGen® Exome Research Panel v2.0 (Integrated DNA Technologies, USA) for liquid-phase capture to build a library targeting the human whole-exome region. The prepared DNA library was then used to perform 150-bp paired-end sequencing of the whole exome on the MGISEQ-T7 sequencing platform (Mgi Tech Co., Ltd., China).

Bioinformatics analysis

Bioinformatics analysis revealed that each sample yielded over 10 Gb of sequencing data, with coverage over 20× at more than 99% and an average read depth exceeding 120× across all analyzed genomic regions. Sequencing reads were aligned to the GRCh38 human reference genome using Burrows–Wheeler Aligner (BWA, v0.7.16a, http://bio-bwa.sourceforge.net/). To ensure alignment accuracy, low-quality reads and adapter sequences were removed using fastp (fastp v0.21.0). Polymerase chain reaction (PCR) duplicates were eliminated with sambamba (sambamba 0.8.0). Variant calling was performed on the generated BAM file using the Genome Analysis Toolkit (GATK, v4.0.8.1, https://gatk.broadinstitute.org/). Detected variants included single-nucleotide variants (SNVs) and small insertions/deletions (Indels) less than 50 bp in length. Variants were filtered based on sequencing depth (>20X) and variant quality (QualByDepth>2). ANNOVAR (v20200608, https://annovar.openbioinformatics.org/) was utilized for variant annotation. Annotated variants were further filtered and prioritized by minor allele frequency (ExAC, gnomAD, 1000 Genomes, and in-house variant frequency database), protein effect (missense, synonymous, frameshift, in-frame indel, and splicing variants), various in silico prediction scores (e.g., CADD_PHRED, REVEL, SIFT, Polyphen-2, MutationTaster, and SpliceAI), inheritance pattern of related disease (autosomal dominant, autosomal recessive, or X-linked), and clinical phenotype concordance. Exonic copy number variants (CNVs) were calculated using an in-house algorithm by analyzing target region read depths relative to an internal control dataset. Cryptic splice sites were predicted using Human Splice Finder (HSF, v3.1, http://www.umd.be/HSF3/HSF.shtml), VarSEAK (https://varseak.bio/index.php), and SpliceAI (https://spliceailookup.broadinstitute.org/). The pathogenicity of candidate variants was interpreted according to the American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) (; ) and ClinGen sequence variant interpretation (SVI) guidelines (https://clinicalgenome.org/working-groups/sequence-variant-interpretation/). The structural domain annotation of identified variants was performed using InterPro (https://www.ebi.ac.uk/interpro). Variants were mapped to functional domains based on the InterPro and classified according to their positions within key domains. To assess the evolutionary conservation of amino acid residues affected by missense variants, multiple sequence alignment was performed using ClustalW (https://www.genome.jp/tools-bin/clustalw).

Sanger sequencing

Specific primer pairs for the candidate variants were designed using the Primer3 algorithm (Primer3web version 4.1.0, https://bioinfo.ut.ee/primer3/). The designed primers were used to amplify regions of DNA sequence containing the variants by PCR. The PCR amplicons were purified, and then Sanger sequencing was performed for PCR products on the ABI 3730XL DNA analyzer platform (Thermo Fisher Scientific, Waltham, United States). The nucleotide primer sequences of candidate variants are listed in Table 1.

TABLE 1

VariantPrimerPrimer sequence (5′–3′)Melts degreeProduct size (bp)
c.6386 + 4A>GABCA4-chr1-94000998-Fatc​tgc​tgc​tgg​aga​atg​ag60°C371
ABCA4-chr1-94000998-Racg​tca​tcg​tga​gca​tca​tc60°C371
c.3286_3287delABCA4-chr1-94042802-Fgtc​act​gat​aaa​ccc​cct​tc60°C380
ABCA4-chr1-94042802-Raaa​cca​ctg​ctg​ggt​taa​gg60°C380
c.2261T>CABCA4-chr1-94056722-Fctt​tcc​agg​gca​gat​gaa​tc60°C379
ABCA4-chr1-94056722-Rtca​tgc​tgt​gct​ttc​tgc​tc60°C379
c.2424C>GABCA4-chr1-94055274-Fggt​acc​aag​cga​gta​agc​cat60°C307
ABCA4-chr1-94055274-Raca​ctg​aga​gac​aat​ggg​gac59°C307
c.1901T>CABCA4-chr1-94062613-Fgag​ctt​ggg​aga​agc​agg​tt60°C374
ABCA4-chr1-94062613-Rtga​gct​atc​caa​gcc​cgt​tc59°C374

The primer sequences of candidate variants.

Minigene splicing assay

The c.3286_3287del variant has been previously reported; therefore, functional studies were not conducted on it. We utilized a minigene splicing assay to assess the impact of a novel intronic ABCA4 variant (c.6386 + 4A>G) on mRNA splicing. In vitro minigene functional verification was performed by Bioeagle Biotech Co., Ltd. (Wuhan, China). Wild-type and mutant-type ABCA4 genomic segments encompassing exons 45 to 47 were cloned and inserted into the splicing pcDNA3.1 vector to construct pcDNA3.1-ABCA4-wt and pcDNA3.1-ABCA4-mut minigenes, respectively. First, two pairs of nested primers (118564-ABCA4-F/120474-ABCA4-R and 118799-ABCA4-F/120730-ABCA4-R) were designed to PCR amplify the ABCA4 genomic region encompassing exons 45–47. The 1181-bp wild-type fragment was amplified from the nested PCR product using primers pcDNA3.1-ABCA4-HindIII-F and pcDNA3.1-ABCA4-XhoI-R. The anterior and posterior segments of the mutant fragment were separately amplified using the primers pcDNA3.1-ABCA4-HindⅢ-F ABCA4-mut-R and pcDNA3.1-ABCA4-XhoI-R ABCA4-mut-F, respectively, followed by equimolar mixing of the products. The full-length 1181-bp mutant fragment was PCR-amplified from this mixture using the primers pcDNA3.1-ABCA4-HindIII-F and pcDNA3.1-ABCA4-XhoI-R. The wild-type and mutant fragments were cloned and inserted into the empty pcDNA3.1 vector (Figure 1). MCF-7 and HEK293T cells were transfected with the minigene constructs using Hieff Trans liposomal transfection reagent (Yeasen Biotechnology Shanghai Co., Ltd., Cat. China) per the manufacturer’s protocol. After 48-h incubation, total cellular RNA was isolated using TRIzol reagent and reverse-transcribed. The cDNA was PCR-amplified with vector primers pcDNA3.1-F and pcDNA3.1-R. Sanger sequencing was used to analyze the splicing pattern of the wild-type and mutant minigene vectors (with parallel validation in the pcMINI-N vector; see Supplementary Figure S1). The primer sequences are listed in Table 2.

FIGURE 1

TABLE 2

PrimerPrimer sequence (5′–3′)Melts degreeProduct size (bp)
118,564-ABCA4-Fgca​gca​atg​gct​gtg​ttc​ac592,167
120,730-ABCA4-Raca​atg​tgt​ggg​tac​agg​at602,167
118,799-ABCA4-Ftcc​att​tcc​cag​gcg​gct​ac601,676
120,474-ABCA4-Rcag​gct​ctg​cat​gca​tcc​tg561,676
pcDNA3.1-ABCA4-HindIII-Fact​taa​gct​tat​ggt​tgc​aaa​ctg​gag​tat​taa591,181
pcDNA3.1-ABCA4-XhoI-Rtag​act​cga​gtt​tgg​act​tga​gat​gct​gaa​t621,181
ABCA4-mut-Fcac​atc​cca​cag​gca​gga​gat​tcc​cag​ggc​t68-
ABCA4-mut-Ragc​cct​ggg​aat​ctc​ctg​cct​gtg​gga​tgt​g66-
pcDNA3.1-Fcta​gag​aac​cca​ctg​ctt​ac53434
pcDNA3.1-Rtag​act​cga​gtt​tgg​act​tga​gat​gct​gaa​t62434

The primer sequences of the minigene splicing assay.

Results

Clinical features

Three unrelated Han Chinese patients (two males and one female) with RP or Stargardt Disease. underwent WES. The basic clinical phenotypic features and molecular diagnoses of these three patients are shown in Table 3. Bilateral visual acuity in each patient showed a progressive decline. No family history of RP, macular degeneration, or other retinal dystrophies was reported in any of the patients.

TABLE 3

Patient IDAge (years)SexAge of onset (years)BCVA (od/os)Exon/IntronVariantsVariant typeZygote typePathogenicityMutation statusAllele segregation testPhenotypeInitial clinical diagnosis
Patient 138Male120.02/0.02Intron46c.6386 + 4A>GSplicingHomLPNovelFalseVisual loss, visual field loss, decreased RNFL thickness, retinal vessel attenuationRP19
Patient 228Male240.03/0.04Exon15c.2261T>C
p. (Phe754Ser)
MissenseHetLPReportedTrueVisual loss, posterior pigment deposition, hypofluorescence in the macular area and posterior poleRP19
Exon22c.3286_3287del
p. (Ser1096Glufs*78)
FrameshiftHetPNovel
Patient 332Female100.1/0.02Exon16c.2424C>G
p. (Tyr808Ter)
NonsenseHetPReportedTruePunctate hyperfluorescent spots, decreased RNFL thickness, reduced visual sensitivitySTGD1
Exon13c.1901T>C
p. (Leu634Pro)
MissenseHetVUSReported

The clinical features of three patients and ABCA4 variants identified in this study.

Abbreviations: BCVA, best corrected visual acuity; od, right eye; os, left eye; Hom, homozygous; Het, heterozygous; P, pathogenic; LP, likely pathogenic; VUS, variant of uncertain significance; RP19, retinitis pigmentosa 19; STGD1, Stargardt disease 1.

Patient 1

Fundus photography of Patient 1, performed in another hospital, showed attenuated blood vessels, bone spicule-shaped pigment deposits, and optic disc pallor. Visual field perimetry disclosed significant peripheral visual field loss in both eyes. High-resolution SD-OCT of the retinal nerve fiber layer (RNFL) showed a decrease in the nasal inner retinal thickness (Figure 2).

FIGURE 2

Patient 2

Fundus photography of Patient 2 showed normal-appearing optic disks and hyperpigmentation in the macula in both eyes. Fluorescein fundus angiography (FFA) demonstrated diffuse hypofluorescence in the macular area and posterior pole bilaterally. Macular OCT revealed degeneration and thickening of the outer retina (Figure 3).

FIGURE 3

Patient 3

Fundus images of Patient 3 revealed atrophy and pigmentation of the posterior pole of the retina. FFA displayed multiple hyperfluorescent punctate lesions within the macular area. Macular OCT showed a marked decrease in the thickness of the bilateral RPE, retinal neurosensory layer, and choroid layer in the macula. The RNFL showed abnormal thinning in the superotemporal and inferotemporal quadrants around both optic disks. The visual field showed reduced visual sensitivity and defects in the central 30° bilaterally (Figure 4).

FIGURE 4

Molecular analysis and minigene assay

Patient 1 revealed a novel homozygous variant c.6386 + 4A>G in the intron 46 of the ABCA4 gene adjacent to the donor splice site (Figure 5A). This variant was absent from public databases (gnomAD v4.1.0, https://gnomad.broadinstitute.org/) and the literature. In silico analyses using SpliceAI and HSF predicted disruption of the wild-type donor splice site. The minigene splicing assay results showed that both pcDNA3.1-ABCA4-wt and pcDNA3.1-ABCA4-mut plasmids were successfully transfected into MCF-7 and HEK293T cells (Figure 5B). Agarose gel electrophoresis revealed that the mutant type had two bands. One band was the same as the wild-type band in length, and the other band was smaller with faster migration (Figure 5C). Sanger sequencing confirmed that the wild-type minigene transcribed the normal mRNA product composed of exon 45, exon 46, and exon 47. However, the mutant minigene mRNA transcribed product lacked 47 bp of exon 46 (c.6340_6386del), causing a frameshift to p. (Val2114Hisfs*5) (Figure 5D). The transcriptional results of the pc-MINI-N vector are shown in Supplementary Figure S2. The variant was classified as likely pathogenic based on ACMG/AMP and ClinGen SVI guidelines with evidence from a minigene assay (PS3+PM2_Supporting + PM3_Supporting).

FIGURE 5

Patient 2 revealed biallelic ABCA4 variants: c.3286_3287del, p.(Ser1096Glufs*78), and c.2261T>C, p.(Phe754Ser). Segregation analysis confirmed that the c.3286_3287del variant was paternally inherited, and the c.2261T>C variant was maternally inherited (Figures 6A,B). c.2261T>C is a phenylalanine-to-serine missense substitution, predicted to be deleterious by SIFT (https://www.jcvi.org/research/provean) and MutationTaster (https://sites.google.com/site/revelgenomics/). Conservation analysis showed that the mutated phenylalanine residue was conserved across seven species, spanning vertebrates from mammals to amphibians (Figure 7). The variant was previously reported in three Chinese patients with STGD1 or cone–rod dystrophy (; ; ) and was categorized as likely pathogenic (PM2_Supporting + PM3_Strong + PP3). The c.3286_3287del frameshift variant was predicted to undergo nonsense-mediated mRNA decay by the AutoPVS1 tool (https://autopvs1.genetics.bgi.com). It was classified as pathogenic (PVS1+PM2_Supporting + PM3).

FIGURE 6

FIGURE 7

Patient 3 inherited the c.2424C>G, p.(Tyr808Ter) variant from her mother and the c.1901T>C, p.(Leu634Pro) variant from her father (Figures 6C,D). The c.2424C>G variant creates a premature stop codon at position 808. Sanger sequencing confirmed that the nonsense variant was also present in the patient’s unaffected son. This variant was previously reported in over 10 STGD1 or cone-rod dystrophy patients with compound heterozygous genotypes in the ABCA4 gene (; ; ; ). It was classified as a pathogenic variant (PVS1+PM3_Verystrong + PM2_Supporting). The c.1901T>C variant results in a leucine-to-proline missense substitution. Multiple sequence alignment showed that codon 634 was conserved across seven species (Figure 7), although the variant was classified as a variant of uncertain significance (PM2_Supporting + PM3+PP3) based on limited evidence to date. Population frequencies and in silico prediction results of five ABCA4 variants in the three patients are provided in Supplementary Table S1.

Discussion

ABCA4 is associated with STGD1, CORD3, RP19, and AMD2 (; ; ). Large cohort studies across diverse ethnic populations have defined the ABCA4 mutational spectrum, associated phenotypes, and genotype–phenotype correlations in these IRDs (; ; ; ; ; ). Pathogenic variants encompass various types, including missense and nonsense variants, splicing (canonical and non-canonical) variants, deep-intronic variants, in-frame indels, large indels, and complex rearrangements (Human Gene Mutation Database; https://www.hgmd.cf.ac.uk/ac/search.php). Several pathogenic variants exhibit marked ethnic predilection, with notably elevated frequencies in specific populations: c.5882G>A (20.5% allele frequency in Europeans) (), c.6320G>A (19.32% in African Americans) (), and c.2424C>G (5.2% in Chinese) (). Further mechanistic studies of these variants may elucidate the molecular basis of the relationship between genotypic and phenotypic manifestations.

In this study, we describe three middle-aged patients with ABCA4-associated IRDs. Patient 1 had the homozygous c.6386 + 4A>G splicing variant and was born to consanguineous parents. Although parental samples were unavailable, we infer that his father and mother were both heterozygous carriers of this variant. Most disease-causing splicing variants occur at canonical +1/+2 or −1/−2 positions (); however, non-canonical splice-site variants can also disrupt mRNA processing and result in a dysfunctional protein. Several studies have reported different intronic variants in the ABCA4 gene associated with retinal degeneration. For example, Cremers et al. reported five cone-rod dystrophy (CRD) patients in a consanguineous family who were compound heterozygotes for the IVS30+1G→T variant and the IVS40+5G→A variant, which would lead to skipping of exon 40 with a reading frame shift in exon 41 (). Sangermano et al. performed in vitro minigene splice assays to analyze the effect of 47 non-canonical splice-site variants (NCSS) on ABCA4 pre-mRNA splicing. Among these variants, 29 variants (62%) induced single or multi-exon skipping (e.g., c.302 + 4A>C, c.2919-10T>C), eight variants (17%) showed multiple splice defects like skipping/elongation or skipping/retention (e.g., c.160 + 5G>C, c.303-3C>G, and c.5898+5del), six variants (13%) resulted in exon elongation (e.g., c.1937 + 13T>G, c.5313-3C>G) (). Fadaie et al. found that deep-intronic variants, such as c.193-619A>G, c.2919-826T>A, and c.3050 + 370C>T, caused significant splicing defects in minigene splice assays. These variants generated pseudoexons containing premature stop codons, leading to predicted truncated ABCA4 proteins ().

The variants c.3286_3287del, p.(Ser1096Glufs*78), and c.2261T>C, p.(Phe754Ser) in trans were identified in Patient 2, while c.2424C>G, p.(Tyr808Ter), and c.1901T>C, p.(Leu634Pro) in trans were identified in Patient 3 in ABCA4. Both patients had a loss-of-function variant that was predicted to cause protein truncation and a missense variant that was predicted to be likely to damage protein function according to sequence conservation analysis. Our structural analysis revealed that p. (Leu634Pro), p.(Phe754Ser), and p.(Tyr808Ter) were located within the ABC-2 family transporter protein domain, while p.(Ser1096Glufs*78) was mapped to the ATP-binding cassette domain.

A genotype-phenotype correlation model elucidated that residual activity of the ABCA4 protein was correlated with the severity of retinal dystrophy (van Driel et al., 1998; Maugeri et al., 1999). According to the activity of the ABCA4 protein, variants in the ABCA4 gene are classified as deleterious (no activity), severe, moderately severe, or mild. Patients harboring two severe variants or null alleles typically present with pan-retinal cone–rod dystrophy (CRD). The combination of a severe variant with a moderately severe variant leads to CRD, while a severe variant with a mild variant or two moderately severe variants results in classic Stargardt disease type 1 (STGD1) (). In our study cohort, Patient 1, who carried the homozygous c.6386 + 4A>G, p. (Val2114Hisfs*5) variant, presented with retinitis pigmentosa, a type of rod-cone dystrophy. Our minigene assay demonstrated that the variant c.6386 + 4A>G led to disruption of RNA splicing, resulting in a predicted truncated protein. Patient 2, who carried the severe variant c.3286_3287del, p.(Ser1096Glufs*78) and the moderately severe variant c.2261T>C, p.(Phe754Ser), showed a milder phenotype of retinitis pigmentosa. Patient 3, harboring two variants c.2424C>G, p.(Tyr808Ter) and c.1901T>C, p.(Leu634Pro), was diagnosed with STGD1. The variant p.(Tyr808Ter), in combination with one of the missense variants p.(Arg24His), p.(Phe2188Ser), p.(Ser100Tyr), or p.(Asn965Ser) in trans, has been previously associated with more than seven STGD1 patients (; ).

WES is a powerful and widely used technique to diagnose genetic abnormalities, yet it has limitations. Because WES primarily targets exonic regions and their immediate flanking sequences, it frequently misses deep-intronic variants that may affect splicing. This is particularly relevant for genes like ABCA4, where numerous pathogenic deep-intronic variants have been identified (; ). Furthermore, WES often detects a substantial number of variants of uncertain significance (VUS) during genetic screening of human exons and their flanking regions (; ). To interpret the clinical significance of VUS, further research and evidence are needed. Functional studies play a pivotal role in variant assessment, with minigene assays providing a rapid tool to evaluate effects on splicing compared to wild-type sequences (). Minigene assays have been applied across multiple genetic disorders. For example, Valeria et al. demonstrated efficient characterization of NF1 splicing alterations using minigene assays (). María et al. performed functional studies of ALMS1 variants with minigene assays (). Additionally, Ruixiao et al. identified seven exonic variants in SLC4A1, ATP6V1B1, and ATP6V0A4 that altered RNA splicing based on minigene assays (). Wang et al. proved that ABCA4 deep-intronic variants contributed to nearly half of unsolved Stargardt cases through minigene experiments (). While in silico predictions provide clues, experimental validation is key for the clinical interpretation of VUS (). Integrating evidence from sequencing, predictive models, and functional assays will advance the identification of disease-causing variants and elucidate underlying mechanisms of pathogenesis. These findings may contribute to the interpretation of variants of uncertain significance and highlight the potential of functional assays, such as minigene systems, in supporting variant classification. While further validation in larger cohorts is needed, this approach could aid in improving genetic counseling and personalized management.

In summary, we described three patients with ABCA4-related IRD phenotypes harboring five different variants in the ABCA4 gene. The generated data may have implications for genetic counseling, disease management, and support early intervention for a better prognosis in ABCA4-related IRD.

Statements

Data availability statement

The whole exome sequencing data reported in this paper have been deposited in the Genome Sequence Archive (Genomics, Proteomics & Bioinformatics 2021) in National Genomics Data Center (Nucleic Acids Res 2022), China National Center for Bioinformation / Beijing Institute of Genomics, Chinese Academy of Sciences (GSA-Human: HRA011833) that are publicly accessible at https://ngdc.cncb.ac.cn/gsa-human. The patients' phenotypes and detected variants have been submitted to ClinVar (https://www.ncbi.nlm.nih.gov/clinvar/) under the submission numbers SCV004174835, SCV004174061, SCV004174062, SCV004174063, SCV004174064.

Ethics statement

The studies involving humans were approved by the Human Research Ethics Committee of Hubei Provincial Hospital of Traditional Chinese Medicine. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants’ legal guardians/next of kin. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author contributions

QL: writing – original draft and writing – review and editing. JH: writing – original draft and writing – review and editing. LS: methodology, software, and writing – review and editing. GZ: data curation, investigation, and writing – review and editing. LX: data curation, investigation, and writing – review and editing. KW: formal analysis and writing – review and editing. XL: investigation, software, and writing – review and editing. LY: methodology and writing – review and editing. DL: conceptualization and writing – review and editing. LM: methodology and writing – review and editing. JL: conceptualization, funding acquisition, project administration, resources, supervision, and writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by the Knowledge Innovation Special Dawn Program Project of the Wuhan Science and Technology Bureau (No. 2022020801020510) and the Science and technology research project of the Hubei Provincial Department of Education (No. B2021108).

Acknowledgments

We would like to express our gratitude to Dr. Guohua Yang of the School of Basic Medical Science of Wuhan University for providing technical support. We also thank the patients for participating in this study.

Conflict of interest

The 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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2025.1516872/full#supplementary-material

References

  • 1

    AbramowiczA.GosM. (2018). Splicing mutations in human genetic disorders: examples, detection, and confirmation. J. Appl. Genet.59 (3), 253268. 10.1007/s13353-018-0444-7

  • 2

    Al-KhuzaeiS.BroadgateS.FosterC. R.ShahM.YuJ.DownesS. M.et al (2021). An overview of the genetics of ABCA4 retinopathies, an evolving story. Genes (Basel)12 (8), 1241. 10.3390/genes12081241

  • 3

    AllikmetsR.WassermanW. W.HutchinsonA.SmallwoodP.NathansJ.RoganP. K.et al (1998). Organization of the ABCR gene: analysis of promoter and splice junction sequences. Gene215 (1), 111122. 10.1016/s0378-1119(98)00269-8

  • 4

    Alvarez-SattaM.Castro-SanchezS.PousadaG.ValverdeD. (2017). Functional analysis by minigene assay of putative splicing variants found in Bardet-Biedl syndrome patients. J. Cell Mol. Med.21 (10), 22682275. 10.1111/jcmm.13147

  • 5

    BeharryS.ZhongM.MoldayR. S. (2004). N-retinylidene-phosphatidylethanolamine is the preferred retinoid substrate for the photoreceptor-specific ABC transporter ABCA4 (ABCR). J. Biol. Chem.279 (52), 5397253979. 10.1074/jbc.M405216200

  • 6

    BernardisI.ChiesiL.TenediniE.ArtusoL.PercesepeA.ArtusiV.et al (2016). Unravelling the complexity of inherited retinal dystrophies molecular testing: added value of targeted next-generation sequencing. Biomed. Res. Int.2016, 6341870. 10.1155/2016/6341870

  • 7

    BlackshawS.FraioliR. E.FurukawaT.CepkoC. L. (2001). Comprehensive analysis of photoreceptor gene expression and the identification of candidate retinal disease genes. Cell107 (5), 579589. 10.1016/s0092-8674(01)00574-8

  • 8

    BouzidiA.CharouteH.CharifM.AmalouG.KandilM.BarakatA.et al (2022). Clinical and genetic spectrums of 413 North African families with inherited retinal dystrophies and optic neuropathies. Orphanet J. Rare Dis.17 (1), 197. 10.1186/s13023-022-02340-7

  • 9

    ChenL.WangN.LaiM.HouF.HeJ.FanX.et al (2022). Clinical and genetic investigations in Chinese families with retinitis pigmentosa. Exp. Biol. Med. (Maywood)247 (12), 10301038. 10.1177/15353702221085711

  • 10

    CremersF. P.van de PolD. J.van DrielM.den HollanderA. I.van HarenF. J.KnoersN. V.et al (1998). Autosomal recessive retinitis pigmentosa and cone-rod dystrophy caused by splice site mutations in the Stargardt's disease gene ABCR. Hum. Mol. Genet.7 (3), 355362. 10.1093/hmg/7.3.355

  • 11

    CremersF. P. M.LeeW.CollinR. W. J.AllikmetsR. (2020). Clinical spectrum, genetic complexity and therapeutic approaches for retinal disease caused by ABCA4 mutations. Prog. Retin Eye Res.79, 100861. 10.1016/j.preteyeres.2020.100861

  • 12

    Del Pozo-ValeroM.Riveiro-AlvarezR.Blanco-KellyF.Aguirre-LambanJ.Martin-MeridaI.IancuI. F.et al (2020). Genotype-phenotype correlations in a Spanish cohort of 506 families with biallelic ABCA4 pathogenic variants. Am. J. Ophthalmol.219, 195204. 10.1016/j.ajo.2020.06.027

  • 13

    FahimA. T.DaigerS. P.WeleberR. G. (1993). “Nonsyndromic retinitis pigmentosa overview,” in GeneReviews((R)). Editors AdamM. P.FeldmanJ.MirzaaG. M.PagonR. A.WallaceS. E.AmemiyaA. (Seattle (WA)).

  • 14

    FujinamiK.ZernantJ.ChanaR. K.WrightG. A.TsunodaK.OzawaY.et al (2013). ABCA4 gene screening by next-generation sequencing in a British cohort. Invest. Ophthalmol. Vis. Sci.54 (10), 66626674. 10.1167/iovs.13-12570

  • 15

    Haer-WigmanL.van Zelst-StamsW. A.PfundtR.van den BornL. I.KlaverC. C.VerheijJ. B.et al (2017). Diagnostic exome sequencing in 266 Dutch patients with visual impairment. Eur. J. Hum. Genet.25 (5), 591599. 10.1038/ejhg.2017.9

  • 16

    HarrisonS. M.BieseckerL. G.RehmH. L. (2019). Overview of specifications to the ACMG/AMP variant interpretation guidelines. Curr. Protoc. Hum. Genet.103 (1), e93. 10.1002/cphg.93

  • 17

    HuF. Y.LiJ. K.GaoF. J.QiY. H.XuP.ZhangY. J.et al (2019). ABCA4 gene screening in a Chinese cohort with Stargardt disease: identification of 37 novel variants. Front. Genet.10, 773. 10.3389/fgene.2019.00773

  • 18

    JaffalL.JoumaaH.NoureldineJ.BanjakM.IbrahimM.MradZ.et al (2023). The genetic landscape of inherited retinal dystrophies in Arabs. BMC Med. Genomics16 (1), 89. 10.1186/s12920-023-01518-7

  • 19

    JiangF.PanZ.XuK.TianL.XieY.ZhangX.et al (2016). Screening of ABCA4 gene in a Chinese cohort with Stargardt disease or cone-rod dystrophy with a report on 85 novel mutations. Invest. Ophthalmol. Vis. Sci.57 (1), 145152. 10.1167/iovs.15-18190

  • 20

    KaraliM.BanfiS. (2015). Inherited Retinal Dystrophies: the role of gene expression regulators. Int. J. Biochem. Cell Biol.61, 115119. 10.1016/j.biocel.2015.02.007

  • 21

    KhanM.CornelisS. S.KhanM. I.ElmelikD.MandersE.BakkerS.et al (2019). Cost-effective molecular inversion probe-based ABCA4 sequencing reveals deep-intronic variants in Stargardt disease. Hum. Mutat.40 (10), 17491759. 10.1002/humu.23787

  • 22

    KrawczakM.ReissJ.CooperD. N. (1992). The mutational spectrum of single base-pair substitutions in mRNA splice junctions of human genes: causes and consequences. Hum. Genet.90 (1-2), 4154. 10.1007/BF00210743

  • 23

    LiC.LuoY.XieY.ZhangZ.LiuY.ZouL.et al (2024). Structural and functional prediction, evaluation, and validation in the post-sequencing era. Comput. Struct. Biotechnol. J.23, 446451. 10.1016/j.csbj.2023.12.031

  • 24

    LiW.HeX. D.YangZ. T.HanD. M.SunY.ChenY. X.et al (2023). De novo mutations contributes approximately 7% of pathogenicity in inherited eye diseases. Invest. Ophthalmol. Vis. Sci.64 (2), 5. 10.1167/iovs.64.2.5

  • 25

    ManleyA.MeshkatB. I.JablonskiM. M.HollingsworthT. J. (2023). Cellular and molecular mechanisms of pathogenesis underlying inherited retinal dystrophies. Biomolecules13 (2), 271. 10.3390/biom13020271

  • 26

    MaugeriA.van DrielM. A.van de PolD. J.KleveringB. J.van HarenF. J.TijmesN.et al (1999). The 2588G↠C mutation in the ABCR gene is a mild frequent founder mutation in the Western European population and allows the classification of ABCR mutations in patients with Stargardt disease. Am J Hum Genet. 64 (4), 102435. 10.1086/302323

  • 27

    MoldayR. S.GarcesF. A.ScortecciJ. F.MoldayL. L. (2022). Structure and function of ABCA4 and its role in the visual cycle and Stargardt macular degeneration. Prog. Retin Eye Res.89, 101036. 10.1016/j.preteyeres.2021.101036

  • 28

    MorbidoniV.BaschieraE.ForzanM.FuminiV.AliD. S.GiorgiG.et al (2021). Hybrid minigene assay: an efficient tool to characterize mRNA splicing profiles of NF1 variants. Cancers (Basel)13 (5), 999. 10.3390/cancers13050999

  • 29

    MullinsR. F.KuehnM. H.RaduR. A.EnriquezG. S.EastJ. S.SchindlerE. I.et al (2012). Autosomal recessive retinitis pigmentosa due to ABCA4 mutations: clinical, pathologic, and molecular characterization. Invest. Ophthalmol. Vis. Sci.53 (4), 18831894. 10.1167/iovs.12-9477

  • 30

    Olivares-GonzalezL.VelascoS.CampilloI.RodrigoR. (2021). Retinal inflammation, cell death and inherited retinal dystrophies. Int. J. Mol. Sci.22 (4), 2096. 10.3390/ijms22042096

  • 31

    Perea-RomeroI.GordoG.IancuI. F.Del Pozo-ValeroM.AlmogueraB.Blanco-KellyF.et al (2021). Genetic landscape of 6089 inherited retinal dystrophies affected cases in Spain and their therapeutic and extended epidemiological implications. Sci. Rep.11 (1), 1526. 10.1038/s41598-021-81093-y

  • 32

    QuL. H.JinX.ZengC.ZhouN. G.LiuY. H.LinY. (2021). Targeted next-generation sequencing identifies ABCA4 mutations in Chinese families with childhood-onset and adult-onset Stargardt disease. Biosci. Rep.41 (6). 10.1042/BSR20203497

  • 33

    RiveraA.WhiteK.StohrH.SteinerK.HemmrichN.GrimmT.et al (2000). A comprehensive survey of sequence variation in the ABCA4 (ABCR) gene in Stargardt disease and age-related macular degeneration. Am. J. Hum. Genet.67 (4), 800813. 10.1086/303090

  • 34

    RomanoF.LamannaF.BoonC. J. F.SiligatoA.KalraG.AgarwalA.et al (2024). Clinical, genotypic, and imaging characterization of the spectrum of ABCA4 retinopathies. Ophthalmol. Retina8 (5), 509519. 10.1016/j.oret.2023.10.023

  • 35

    RudolphG.KalpadakisP.HaritoglouC.RiveraA.WeberB. H. (2002). Mutations in the ABCA4 gene in a family with Stargardt's disease and retinitis pigmentosa (STGD1/RP19). Klin. Monbl Augenheilkd219 (8), 590596. 10.1055/s-2002-34425

  • 36

    SangermanoR.KhanM.CornelisS. S.RichelleV.AlbertS.GarantoA.et al (2018). ABCA4 midigenes reveal the full splice spectrum of all reported noncanonical splice site variants in Stargardt disease. Genome Res.28 (1), 100110. 10.1101/gr.226621.117

  • 37

    SchulzH. L.GrassmannF.KellnerU.SpitalG.RutherK.JagleH.et al (2017). Mutation spectrum of the ABCA4 gene in 335 Stargardt disease patients from a multicenter German cohort-impact of selected deep intronic variants and common SNPs. Invest. Ophthalmol. Vis. Sci.58 (1), 394403. 10.1167/iovs.16-19936

  • 38

    SheremetN. L.GrushkeI. G.ZhorzholadzeN. V., A STanasA. S.StrelnikovV. V. (2018). Inherited retinal diseases in patients with ABCA4 gene mutations. Inherited retinal diseases in patients with ABCA4 gene mutations. 134(4):6873. 10.17116/oftalma201813404168

  • 39

    SunZ.YangL.LiH.ZouX.WangL.WuS.et al (2021). Clinical and genetic analysis of the ABCA4 gene associated retinal dystrophy in a large Chinese cohort. Exp. Eye Res.202, 108389. 10.1016/j.exer.2020.108389

  • 40

    TavtigianS. V.GreenblattM. S.HarrisonS. M.NussbaumR. L.PrabhuS. A.BoucherK. M.et al (2018). Modeling the ACMG/AMP variant classification guidelines as a Bayesian classification framework. Genet. Med.20 (9), 10541060. 10.1038/gim.2017.210

  • 41

    TsangS. H.SharmaT. (2018). Stargardt disease. Adv. Exp. Med. Biol.1085, 139151. 10.1007/978-3-319-95046-4_27

  • 42

    TsybovskyY.MoldayR. S.PalczewskiK. (2010). The ATP-binding cassette transporter ABCA4: structural and functional properties and role in retinal disease. Adv. Exp. Med. Biol.703, 105125. 10.1007/978-1-4419-5635-4_8

  • 43

    UmairM. (2023). Rare genetic disorders: beyond whole-exome sequencing. J. Gene Med.25 (10), e3503. 10.1002/jgm.3503

  • 44

    van DrielM. A.MaugeriA.KleveringB. J.HoyngC. B.CremersF. P. (1998). ABCR unites what ophthalmologists divide(s). Ophthalmic Genet.19(3), 11722. 10.1076/opge.19.3.117.2187

  • 45

    VerbakelS. K.van HuetR. A. C.BoonC. J. F.den HollanderA. I.CollinR. W. J.KlaverC. C. W.et al (2018). Non-syndromic retinitis pigmentosa. Prog. Retin. Eye Res.66, 157186. 10.1016/j.preteyeres.2018.03.005

  • 46

    WangY.WangP.YiZ.OuyangJ.JiangY.LiS.et al (2025). ABCA4 deep intronic variants contributed to nearly half of unsolved Stargardt cases with a milder phenotype. Invest. Ophthalmol. Vis. Sci.66 (1), 65. 10.1167/iovs.66.1.65

  • 47

    YangY.MuznyD. M.ReidJ. G.BainbridgeM. N.WillisA.WardP. A.et al (2013). Clinical whole-exome sequencing for the diagnosis of mendelian disorders. N. Engl. J. Med.369 (16), 15021511. 10.1056/NEJMoa1306555

  • 48

    ZernantJ.CollisonF. T.LeeW.FishmanG. A.NoupuuK.YuanB.et al (2014). Genetic and clinical analysis of ABCA4-associated disease in African American patients. Hum. Mutat.35 (10), 11871194. 10.1002/humu.22626

  • 49

    ZhangR.ChenZ.SongQ.WangS.LiuZ.ZhaoX.et al (2021). Identification of seven exonic variants in the SLC4A1, ATP6V1B1, and ATP6V0A4 genes that alter RNA splicing by minigene assay. Hum. Mutat.42 (9), 11531164. 10.1002/humu.24246

  • 50

    ZhuQ.RuiX.LiY.YouY.ShengX. L.LeiB. (2021). Identification of four novel variants and determination of genotype-phenotype correlations for ABCA4 variants associated with inherited retinal degenerations. Front. Cell Dev. Biol.9, 634843. 10.3389/fcell.2021.634843

Summary

Keywords

ABCA4 gene, inherited retinal dystrophy, whole-exome sequencing, minigene assay, Stargardt disease, retinitis pigmentosa

Citation

Luo Q, Huang J, Shi L, Zhang G, Xue L, Wu K, Li X, Yang L, Li D, Mao L and Luo J (2025) Identification and functional characterization of ABCA4 gene variants in three patients with Stargardt disease or retinitis pigmentosa. Front. Genet. 16:1516872. doi: 10.3389/fgene.2025.1516872

Received

25 October 2024

Accepted

13 May 2025

Published

18 June 2025

Volume

16 - 2025

Edited by

Akio Oishi, Nagasaki University, Japan

Reviewed by

Ivana Mihalek, University of Rijeka, Croatia

Zelia Corradi, Radboud University Medical Centre, Netherlands

Updates

Copyright

*Correspondence: Jihong Luo,

† These authors share first authorship

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics