Skip to main content

Clinical exome sequencing facilitates the understanding of genetic heterogeneity in Leber congenital amaurosis patients with variable phenotype in southern India

Abstract

Background

Leber congenital amaurosis (LCA), primarily characterized by retinal degeneration is the most severe form of inherited retinal dystrophy (IRD) responsible for congenital blindness. The presence of phenotypic heterogeneity makes the diagnosis of LCA challenging, especially in the absence of pronounced disease pathognomonic, yet it can be well comprehended by employing molecular diagnosis. Therefore, the present study aimed to reveal the causative mutations in ten LCA patients with variable phenotypes using clinical exome sequencing (CES).

Methods

CES was performed in ten unrelated LCA patients. Ophthalmic information and family history of all patients were obtained to make a meaningful interpretation. The clinical exome data was analyzed and prioritized using a bioinformatics pipeline to identify mutations, which was further validated by Sanger sequencing. Segregation analysis was also performed on available family members.

Results

CES led to the identification of causative mutations in nine LCA patients. Seven patients harbored a mutation in six LCA candidate genes, including RPE65, LCA5 (n = 2), CRX, PRPH2, CEP290, and ALMS1, while two patients possess a mutation in IFT80 and RP1, known to cause other diseases. Three novel mutations in LCA5 (c.1823del), CRX (c.848del) and CEP290 (c.2483G > T) were identified. The current study reports for the first time, a mutation in PRPH2, CEP290, and ALMS1 from the Indian population. Additionally, we observed a novel association of LCA phenotype with IFT80 known to cause Jeune syndrome. Based on the genetic finding, the patient AS09, who harbored a mutation in the RP1 gene, was re-diagnosed with early-onset retinitis pigmentosa.

Conclusion

In conclusion, the results underline the importance of CES in clinically diagnosed LCA patients with variable phenotypes. The correlation between mutations in candidate genes and clinical phenotypes, helps to refine the clinical diagnosis. However, molecular evaluation with a larger cohort of LCA patients is needed for better understanding of the mutational spectrum in southern India.

Background

Inherited retinal dystrophies (IRD) are a heterogeneous group of diseases that cause significant vision loss due to irreversible retinal degeneration. Leber congenital amaurosis (LCA) is one of the most severe and earliest forms of IRD responsible for infantile blindness with an estimated prevalence of 3 in 100,000 worldwide [1]. The incidence of LCA in the South Indian population is quite often due to consanguineous marriages and genetically isolated communities [2]. The disease begins in the first year of life with photoreceptors degeneration (rod and cone cells) and progresses through serious visual defects. The clinical hallmark of LCA includes decreased visual acuity, non-recordable electroretinogram (ERG), and sluggish pupillary response. In addition, other commonly observed clinical signs are nystagmus, Franceschetti’s oculo-digital sign, strabismus, high hyperopia, cataract, and keratoconus [3].

Although LCA is a monogenic disease, mutations in more than 29 genes have been implicated. Among these, twenty-six genes follow the autosomal recessive pattern, the classical mode of inheritance in LCA. Two genes IMPDH1 and OTX2, inherit the disease in an autosomal dominant manner, while the CRX gene is inherited either in an autosomal dominant or recessive pattern [3, 4]. Approximately 70% of these genes contribute to non-syndromic LCA cases. Previous studies from the Indian population have observed mutations in fourteen candidate genes, including GUCY2D, RPE65, AIPL1, RPGRIP1, LCA5, IQCB1, CRB1, SPATA7, RDH12, NMNAT1, KCNJ13, CRX, RD3, and TULP1 [5,6,7]. Few LCA-associated genes like CEP290, ALMS1, IFT140, and IQCB1 also contribute to other syndromes such as Joubert syndrome, peroxisomal disease, Alstrom syndrome, Batten disease, and Senior Loken syndrome with similar ocular manifestations as observed in LCA. It mostly overlaps with a milder form of the same disease called severe early childhood-onset retinal dystrophy or early-onset retinitis pigmentosa [3, 8]. Hence, the existence of various clinical phenotypes necessitates molecular genetic testing to identify causative mutations for accurate diagnosis at the earliest instance.

Previous studies from the Indian cohort have screened very few LCA candidate genes through Sanger sequencing, homozygosity mapping, micro-array, and disease-specific targeted sequencing [5,6,7]. Hence, the present study employs clinical exome sequencing (CES), which targets ~ 8000 genes with known clinical implications, and thus bypass the shortcoming of other previous techniques. To our knowledge, this is the first study to use the CES approach for the molecular diagnosis of LCA patients in South India.

Therefore, the current study aimed to identify the underlying disease mutation in ten LCA patients with variable phenotypes using CES. Through the comprehensive analysis of clinical and genetic datasets, this study will contribute to the existing knowledge of genotype-phenotype associations towards LCA.

Methods

Ethics statement

This study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Ethics Committee of Aravind Eye Hospital, Madurai, Tamil Nadu, India (IRB2016017BAS). Written informed consent was obtained from all study patients or guardians in the case of minors or children.

Patient recruitment

Ten unrelated individuals (AS01 – AS10) diagnosed with LCA were recruited from the Paediatric Clinic, Aravind Eye Hospital, Madurai, Tamil Nadu, India. All the patients are of South Indian origin (Tamil Nadu, Kerala, and Andhra Pradesh). Comprehensive ophthalmic examinations including visual acuity, cycloplegic refraction, color fundus photography (Topcon, Inc., Tokyo, Japan), spectral-domain optical coherence tomography (SD-OCT) and autofluorescence were performed for all study patients. ERG was recorded through the UTAS Ganzfeld-LKC technology system and Burien-Allen bipolar electrodes based on the standards of the International Society for Clinical Electrophysiology of Vision. In children below 6 years of age, ERG was performed under ketamine anesthesia.

Clinical diagnosis of LCA was based on the following criteria: i) severe visual impairment during the first year of life, especially with Franceschetti’s oculo-digital signs (eye-poking, rubbing and pressing); ii) non-recordable ERG; iii) Nystagmus or roving eye movement. A detailed pedigree was obtained as well as other particulars such as ethnic predisposition, family history and consanguinity.

DNA isolation and CES

Genomic DNA was extracted from the peripheral blood of patients and available family members using the modified salt precipitation method [9]. CES of 10 LCA patients was performed at Medgenome, Bangalore, India. Sequencing libraries were prepared using clinical exome panel (Cev3), which covers approximately 8332 diseases causing genes, including 29 known LCA genes. Paired-end sequencing was performed to generate 2 × 150 bp reads at 100× sequencing depth using the HiSeq X Ten platform.

CES data analysis

The pre-processing of the fastq file was performed using Cutadapt (v1.8) to exclude low-quality reads, adapters and primer sequences [10]. The pre-processed reads were aligned against the human genome reference sequence hg19 (GRCh37) using Burrows-Wheeler Aligner (BWA)-MEM (v.0.7.12) [11]. Picard tool (v.1.140; https://broadinstitute.github.io/picard/) was employed to remove PCR-duplicates. The IndelRealigner and Base Recalibrator from the Genome Analysis Toolkit (GATK, v.3.6) were used for local realignment in regions containing potential indels and recalibrating the base quality scores of all reads. Both the GATK Haplotypecaller and UnfiedGenotyper were used for variant calling [12]. Then, the variants were annotated by the Variation and Mutation Annotation Toolkit (VariMAT, v.2.4.1; https://omictools.com/varimat-tool).

Variant prioritization

Pathogenic variants were prioritized as per the ACMG (American College of Medical Genetics and Genomics) standards and guidelines [13]. Briefly, Nonsense, frameshift, canonical ±1 or 2 splice sites, non-synonymous and in-frame variants located in the exonic region were considered for variant prioritization. The annotated variants were screened against databases such as Exome Aggregation Consortium (ExAC) and 1000 Genomes to exclude variants with allele frequency > 0.01.

The non-synonymous variants were considered as pathogenic only when at least four out of five insilico functional prediction algorithms such as SIFT (https://sift.bii.a-star.edu.sg/), PolyPhen2 (http://genetics.bwh.harvard.edu/pph2/), Mutation Taster (http://www.mutationtaster.org/), Mutation Assessor (http://mutationassessor.org/r3/) and FATHMM (http://fathmm.biocompute.org.uk/) were predicted to be deleterious. The conservation tools including GERP, SiPhy and PhastCons were used to predict the impact of non-synonymous variants. HOPE (Have (y) Our Protein Explained) was used to predict the structural impact of non-synonymous variants [14].

Validation of variants and segregation analysis

The identified pathogenic variants were validated in probands using Sanger sequencing. Forward and reverse primers for Sanger sequencing was designed using Primer-BLAST (Supplementary Table S1). The chromatograms were visualized using chromas v.2.6.6 software (http://www.technelysium.com.au/chromas.html) and the nucleotide sequence was analyzed using BLAST. The pathogenic variants were also investigated in ethnically matched control samples. Segregation analysis of variants was also performed in available family members.

Results

Clinical characteristics

Nystagmus or roving eye movement and the oculo-digital sign was the consistent feature observed in all patients. Based on the classification of visual acuity according to the World Health Organization’s ICD-11 (International Classification of Disease 11) 2018, seven patients were legally blind while three had severe visual impairment. Scotopic and photopic responses of ERG were non-recordable in all patients. The fundus, SD-OCT, and autofluorescence of all patients are shown in Figs. 1 and 2. Two patients (AS06 and AS10) had systemic features such as head nodding, secondary behavioral changes, kidney cyst, and ichthyosis. Table 1 summarizes the clinical findings of each patient with a respective genotype.

Fig. 1
figure1

Fundus and SD-OCT of study patients. Fundus presentation ranges from greyish desaturated background to pigmentary retinopathy. d, g, i and j Patients AS03, AS06, AS08 II:2 and AS08 II:3 affected by LCA5, CEP290 and IFT80 were noted with marbled fundus. i and j AS08 II:2 and AS08 II:3 also had macular coloboma indicated by an arrow. f Yellow vitelliform, egg yolk like well-circumscribed lesion centered at the fovea less than 1/3 of disc with central hyperpigmented spot was observed in patient AS05 with PRPH2 mutation indicated by an arrow. k white arrow indicates the Bull’s eye macula in patient AS09 carrying the RP1 mutation. The accompanying SD-OCT revealed normal retinal architecture in patients with RPE65, LCA5 (AS02) and PRPH2, whereas other patients had a lack of lamination or distorted retina of variable thickness resembling an immature retina. bd and h Patients AS01, AS02, AS03 and AS07 affected by RPE65, LCA5 and ALMS1 mutations showed preserved outer retinal layer outlined by white lines. f White upper arrow specifies the vitelliform lesions at the macula in patient AS05 affected by PRPH2 mutation. i and j AS08 II:2 and the twin AS08 II:3 presented with crater-like depression indicated by a white down arrow. Further information on these patients are described in Table 1

Fig. 2
figure2

Fundus Autofluorescence photographs of some patients. Autofluorescence imaging was performed in five patients. Deviation from normal was noted in all patients

Table 1 Clinical characteristics of ten South Indian patients involved in this study

Pedigree

Except for patient AS04, all other patients had suspected autosomal recessive pedigree due to the consanguinity and absence of consecutive generation disease. Patient AS04 may have an autosomal dominant pattern as the patient’s mother had a history of impaired vision while the father and younger sibling are normally sighted. The pedigree of all patients is shown in Fig. 3.

Fig. 3
figure3

Pedigree of ten unrelated patients involved in the study. Solid symbols with an arrow – probands; Solid symbols (M/M) – homozygous affected; half-filled symbols (M/+) – heterozygous affected; dotted symbols (M/+) – carrier; unfilled symbol (+/+) – Wild type; Slash through the symbol – deceased, Consanguinity represented by the double line

Data analysis

On average, 7.5 GB of data were generated per exome, of which 95% is above Q30. Each exome contains approximately 49,000,000 reads, while an average of 47,000,000 reads remained after adapter trimming. Overall, 99.97% of reads aligned with human reference genome hg19, and 94% passed alignment. Among these mapped reads, around 86% of reads located in the target region with a sequence depth ranged from 99 to 139.09X.

Variant calling attained a total of 31,000 variants, including 30,700 substitutions and 300 indels per sample. Almost 17,000 variants were filtered by excluding variants in noncoding regions. Except for synonymous and UTR (untranslated region), nearly 6000 variants were considered for further analysis. We retained around 600 variants (100 indels and 500 substitutions) with MAF ≤ 0.01. LCA can inherit by either autosomal recessive or dominant patterns, so both homozygous and heterozygous variants were considered for analysis. Among the 500 substitutions, 50 variants were predicted to be potentially pathogenic by functional prediction and conservation tools. From the 150 variants (100 indels and 50 substitutions), putative pathogenic variants were identified. The overall summary of exome data and variant prioritization are shown in Supplementary Tables S2 and S3.

Pathogenic mutations and segregation analysis

CES identified disease-causing mutations in nine (AS01 to AS09) out of ten patients, of which seven (AS01 to AS07) patients harbored mutations in the LCA-associated genes, including RPE65, LCA5, CRX, PRPH2, CEP290, and ALMS1. Among these, three mutations (AS03 - LCA5: c.1823del, p.Leu608TyrfsTer30; AS04 - CRX: c.848del, p.Met283ArgfsTer88; AS06 - CEP290: c.2483G > T, p.Ser828Ile) were novel. To our knowledge, this is the first report of a mutation of these LCA candidate genes, including PRPH2: c.629C > T, p.Pro210Leu; CEP290: c.2483G > T, p.Ser828Ile; and ALMS1: c.11310_11313del, p.Glu3771TrpfsTer18 from an Indian population. The remaining two patients AS08 and AS09 were identified with other retinal disease genes IFT80 and RP1. The mutation (AS08 - c.1936G > T, p.Val646Phe) in the IFT80 gene has not yet been reported to be associated with the LCA phenotype. Molecular diagnosis of AS09 with RP1 mutation (c.3751_3752del, p.Val1251PhefsTer9) led to a revision of the clinical diagnosis as early-onset retinitis pigmentosa. The details of the mutations identified in this study have been summarized in Table 2.

Table 2 Putative variants of nine South Indian patients identified by clinical exome sequencing

As per the ACMG guideline, six of them (AS01 to AS04, AS07, and AS09 in Table 2) were classified to have pathogenic mutations based on the combination of the following criteria’s: PVS1: Null variant in a gene where the loss of function is a known mechanism of disease, PM2: Absent or low frequency in population databases, PP3: Several computational evidence for the mutation’s deleterious effect, PP5: Reported as pathogenic by a reputable source and PP4: Patients phenotype or family history supports variant. Patient AS05 carried a likely pathogenic mutation, which follows PM1: Variant at hotspots or functional domains, PM2, PP3, and PP5. The other two (AS06 and AS08) had variants of uncertain significance.

The identified mutations were validated in patient samples using Sanger sequencing. Chromatograms of the novel and reported mutations are shown in Fig. 4 and Supplementary Figure S1.

Fig. 4
figure4

Sequence chromatogram and segregation analysis of mutations identified by this study. a Autosomal recessive pattern – Chromatogram shows the homozygous peak in the proband and heterozygous peak in carriers (parents and siblings). The affected sibling of AS06 (II:3) has a homozygous peak; b Autosomal dominant pattern of inheritance – Heterozygous peak in affected proband and mother, wild-type peak in control and other family members (father and sibling). The red box and black arrow indicate the altered nucleic acids

Discussion

The existence of high clinical heterogeneity ensuing from intricate genetics has been demonstrated for LCA etiology. Hence, CES was performed for ten LCA patients with variable phenotypes to unravel causative mutations contributing to disease pathogenesis.

RPE65 is one of the most common LCA candidate genes and mutation in this gene contributes to 3–16% of LCA cases worldwide [4]. Patient AS01 was identified with a pathogenic frameshift mutation in the carotenoid oxygenase domain of the RPE65 gene, which creates a stop codon at exon 5, and thus produces a truncated protein of 129 amino acids residues. It might have partial function compared to the wild-type protein of 533 amino acids or may undergo nonsense-mediated mRNA decay. This category of gene mutations causes a deficiency of 11-cis-retinal especially in rods compared to cones, leading to nyctalopia [3, 8]. Patients with RPE65 mutation also exhibit some extent of improvement in visual acuity over the first decade of life, but it will eventually deteriorate in the later decades [3]. Our patient also had a history of nyctalopia and showed transient improvement in visual acuity, consistent with earlier reports. The same mutation was formerly identified in South Indian LCA patients among many other mutations identified at this c.361 codon position, which makes it a hot spot for Indian LCA patients [5, 7, 15].

LCA5 consists of 697 amino acids, which encodes highly conserved ciliary protein Lebercilin. Despite its wide expression in human tissues, LCA5 mutations are restricted to only cause retinal dysfunction with a prevalence rate of 1–2% [4]. The majority of mutations reported with the LCA5 phenotype are null mutations [16]. Our study identified different homozygous null mutations in two unrelated patients AS02 [17, 18] and AS03. Patient AS02 below 10 years of age showed mild improvement in vision, while patient AS03 was legally blind at 11 years of age. It was reported that patients with LCA5 mutations showed improvement in vision and eventually decline after the first decade of life [19]. Even though Patient AS02 showed normal macula in fundoscopy and preserved central macular outer retinal layer in optic coherence tomography, increased autofluorescence at the macula was observed, suggesting increased metabolic activity of the retinal pigment epithelium. A marbled fundus was observed in patient AS03, as seen in the CEP290-related LCA phenotype. Patient AS03 developed a posterior subcapsular cataract, reported as a common feature of patients with LCA5 mutations [16].

CRX (cone-rod homeobox) is often reported to cause LCA in an autosomal dominant pattern [20]. The present study also provides evidence by identifying a novel heterozygous frameshift mutation in patient AS04. Segregation analysis revealed the same mutation in the affected mother, while absent in normally sighted father and younger sibling, indicating that the patient inherited a disease-related mutation from his mother. CRX consists of 299 amino acid residues. Whereas in the proband AS04, frameshift shift mutation in exon 4 removes the native stop codon, as a result of which a larger open reading frame consisting of 369 amino acids is produced. Earlier reports on CRX mutations in the index cases have shown thinned, abnormal lamellar structure, and macular atrophy without noticeable signal of inner and outer segments junction in SD-OCT [21]. Similarly, patient AS04 had a lack of lamination and complete loss of the outer retinal layer.

Individuals affected by PRPH2 mutation are known to have pattern dystrophy (butterfly-shaped pigment dystrophy and Adult-onset foveomacular vitelliform dystrophy) with a broad spectrum of clinical appearance, LCA and retinitis pigmentosa [22, 23]. Patient AS05 carried a reported missense mutation in the cytoplasmic domain of PRPH2 [24]. HOPE predicted that the wild-type residue proline is very rigid and required special conformation of the protein backbone. Thus, the mutation at that position might affect protein function by disrupting special conformation [12]. The vitelliform lesion is the most commonly encountered clinical presentation of adult-onset foveomacular vitelliform dystrophy and LCA, due to PRPH2-mediated phenotypes [22, 23]. In this study, patient AS05 also demonstrated early retinal defects with vitelliform lesions, consistent with the abovementioned studies.

Interestingly, CES analysis of patient AS06 revealed a novel homozygous missense mutation in the CEP290 gene. The HOPE tool has predicted that the mutation is present in the conserved region required for interaction with IQCB1 (IQ Motif Containing B1). This mutation introduces more hydrophobic residues, which may affect the hydrogen-bond formation and results in loss of interactions with other molecules [12]. Mutation in CEP290 leads to LCA, Bardet–Biedl syndrome, Senior–Loken syndrome and Joubert syndrome [25]. Unlike the Western population, where CEP290 mutations are the most common cause for LCA with an estimated prevalence of 15–20% [4, 25], to date no CEP290 mutations have been reported from the Indian population [13, 26]. At the age of 11 years, patient AS06 developed a cataract and presented with marbled fundus as reported in earlier studies [27, 28]. However, the presence of systemic anomalies like head nodding and secondary behavioral changes in the patient indicates an association with the abovementioned syndromes, therefore AS06 was advised to undergo complete systemic evaluation including an MRI scan and ultrasound to better understand the disease pathogenesis in this patient. However, due to kidney failure, the patient was deceased. The patient’s younger sibling (2 years) also presented with similar ocular features, head nodding and secondary behavioral changes. Since the sibling was identified with the same mutation, being strictly monitored, we hope to identify related syndromes and offer treatment early.

Mutation in ALMS1 was reported to cause Alstrom syndrome and LCA [29, 30]. One of our patients AS07 carried a homozygous frameshift mutation in ALMS1. Indeed, so far ALMS1 mutation-specific phenotypes have not been described extensively for LCA patients. However, Xu and colleagues have reported homozygous ALMS1 null mutation in six LCA cases with early-onset retinal degeneration, visual acuity from light perception to no light, high hyperopia, roving eye movement, oculo-digital sign, undetectable ERG and tapetal fundus [30]. In agreement with this study, our patient also developed these clinical features, but the fundus revealed a pink disc, thread-like arteriolar attenuation, greyish desaturated appearance and wrinkling of the inner limiting membrane at the macula. The patient did not demonstrate any typical signs of Alstrom syndrome until now, and thus was advised to return for regular clinical assessments for any systemic abnormalities.

IFT80 is a component of Intra Flagellar Transport complex B, which is essential for the assembly and maintenance of motile and sensory cilia [31]. IFT80 mutants underlie Jeune syndrome, an autosomal recessive disease characterized by the constricted thoracic cage, respiratory insufficiency, cystic renal disease, polydactyl disease and retinal degeneration. Patients with Jeune syndrome are likely to develop retinal dystrophies within a few months of birth [32, 33]. In the present work, one of the twins AS08 subjected to CES was identified with a missense mutation in the IFT80 gene. Later, the same mutation was found in another twin through segregation analysis. Both twins were presented with typical LCA clinical features as poor vision at an early stage, multiplanar nystagmus and non-recordable ERG. At 11 years, the patient observed with marbled fundus, similar to CEP290-associated LCA patients. The patient has not yet developed any other symptoms of Jeune syndrome except for the LCA phenotype. Many ciliopathy genes, such as CEP290, IFT140, and IQCB1 initially reported to cause syndromic retinal degeneration, but were later identified as having vital roles in LCA pathogenesis as well [25, 34, 35]. Similarly, the IFT80 gene might also be contributing to the LCA phenotype.

“The patient AS09 was initially diagnosed with LCA as she presented with poor visual acuity, multiplanar nystagmus, oculodigital sign, and non-recordable ERG at the age of 1.75 years. Interestingly, molecular diagnosis identified a homozygous frameshift mutation in exon 4 of the RP1 gene, documented to cause retinitis pigmentosa. Previous studies have shown that the mutation in exon 4 of RP1 causes retinitis pigmentosa in early life, as it encodes 85% of the protein [36]. Similarly, the null mutations in exon 4 of RP1 might cause retinitis pigmentosa in the early stages of life. Thus, based on the genetic finding, the patient was re-defined with clinical diagnosis of autosomal recessive severe early-onset retinitis pigmentosa.”

The patient AS10, who was clinically diagnosed with LCA, was not identified with any underlying variant. In addition, the patient was also noted with dry, scaly skin and hyperpigmented knuckles. Whole exome or genome sequencing is required to understand the genetic cause of this disease phenotype.

Conclusion

CES delineates the causal mutations of the LCA patients involved in this study. Our study established the involvement of new genes in LCA pathogenesis. The molecular finding also helps to assess the risk stratification due to LCA-associated syndromes. Furthermore, our results provide insight into the phenotypic features of study patients, which aid in more accurate clinical diagnosis. However, to obtain a complete picture of the disease, a large cohort of LCA patients must be evaluated and correlated with the corresponding phenotype.

Availability of data and materials

All data generated during this study are included in this article and its additional files.

Abbreviations

LCA:

Leber Congenital Amaurosis

IRD:

Inherited Retinal Dystrophy

CES:

Clinical Exome Sequencing

ERG:

Electroretinogram

SD-OCT:

Spectral-domain optical coherence tomography

References

  1. 1.

    Koenekoop RK. An overview of Leber congenital amaurosis: a model to understand human retinal development. Surv Ophthalmol. 2004;49(4):379–98.

  2. 2.

    Kumaramanickavel G, Joseph B, Vidhya A, Arokiasamy T, Shridhara Shetty N. Consanguinity and ocular genetic diseases in South India: analysis of a five-year study. Community Genet. 2002;5(3):182–5. 

  3. 3.

    Chacon-Camacho OF, Zenteno JC. Review and update on the molecular basis of Leber congenital amaurosis. World J Clin Cases. 2015;3(2):112–24. 

  4. 4.

    Coussa RG, Lopez Solache I, Koenekoop RK. Leber congenital amaurosis, from darkness to light: an ode to Irene Maumenee. Ophthalmic Genet. 2017;38(1):7–15.

  5. 5.

    Verma A, Perumalsamy V, Shetty S, Kulm M, Sundaresan P. Mutational screening of LCA genes emphasizing RPE65 in south Indian cohort of patients. PLoS One. 2013;8(9):e73172.

  6. 6.

    Srilekha S, Arokiasamy T, Srikrupa NN, Umashankar V, Meenakshi S, Sen P, et al. Homozygosity mapping in Leber congenital amaurosis and autosomal recessive retinitis pigmentosa in south Indian families. PLoS One. 2015;10(7):e0131679. 

  7. 7.

    Srikrupa NN, Srilekha S, Sen P, Arokiasamy T, Meenakshi S, Bhende M, et al. Genetic profile and mutation spectrum of Leber congenital amaurosis in a larger Indian cohort using high throughput targeted re-sequencing. Clin Genet. 2018;93(2):329–39.

  8. 8.

    Kumaran N, Moore AT, Weleber RG, Michaelides M. Leber congenital amaurosis/early-onset severe retinal dystrophy: clinical features, molecular genetics and therapeutic interventions. Br J Ophthalmol. 2017;101(9):1147–54. 

  9. 9.

    Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988;16(3):1215.

  10. 10.

    Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10. https://doi.org/10.14806/ej.17.1.200.

    Article  Google Scholar 

  11. 11.

    Li H, Durbin R. Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics. 2010;26(5):589–95.

  12. 12.

    McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.

  13. 13.

    Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17(5):405–24.

  14. 14.

    Venselaar H, Te Beek TA, Kuipers RK, Hekkelman ML, Vriend G. Protein structure analysis of mutations causing inheritable diseases. An e-Science approach with life scientist friendly interfaces. BMC Bioinformatics. 2010;11:548.

    Article  Google Scholar 

  15. 15.

    Glen WB Jr, Peterseim MMW, Badilla R, Znoyko I, Bourg A, Wilson R, et al. A high prevalence of biallelic RPE65 mutations in Costa Rican children with Leber congenital amaurosis and early-onset retinal dystrophy. Ophthalmic Genet. 2019;40(2):110–7. 

  16. 16.

    Mackay DS, Borman AD, Sui R, van den Born LI, Berson EL, Ocaka LA, et al. Screening of a large cohort of Leber congenital amaurosis and retinitis pigmentosa patients identifies novel LCA 5 mutations and new genotype–phenotype correlations. Hum Mutat. 2013;34(11):1537–46.

  17. 17.

    Beryozkin A, Zelinger L, Bandah-Rozenfeld D, Shevach E, Harel A, Storm T, et al. Identification of mutations causing inherited retinal degenerations in the Israeli and Palestinian populations using homozygosity mapping. Invest Ophthalmol Vis Sci. 2014;55(2):1149–60. 

  18. 18.

    Tajiguli A, Xu M, Fu Q, Yiming R, Wang K, Li Y, et al. Next-generation sequencing-based molecular diagnosis of 12 inherited retinal disease probands of Uyghur ethnicity. Sci Rep. 2016;6(1):21384. 

  19. 19.

    den Hollander AI, Roepman R, Koenekoop RK, Cremers FP. Leber congenital amaurosis: genes, proteins and disease mechanisms. Prog Retin Eye Res. 2008;27(4):391–419.

  20. 20.

    Arcot Sadagopan K, Battista R, Keep RB, Capasso JE, Levin AV. Autosomal-dominant Leber congenital amaurosis caused by a heterozygous CRX mutation in a father and son. Ophthalmic Genet. 2015;36(2):156–9.

  21. 21.

    Zou X, Yao F, Liang X, Xu F, Li H, Sui R, et al. De novo mutations in the cone-rod homeobox gene associated with Leber congenital amaurosis in Chinese patients. Ophthalmic Genet. 2015;36(1):21–6. 

  22. 22.

    Boon CJ, den Hollander AI, Hoyng CB, Cremers FP, Klevering BJ, Keunen JE. The spectrum of retinal dystrophies caused by mutations in the peripherin/RDS gene. Prog Retin Eye Res. 2008;27(2):213–35.

  23. 23.

    Khan AO, Al Rashaed S, Neuhaus C, Bergmann C, Bolz HJ. Peripherin mutations cause a distinct form of recessive Leber congenital amaurosis and dominant phenotypes in asymptomatic parents heterozygous for the mutation. Br J Ophthalmol. 2016;100(2):209–15.

  24. 24.

    Koyanagi Y, Akiyama M, Nishiguchi KM, Momozawa Y, Kamatani Y, Takata S, et al. Genetic characteristics of retinitis pigmentosa in 1204 Japanese patients. J Med Genet. 2019;56(10):662–70.

  25. 25.

    Coppieters F, Lefever S, Leroy BP, De Baere E. CEP290, a gene with many faces: mutation overview and presentation of CEP290base. Hum Mutat. 2010;31(10):1097–108. 

  26. 26.

    Sundaresan P, Vijayalakshmi P, Thompson S, Ko AC, Fingert JH, Stone EM. Mutations that are a common cause of Leber congenital amaurosis in northern America are rare in southern India. Mol Vis. 2009;15:1781–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Littink KW, Pott JW, Collin RW, Kroes HY, Verheij JB, Blokland EA, et al. A novel nonsense mutation in CEP290 induces exon skipping and leads to a relatively mild retinal phenotype. Invest Ophthalmol Vis Sci. 2010;51(7):3646–52.

  28. 28.

    Yzer S, Hollander AI, Lopez I, Pott JW, de Faber JT, Cremers FP, et al. Ocular and extra-ocular features of patients with Leber congenital amaurosis and mutations in CEP290. Mol Vis. 2012;18:412–25.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Marshall JD, Muller J, Collin GB, Milan G, Kingsmore SF, Dinwiddie D, et al. Alström syndrome: mutation spectrum of ALMS1. Hum Mutat. 2015;36(7):660–8.

  30. 30.

    Xu Y, Guan L, Xiao X, Zhang J, Li S, Jiang H, et al. ALMS1 null mutations: a common cause of Leber congenital amaurosis and early-onset severe cone–rod dystrophy. Clin Genet. 2016;89(4):442–7.

  31. 31.

    Yang S, Wang C. The intraflagellar transport protein IFT80 is required for cilia formation and osteogenesis. Bone. 2012;51(3):407–17.

  32. 32.

    Hudak LM, Lunt S, Chang CH, Winkler E, Flammer H, Lindsey M, et al. The intraflagellar transport protein ift80 is essential for photoreceptor survival in a zebrafish model of jeune asphyxiating thoracic dystrophy. Invest Ophthalmol Vis Sci. 2010;51(7):3792–9. 

  33. 33.

    Moran J, Sanderson KG, Maynes J, Vig A, Batmanabane V, Kannu P, et al. IFT80 mutations cause a novel complex ciliopathy phenotype with retinal degeneration. Clin Genet. 2018;94(3–4):368–72.

  34. 34.

    Xu M, Yang L, Wang F, Li H, Wang X, Wang W, et al. Mutations in human IFT140 cause non-syndromic retinal degeneration. Hum Genet. 2015;134(10):1069–78. 

  35. 35.

    Otto EA, Loeys B, Khanna H, Hellemans J, Sudbrak R, Fan S, et al. Nephrocystin-5, a ciliary IQ domain protein, is mutated in Senior-Loken syndrome and interacts with RPGR and calmodulin. Nat Genet. 2005;37(3):282–8. 

  36. 36.

    Silva RS, Salles MV, Motta FL, Sallum JM. Retinitis pigmentosa due to RP1 biallelic variants. Sci Rep. 2020;10(1):1603. 

Download references

Acknowledgements

We thank all the study participants for contributing their medical information and genetic material for this research. We thank Dr. Roopam Duvesh, Research Associate, Department of Molecular Genetics, Aravind Medical Research Foundation, Madurai for her support in editing the manuscript, as well as Ms. M. Gowri, Aravind Eye Care System, Madurai for her help in preparing the figures.

Funding

This work was supported by the Department of Biotechnology under Grant BT/NNT/28/SP18830/2018.

Author information

Affiliations

Authors

Contributions

SV performed the molecular genetic studies, clinical exome data analysis and prioritization, Sanger sequencing and drafted the manuscript. RA and VP recruited study subjects and performed clinical assessments. PS supervised the work and helped to draft and critically review the manuscript. AV and DB were involved in study design. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Periasamy Sundaresan.

Ethics declarations

Ethic approval and consent to participate

This study was approved by the Institutional Ethics Committee of Aravind Eye Hospital, Madurai, Tamil Nadu, India (IRB2016017BAS).

Consent for publication

Not applicable.

Competing interests

None of the authors have any proprietary interests or conflicts of interest related to this submission.

Supplementary Information

Additional file 1: Supplementary Table S1.

Primer sequences used for the mutation validation by Sanger sequencing.

Additional file 2: Supplementary Table S2.

Summary of clinical exome data.

Additional file 3: Supplementary Table S3.

Overview of exclusion and prioritization of variants to obtain pathogenic variant from clinical exome data.

Additional file 4: Supplementary Figure S1.

Sanger validation and Segregation analysis of reported mutations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Viswarubhiny, S., Anjanamurthy, R., Vanniarajan, A. et al. Clinical exome sequencing facilitates the understanding of genetic heterogeneity in Leber congenital amaurosis patients with variable phenotype in southern India. Eye and Vis 8, 20 (2021). https://doi.org/10.1186/s40662-021-00243-5

Download citation

Keywords

  • Leber congenital amaurosis
  • Clinical exome sequencing
  • Southern India
  • Molecular diagnosis
  • Genotype-phenotype correlation