Molecular Vision 2018; 24:667-678 <http://www.molvis.org/molvis/v24/667>
Received 11 June 2018 | Accepted 16 October 2018 | Published 18 October 2018

Non-homologous recombination between Alu and LINE-1 repeats results in a 91 kb deletion in MERTK causing severe retinitis pigmentosa

Frida Jonsson,1 Marie Burstedt,2 Therese G. Kellgren,3,4 Irina Golovleva1

1Medical Biosciences/Medical and Clinical Genetics, University of Umeå, Umeå, Sweden; 2Clinical Science/Ophthalmology, University of Umeå, Umeå, Sweden; 3Department of Mathematics and Mathematical Statistics, Umeå University, Umeå, Sweden; 4Computational Life Science Cluster (CLiC), Umeå University, Umeå, Sweden

Correspondence to: Irina Golovleva, Clinical Genetics, University Hospital of Umeå, 90185 Umeå, Sweden. ORCID 0000-0001-8741-0616. Phone: +467856820; FAX: +4690128163; email: irina.golovleva@umu.se

Abstract

Purpose: Retinitis pigmentosa (RP) represents a large group of inherited retinal diseases characterized by clinical and genetic heterogeneity. Among patients with RP in northern Sweden, we identified two severely affected siblings and aimed to reveal a genetic cause underlying their disease.

Methods: Whole exome sequencing (WES) was performed on both affected individuals. Sequence variants were filtered using a custom pipeline to find a rare or novel variant predicted to affect protein function. Genome-wide genotyping was used to identify copy number variants (CNVs) and homozygous regions with potential disease causative genes.

Results: WES uncovered a novel heterozygous variant in the MER proto-oncogene, tyrosine kinase (MERTK) gene, c.2309A>G, p.Glu770Gly located in the tyrosine kinase domain and predicted to be likely pathogenic. The second variant, a large heterozygous deletion encompassing exons 1 to 7 of the MERTK gene, was revealed with genome-wide genotyping. The CNV analysis suggested breakpoints of the deletion, in the 5′-untranslated region and in intron 7. We identified genomic sequences at the site of the deletion as part of L1ME4b (LINE/L1) and AluSx3 that indicated a non-homologous recombination as a mechanism of the deletion evolvement.

Conclusions: Patients with RP in this study were carriers of two novel allelic mutations in the MERTK gene, a missense variant in exon 17 and an approximate 91 kb genomic deletion. Mapping of the deletion breakpoints allowed molecular testing of a cohort of patients with RP with allele-specific PCR. These findings provide additional information about mutations in MERTK for molecular testing of unsolved recessive RP cases and highlight the necessity for analysis of large genomic deletions.

Introduction

Inherited retinal diseases (IRDs) represent genetic conditions leading to visual deficiency and in some instances, to blindness [1]. The most common among IRDs is retinitis pigmentosa (RP). RP is a progressive non-syndromic rod-cone disease with clinical and genetic heterogeneity [2]. The genetic heterogeneity and founder mutations make some genes more common in isolated geographic areas. The population of northern Sweden has been previously described as ideal for genetic research because the country’s rivers and landscape features made it difficult for the inhabitants to move around hundreds of years ago, contributing to a homogenous population [3,4]. A genetic cause of recessive RP in this region is known in less than 50% of the patients; thus, the search for pathogenic disease causative variants continues.

There are at least 60 genes and three loci associated with only autosomal recessive RP (arRP) according to the Retinal Information Network. Many mutations are private or family-specific, responsible for a very small percentage of arRP cases. One of the RP genes is MER proto-oncogene, tyrosine kinase (MERTK; Gene ID ENS00000153208, OMIM 604705), which was mutated in patients with RP from consanguineous families originating in the Middle East, Saudi Arabia, Spain, and Morocco [5-8]. In general, mutations in MERTK are rare affecting less than 1% of patients with RP [6-9] although in a large French cohort the mutations account for 2% [10], and in the isolated population of the Faroe Islands, a large MERTK deletion is responsible for 30% of the RP cases [11].

In this study, we describe two siblings with a severe form of recessive RP who carried a novel missense variant in exon 17 of the MERTK gene on one allele and a novel large genomic deletion encompassing exons 1 to 7 on the second allele. We also consider a non-homologous recombination between genomic repetitive sequences such as LINE1/L1 and Alu as a plausible molecular mechanism in the appearance of large genomic deletions in the MERTK gene.

Methods

Patients and clinical examination

The now elderly siblings from a family originating from Jämtland County in northern Sweden (Figure 1) were interviewed and examined, and their blood was collected at the Östersund Eye Clinic in the early 1990s. After recent genetic findings, their medical records were reanalyzed at the Eye Clinic of the University Hospital of Umeå. The Regional Ethical Review Board in Umeå approved this research project, and the study adhering to the ARVO statement on human subjects was performed according to the Declaration of Helsinki with informed consent obtained from both patients.

DNA and RNA were available from both affected individuals, and DNA of 100 healthy individuals including blood donors and men drafted at 18 years of age from a matched geographic population was included in the study. A cohort of patients with IRDs (n=145) without known genetic cause of their visual impairment with presumed arRP was also included in this study.

High-resolution SNP array

Genomic DNA was isolated from peripheral venous blood as previously described [12]. In short, 10 ml peripheral blood was obtained by venipuncture and collected in EDTA tubes. Genomic DNA was extracted from blood leukocytes using salting out procedure [13]. Genome-wide genotyping was performed using the HumanOmniExpress-24 BeadChip with more than 715,000 single nucleotide polymorphisms (SNPs), in accordance with the manufacturer’s instructions (Illumina Inc., San Diego, CA). To identify regions of homozygosity (ROH), both siblings were genotyped. The data were analyzed using GenomeStudio software (Illumina), and cnvPartition 3.2.0 was applied for copy number variant (CNV) detection by retrieving the log R ratio (LRR, the ratio between the observed and expected probe intensities) and the B allele frequency (BAF). When a CNV is absent, the LRR is around zero, and the BAF is 0, 0.5, or 1 depending on genotypes AA, AB, and BB. Human genome GRCh38 (hg38) was used to assign all chromosome positions.

Molecular genetic analyses

As the first step, DNA from one patient (RP115) was analyzed in 2007 using the arrayed primer extension assay (APEX) available at Asper Biotech (Tartu, Estonia). The testing for arRP included 518 known variants in 12 genes (PDE6A – Gene ID ENSG00000132915, OMIM 180071; PDE6B – Gene ID ENSG00000133256, OMIM 180072; PNR/NR2E3 – Gene ID ENSG00000278570, OMIM 604485; RDH12 – Gene ID ENSG00000139988, OMIM 608830; RGR – Gene ID ENSG00000148604, OMIM 600342; RLBP1 – Gene ID ENSG00000140522, OMIM 180090; SAG – Gene ID ENSG00000130561, OMIM 181031; TULP1 – Gene ID ENSG00000112041, OMIM 602280; CRB1 – Gene ID ENSG00000134376, OMIM 604210; RPE65 – Gene ID ENSG00000116745, OMIM 180069; USH2A – Gene ID ENSG00000042781, OMIM 608400; and USH3A/CLRN1 – Gene ID ENSG00000163646, OMIM 606397).

Whole exome sequencing (WES) was performed in both patients using the exome enrichment kit (Nimblegen SeqCap EZ Exome Library SR) and the HiSeq2000 sequencer (Illumina) using the Genomic Services at Ambry Genetics (San Diego, CA). Generated paired-end reads were aligned to the human reference genome UCSC hg19 using the Illumina CASAVA 1.8.2 software. RTA 1.12.4 (HiSeq Control Software 1.4.5) was used for the initial data processing and base calling. The sequence quality filtering script was executed in the Illumina CASAVA software (version 1.8.2; Illumina, Hayward, CA). A single BAM alignment file was used in SoftGenetics’NextGENe v2.16 for SNP and indel analysis [14]. Additionally, WES data from a control group of six individuals from the same restricted geographic area of northern Sweden (West Bothnia) were available at Computational Life Science Cluster (CLiC), Umeå University. All control individuals were clinically examined, and no signs of any retinal pathology were demonstrated.

PCR on genomic DNA was performed as described elsewhere with gene-specific primers (Table 1). The temperature profile for the PCR consisted of initial denaturing at 95 °C for 10 min, followed by seven cycles of denaturing at 95 °C for 30 s, annealing at 64 °C (−1 °C/cycle) and extension at 72 °C for 30 s; 40 cycles of denaturing at 95 °C for 30 s, annealing at 58 °C for 30 s, extension at 72 °C for 30 s, and a final extension step of 72 °C for 10 min. PCR was performed with the AmpliTaq Gold system (Applied Biosystems, Foster City, CA) with 0.2 mM dNTPs, 1X buffer, 1.5 mM MgCl2, 0.2 mM of each primer, dimethyl sulfoxide (DMSO) 0.02%, 50 ng template, and 0.75 U polymerase in a total volume of 25 µl. The amplified PCR products were purified using Exonuclease I enzyme (Thermo Fisher scientific, Waltham, MA) before sequencing. Cycle sequencing PCR was performed using the BigDye Terminator kit, version 3.1 as recommended by the manufacturer (Applied Biosystems). The products of the sequencing reactions were analyzed on the ABI 3500xL Dx Genetic Analyzer (Applied Biosystems). Sequences were aligned and evaluated using the Sequencher software version 4.9 (Gene Codes Corporation, Ann Arbor, MI).

Allele-specific PCR followed by fragment analysis was performed as previously described with 0.4 mM of fluorescently labeled forward primer, 0.16 mM of a primer specific for wild-type allele, and 0.24 mM of primer specific for mutant alleles (Table 1). The PCR program consisted of denaturing at 94 °C for 10 min, 28 cycles of denaturing at 94 °C for 30 s, annealing at 60 °C for 30 s, extension at 72 °C for 30 s, and final extension steps at 72 °C for 10 min and at 60 °C for 45 min. The fragment analysis based on size measurement was done with capillary electrophoresis performed in the mix of 9.5 μl Hi-DiTM Formamide 3500 Dx series, 0.5 µl of size standard 600, and 1 µl PCR product injected in the ABI 3500xL (Applied Biosystems). The results were visualized by using GeneMapper Software 5.

For reverse transcriptase PCR (RT–PCR), RNA was extracted from whole blood using TRIzol reagent (Life Technologies, Carlsbad, CA) and reverse transcribed into cDNA using a High Capacity cDNA Reverse Transcription kit (Applied Biosystems) according to the manufacturer’s recommendation. RT–PCR was performed with two different enzymes. For amplicon using primers 15F and 18R (Table 1), the conditions and program was 0.2 mM dNTPs, 1X buffer, 1.5 mM MgCl2, 0.2 mM of each primer, 5 µl cDNA and 0.75 U AmpliTaq Gold polymerase in a total volume of 25 µl. The PCR program consisted of denaturing at 95 °C for 10 min, 39 cycles of denaturing at 95 °C for 30s, annealing at 57 °C, and extension at 72 °C for 30 s with a final extension step of 72 °C for 10 min. For amplicons with 7F and 18R primers (Table 1) a mix of 0.4 mM dNTPs, 1X buffer, 1.5 mM MgCl2, 0.2 mM of each primer, 5 µl cDNA, and 0.5 U Platinum Taq DNA polymerase (Invitrogen, Carlsbad, CA) in a total volume of 25 µl was used. The PCR program consisted of denaturing at 95 °C for 10 min, 39 cycles of denaturing 95 °C for 30 s, annealing at 60 °C, and extension at 72 °C for 1 min 20 s, and a final extension step at 72 °C for 7 min. Sequencing of these amplicons was performed as described previously.

In silico analyses of sequence variants were done using computational tools available via the Alamut software version 2.0 (Interactive Biosoftware, Rouen, France) that integrates Sorting Intolerant From Tolerant (SIFT), PolyPhen, Align GVGD, and MutationTaster as pathogenicity prediction tools for missense variants and allows assessment of the variants’ potential impact on splicing. All identified variants were described according to the Sequence Variant Nomenclature recommended by the Human Genome Variation Society (HGVS) and classified according to guidelines approved by the HGVS, the Human Variome Project (HVP), and the Human Genome Organization (HUGO) [15,16]. The variants were interpreted according to 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 by Richards et al. [17].

Results

Clinical findings

Both siblings presented a severe phenotype with night blindness from early childhood and affected visual acuity (VA) at an early school age, 5 and 8 years old (Table 2). The younger brother experienced severe glare problems in his teens. Both patients developed a severe visual handicap with low VA, large scotomas with only small remnants of the visual field in their teens and early adulthood leading to blindness and the need for special training for the visual handicap. Fundus changes and pigmentations were present in the posterior pole when the siblings were in their teens. Dense cataract and nystagmus were observed in RP116 in his 50s. Unfortunately, no fundus images were available from earlier medical records, and no images could be taken due to severe cataract in both patients at their most recent visit to the Eye Clinic. No other eye diseases were present. The absence of retinal disease in the siblings’ parents and other siblings strongly suggested recessive mode of inheritance.

Molecular genetic findings

Genetic testing for known mutations-- The DNA from one affected individual (RP115, Figure 1) was analyzed using the arrayed primer extension (APEX) at Asper Biotech (Tartu, Estonia) aiming to identify reported sequence variants associated with arRP. No disease-causing variants were found among the interrogated 518 variants in 12 genes.

Whole exome sequencing-- DNA from both affected individuals (RP115 and RP116) was subjected to WES along with DNA from six controls from the same geographic region. Coverage for all tested samples ranged from 96.08% to 96.29% with a depth of at least 20X. The apparent autosomal recessive inheritance pattern in the family suggested two strategies for data analysis aiming at identifying either biallelic homozygous or compound heterozygous pathogenic variants. The following steps were included in the analysis of the WES data: (1) identify homozygous or heterozygous variants present in both siblings, (2) exclude variants present in at least one of the control DNA, (3) exclude all synonymous variants and variants located more than 10 bp from exon–intron junctions in all transcripts (MapView), and (4) manually analyze all variants focusing on retinal genes from the RetNet Information Network Database (the number of genes and loci was 305). The data filtering resulted in identification of shared 85 homozygous (80 genes) and 885 (724 genes) heterozygous variants out of 39,119 called variants for RP115 and 39,560 for RP116 (Appendix 1, Appendix 2, and Appendix 3). Only four homozygous variants in four genes were present in the genes associated with IRDs according to RetNet. All of them were excluded as a potential disease cause because their population frequency was 10% or higher (Table 3).

Of the 885 heterozygous variants in 724 genes, 19 were found in 17 genes associated with IRDs according to RetNet (Table 4). The only gene harboring more than one variant was the LRP5 (Gene ID ENSG00000162337, OMIM 603506) associated with aberrant ocular vascularization and loss of vision in genetic disorders such as osteoporosis-pseudoglioma syndrome. In this gene, three variants were present, all with a population frequency of 11.5% or higher. Only three of the 19 variants in USH2A, WFS1 (Gene ID ENSG00000109501, OMIM 606201), and MERTK were predicted to be disease causing by all computational tools; however, USH2A and WFS1 are associated with deafness that these patients do not have. The only novel sequence variant predicted to affect protein function and cause the disease was c.2309A>G p. Glu770Gly (NM_006343.2, NP_006334) in MERTK (Table 4). Unknown by the time of detection, the allele frequency for this variant is now reported at 0.000012 in the genome Aggregation Database. Sanger sequencing of exon 17 confirmed the presence of heterozygous MERTK c.2309A>G variant in both siblings and its absence among healthy individuals (n=100) and in the IRD cohort (n=145) from the same geographic area.

Genome-wide genotyping using the SNP array-- The presence of only one novel likely pathogenic MERTK variant in both patients motivated us to search for a second mutation, such as a large genomic deletion. High-resolution genome-wide genotyping was applied to identify CNVs and homozygosity regions shared by the two patients. The SNP array did not show any regions with neutral loss of heterozygosity (regions of homozygosity, ROH) containing retinal genes or common deletions or duplications when standard cut-off values of 5 Mb for ROH and 100 kb for deletions and duplications were used. Manual analysis revealed only one common ROH on chromosome 2q13 (chr2:110,465,111–110,980,346) and one genomic deletion on the same chromosome (chr2:111,891,895–111,980,985; GRCh38). An 89,091 bp deletion present in both siblings encompassed part of the MERTK gene (Appendix 2).

Breakpoint characterization of the MERTK deletion-- CNV analysis of the genotyping data suggested the deletion breakpoints in the upstream region and in intron 7 of the MERTK gene. Considering this information, we designed PCR primers to enable amplification of an allele containing the deletion and wild-type allele (Figure 2A). The mutant allele was amplified with primers ∆-F-∆ mt-R1 (658 bp, Figure 2B) while the wild-type allele was detected with primers ∆-F- ∆ wt-R (386 bp, Figure 2B). Sanger sequencing of the 658 bp PCR fragment revealed breakpoints chr2:111,890,574–111,981,630 yielding in a 91,057 bp deletion in the MERTK gene (c.-8162_1145–1212del, p.?; Figure 2C). The origin of the GCA sequence at the junction of the two breakpoints is unknown because it is present on both ends.

Population frequency of the deletion analyzed with allele-specific PCR-- To estimate frequency of the approximate 91 kb MERTK deletion in the geographically matched control population, we performed a fragment analysis. PCR was conducted with fluorescently labeled forward primer, ∆-F (Table 1) located immediately upstream of the deletion junction, and two reverse primers, ∆ wt-R in the deleted region detecting a wild-type allele (386 bp, Figure 2B,D), and the other reverse primer ∆ mt-R2, resulting in a deletion-specific PCR product (354 bp, Table 1, Figure 2D). The approximate 91 kb MERTK deletion was not detected in 100 control samples or in 145 samples from a cohort of patients with IRDs.

Approximate 91 kb deletion and c.2309A>G represent allelic variants-- Two novel mutations in MERTK were considered to be likely disease causing in two affected individuals although the lack of parental and offspring samples made segregation analysis impossible. To prove that these mutations are allelic, we extracted RNA from peripheral blood of RP116 and performed RT–PCR with forward primers either in exon 7 or 15 and reverse primer in exon 18 (Table 1, Figure 3A). The PCR products of 1292 bp and 372 bp (Figure 3B) were analyzed with Sanger sequencing. In both fragments, only the G nucleotide was seen in position 2309 in exon 17 (Figure 3C). As the PCR fragments could be amplified only from the allele without the deletion, the result confirmed that the mutations are on different alleles, and both patients are compound heterozygous.

Discussion

Mutations in the MERTK gene as a cause of IRDs are rare and explain about 1% cases of autosomal recessive RP [8]. MERTK is known as a member of the MER/AXL/TYRO3-receptor kinase family encoding a transmembrane protein with two fibronectin domains, two immunoglobulin-like domains, and one tyrosine kinase domain. The gene participates in many physiologic processes, such as cell survival, migration, differentiation, and phagocytosis of apoptotic cells. The protein has a crucial role in maintenance of retinal photoreceptors by mediating phagocytosis of shed photoreceptor outer segments [18,19]. Initially, mutations in the MerTK gene were identified in a natural animal model of recessive retinal degeneration, in the Royal College of Surgeons (RCS) rat [20], and later in human patients with IRDs [9,21]. Mutations in MERTK are described in patients with different phenotypes, such as retinitis pigmentosa, rod-cone and cone-rod dystrophy, Leber congenital amaurosis, and retinal dystrophy [5,22-26].

According to the Human Gene Mutation Database (HGMD® Professional) to date, 68 mutations in MERTK have been reported, including 37 missense, 12 splice-site, ten small deletions, insertions/duplications, and small indels, and six gross deletions. In this study, we discovered one novel variant in the tyrosine kinase domain of MERTK. The variant represents a change in a highly conserved nucleotide (phyloP 4.97), and a highly conserved amino acid preserved in 11 species. Computational tools such as MutationTaster, SIFT, and PolyPhen predict c.2309A>G to be disease causing. We also assessed this variant according to 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 [17]. There is a line of factors for highly likely pathogenicity of c.2309A>G such as variant location in functionally active tyrosine kinase domain (PM1); absence from controls and low frequency in the genome Aggregation Database (PM2); coexistence with the second mutation, a large genomic deletion in trans (PM3); and multiple lines of computational evidence support deleterious effect of the variant (PP3).

Gross deletions presented by a 1.73 kb deletion of exon 15 [27], a 9.86 kb deletion of exon 8 [7], a 25 kb deletion encompassing exons 6 to 8 [28], a 91 kb deletion encompassing exons 1 to 7 [11,29], and a large deletion encompassing exon 3 to 19 [30]. Notably, that homozygous approximate 91 kb deletion in the MERTK gene was a founder mutation in the Faroe Islands responsible for about 30% of RP in that population [11]. The breakpoints of the deletion in the present study and the deletion reported by Ostergaard et al. are very similar: They differ by only the retention of the G nucleotide in these patients at the site of the deletion [11]. In a cohort of patients with RP in northern Sweden, no MERTK deletion or the c.2309A>G mutation was detected indicating that these two MERTK variants do not represent founder mutations in this region.

A common feature of three gross deletions present in retinal dystrophy cases was the absence of exon 8 [7,11,28]. An in-frame MERTK deletion of 25 kb (c.845–1450del; p.Ala282_His483del) removes part of the protein, corresponding to two fibronectin type III domains [28]. In the case of the 91 kb deletion in the present study, both fibronectin type III domains and both immunoglobulin-like domains are deleted, but because of the deletion breakpoint is more than 8,000 bp upstream from the start codon, the gene is unlikely to be transcribed aborting protein expression. It was shown that a start codon mutation c.3G>A causing severe RP annuls protein expression that may lead to loss of MERTK function [26]. The two patients in the present study suffered a severe visual handicap during their whole life span (more than 80 years) similar to other patients with gross deletions.

In this case, at the site of the deletion, there were no sequence deviations, such as an insertion in the case of a 25 kb deletion [28] or duplication [7,27]. The mechanism of the MERTK deletion appearance seems to be complex. One hypothesis is a non-homologous recombination between two AluY elements because a complete AluY element insertion was found in the case of the exon 8 deletion [7]. The finding of homology to three intronic regions on chromosome 2, reverse complement compared to wild-type intron 5, MIR4435–2HG, and LINC00152 did not support this hypothesis [28]. As a possible explanation, the authors suggested the insertion of a transposable element located within the gene [28]. We examined the site of the deletion by performing a BLAST search of the flanking junction sequences and detected that the 5′-end sequence is part of the 286 bp repeat AluSx3 (chr2:111,890,318–111,890,603) and the 3′-end is part of the 108 bp repeat L1ME4b (chr2:111,981,585–111,981,692; Appendix 3). The AluSx3 sequence, 8,162 bp upstream from the start codon, was joined directly to a LINE L1 sequence (L1ME4b) in intron 7 6,160 bp downstream from the last nucleotide in exon 7. Notably that the 3′-breakpoint sequence of L1ME4b contained a matrix attachment region (MAR) AT-rich motif near the endpoint which was highlighted in the study of the recurrent deletion of the RB1 gene caused by a non-homologous recombination between a LINE-1HS and a MER21B element [31]. The first example of a non-homologous recombination between Alu and L1 repeats was described in a dystrophin gene causing a 430 kb deletion [32]. Sequencing of the 430 kb deletion ends revealed the presence of Alu sequence in intron 1 at the 5′-end and the L1 sequence in intron 7 at the 3′-end [32].

In summary, we highlight the possibility that MERTK deletions might be more frequent than estimated and screening for large genomic deletions should be done in IRD cases without known genetic cause. We also show that in recessive cases one should consider compound heterozygosity when one mutation is a single nucleotide change and another mutation is a large CNV. Finally, genetic diagnosis is extremely important considering the progress in the gene therapy for defects in MERTK [33,34].

Appendix 1. Supplementary figure 1.

Appendix 2. Supplementary figure 2.

Appendix 3. Supplementary figure 3.

Acknowledgments

The authors thank our study participants for providing clinical histories and DNA samples. This investigation was supported by grants from Stiftelsen Kronprinsessan Margaretas Arbetsnämnd för synskadade, Ögonfonden and Umeå University and Västerbotten County Council in cooperation in the fields of Medicine, Odontology, and Health.

References

  1. Corton M, Nishiguchi KM, Avila-Fernandez A, Nikopoulos K, Riveiro-Alvarez R, Tatu SD, Ayuso C, Rivolta C. Exome sequencing of index patients with retinal dystrophies as a tool for molecular diagnosis. PLoS One. 2013; 8:e65574 [PMID: 23940504]
  2. Nash BM, Wright DC, Grigg JR, Bennetts B, Jamieson RV. Retinal dystrophies, genomic applications in diagnosis and prospects for therapy. Transl Pediatr. 2015; 4:139-63. [PMID: 26835369]
  3. Burstedt MS, Sandgren O, Holmgren G, Forsman-Semb K. Bothnia dystrophy caused by mutations in the cellular retinaldehyde-binding protein gene (RLBP1) on chromosome 15q26. Invest Ophthalmol Vis Sci. 1999; 40:995-1000. [PMID: 10102298]
  4. Kohn L, Burstedt MS, Jonsson F, Kadzhaev K, Haamer E, Sandgren O, Golovleva I. Carrier of R14W in carbonic anhydrase IV presents Bothnia dystrophy phenotype caused by two allelic mutations in RLBP1. Invest Ophthalmol Vis Sci. 2008; 49:3172-7. [PMID: 18344446]
  5. Brea-Fernandez AJ, Pomares E, Brion MJ, Marfany G, Blanco MJ, Sanchez-Salorio M, Gonzalez-Duarte R, Carracedo A. Novel splice donor site mutation in MERTK gene associated with retinitis pigmentosa. Br J Ophthalmol. 2008; 92:1419-23. [PMID: 18815424]
  6. Charbel Issa P, Bolz HJ, Ebermann I, Domeier E, Holz FG, Scholl HP. Characterisation of severe rod-cone dystrophy in a consanguineous family with a splice site mutation in the MERTK gene. Br J Ophthalmol. 2009; 93:920-5. [PMID: 19403518]
  7. Mackay DS, Henderson RH, Sergouniotis PI, Li Z, Moradi P, Holder GE, Waseem N, Bhattacharya SS, Aldahmesh MA, Alkuraya FS, Meyer B, Webster AR, Moore AT. Novel mutations in MERTK associated with childhood onset rod-cone dystrophy. Mol Vis. 2010; 16:369-77. [PMID: 20300561]
  8. Tschernutter M, Jenkins SA, Waseem NH, Saihan Z, Holder GE, Bird AC, Bhattacharya SS, Ali RR, Webster AR. Clinical characterisation of a family with retinal dystrophy caused by mutation in the Mertk gene. Br J Ophthalmol. 2006; 90:718-23. [PMID: 16714263]
  9. Gal A, Li Y, Thompson DA, Weir J, Orth U, Jacobson SG, Apfelstedt-Sylla E, Vollrath D. Mutations in MERTK, the human orthologue of the RCS rat retinal dystrophy gene, cause retinitis pigmentosa. Nat Genet. 2000; 26:270-1. [PMID: 11062461]
  10. Audo IS, Boulanger-Semama E, Mohand-Said S, Démontant V, Condroyer C, Antonio A, Boyard F, El Shamieh S, Sahel J, Zeitz C. MERTK mutatios account for 2% of cases with inherited retinal dystrophy in a large French cohort. Invest Ophthalmol Vis Sci. 2016; 57:660
  11. Ostergaard E, Duno M, Batbayli M, Vilhelmsen K, Rosenberg T. A novel MERTK deletion is a common founder mutation in the Faroe Islands and is responsible for a high proportion of retinitis pigmentosa cases. Mol Vis. 2011; 17:1485-92. [PMID: 21677792]
  12. Balciuniene J, Johansson K, Sandgren O, Wachtmeister L, Holmgren G, Forsman K. A gene for autosomal dominant progressive cone dystrophy (CORD5) maps to chromosome 17p12-p13. Genomics. 1995; 30:281-6. [PMID: 8586428]
  13. Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988; 16:1215 [PMID: 3344216]
  14. Liu X, Jian X, Boerwinkle E. dbNSFP: a lightweight database of human nonsynonymous SNPs and their functional predictions. Hum Mutat. 2011; 32:894-9. [PMID: 21520341]
  15. den Dunnen JT. Sequence Variant Descriptions: HGVS Nomenclature and Mutalyzer. Curr Protoc Hum Genet. 2016; 90:7
  16. den Dunnen JT, Dalgleish R, Maglott DR, Hart RK, Greenblatt MS, McGowan-Jordan J, Roux AF, Smith T, Antonarakis SE, Taschner PE. HGVS Recommendations for the Description of Sequence Variants: 2016 Update. Hum Mutat. 2016; 37:564-9. [PMID: 26931183]
  17. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, Grody WW, Hegde M, Lyon E, Spector E, Voelkerding K, Rehm HL, Committee ALQA. 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:405-24. [PMID: 25741868]
  18. Qingxian L, Qiutang L, Qingjun L. Regulation of phagocytosis by TAM receptors and their ligands. Front Biol (Beijing). 2010; 5:227-37. [PMID: 21057587]
  19. Strick DJ, Vollrath D. Focus on molecules: MERTK. Exp Eye Res. 2010; 91:786-7. [PMID: 20488176]
  20. D’Cruz PM, Yasumura D, Weir J, Matthes MT, Abderrahim H, LaVail MM, Vollrath D. Mutation of the receptor tyrosine kinase gene Mertk in the retinal dystrophic RCS rat. Hum Mol Genet. 2000; 9:645-51. [PMID: 10699188]
  21. Parinot C, Nandrot EF. A Comprehensive Review of Mutations in the MERTK Proto-Oncogene. Adv Exp Med Biol. 2016; 854:259-65. [PMID: 26427420]
  22. Thompson DA, McHenry CL, Li Y, Richards JE, Othman MI, Schwinger E, Vollrath D, Jacobson SG, Gal A. Retinal dystrophy due to paternal isodisomy for chromosome 1 or chromosome 2, with homoallelism for mutations in RPE65 or MERTK, respectively. Am J Hum Genet. 2002; 70:224-9. [PMID: 11727200]
  23. Shahzadi A, Riazuddin SA, Ali S, Li D, Khan SN, Husnain T, Akram J, Sieving PA, Hejtmancik JF, Riazuddin S. Nonsense mutation in MERTK causes autosomal recessive retinitis pigmentosa in a consanguineous Pakistani family. Br J Ophthalmol. 2010; 94:1094-9. [PMID: 20538656]
  24. Coppieters F, Van Schil K, Bauwens M, Verdin H, De Jaegher A, Syx D, Sante T, Lefever S, Abdelmoula NB, Depasse F, Casteels I, de Ravel T, Meire F, Leroy BP, De Baere E. Identity-by-descent-guided mutation analysis and exome sequencing in consanguineous families reveals unusual clinical and molecular findings in retinal dystrophy. Genet Med. 2014; 16:671-80. [PMID: 24625443]
  25. Srilekha S, Arokiasamy T, Srikrupa NN, Umashankar V, Meenakshi S, Sen P, Kapur S, Soumittra N. Homozygosity Mapping in Leber Congenital Amaurosis and Autosomal Recessive Retinitis Pigmentosa in South Indian Families. PLoS One. 2015; 10:e0131679 [PMID: 26147992]
  26. Jinda W, Poungvarin N, Taylor TD, Suzuki Y, Thongnoppakhun W, Limwongse C, Lertrit P, Suriyaphol P, Atchaneeyasakul LO. A novel start codon mutation of the MERTK gene in a patient with retinitis pigmentosa. Mol Vis. 2016; 22:342-51. [PMID: 27122965]
  27. Siemiatkowska AM, Arimadyo K, Moruz LM, Astuti GD, de Castro-Miro M, Zonneveld MN, Strom TM, de Wijs IJ, Hoefsloot LH, Faradz SM, Cremers FP, den Hollander AI, Collin RW. Molecular genetic analysis of retinitis pigmentosa in Indonesia using genome-wide homozygosity mapping. Mol Vis. 2011; 17:3013-24. [PMID: 22128245]
  28. Evans DR, Green JS, Johnson GJ, Schwartzentruber J, Majewski J, Beaulieu CL, Qin W, Marshall CR, Paton TA, Roslin NM, Paterson AD, Fahiminiya S, French J, Boycott KM, Woods MO, Consortium FC. Novel 25 kb Deletion of MERTK Causes Retinitis Pigmentosa With Severe Progression. Invest Ophthalmol Vis Sci. 2017; 58:1736-42. [PMID: 28324114]
  29. Ellingford JM, Barton S, Bhaskar S, O’Sullivan J, Williams SG, Lamb JA, Panda B, Sergouniotis PI, Gillespie RL, Daiger SP, Hall G, Gale T, Lloyd IC, Bishop PN, Ramsden SC, Black GCM. Molecular findings from 537 individuals with inherited retinal disease. J Med Genet. 2016; 53:761-7. [PMID: 27208204]
  30. Boulanger-Scemama E, El Shamieh S, Demontant V, Condroyer C, Antonio A, Michiels C, Boyard F, Saraiva JP, Letexier M, Souied E, Mohand-Said S, Sahel JA, Zeitz C, Audo I. Next-generation sequencing applied to a large French cone and cone-rod dystrophy cohort: mutation spectrum and new genotype-phenotype correlation. Orphanet J Rare Dis. 2015; 10:85 [PMID: 26103963]
  31. Albrecht P, Bode J, Buiting K, Prashanth AK, Lohmann DR. Recurrent deletion of a region containing exon 24 of the RB1 gene caused by non-homologous recombination between a LINE-1HS and MER21B element. J Med Genet. 2004; 41:e122 [PMID: 15591264]
  32. Suminaga R, Takeshima Y, Yasuda K, Shiga N, Nakamura H, Matsuo M. Non-homologous recombination between Alu and LINE-1 repeats caused a 430-kb deletion in the dystrophin gene: a novel source of genomic instability. J Hum Genet. 2000; 45:331-6. [PMID: 11185740]
  33. Ghazi NG, Abboud EB, Nowilaty SR, Alkuraya H, Alhommadi A, Cai H, Hou R, Deng WT, Boye SL, Almaghamsi A, Al Saikhan F, Al-Dhibi H, Birch D, Chung C, Colak D, LaVail MM, Vollrath D, Erger K, Wang W, Conlon T, Zhang K, Hauswirth W, Alkuraya FS. Treatment of retinitis pigmentosa due to MERTK mutations by ocular subretinal injection of adeno-associated virus gene vector: results of a phase I trial. Hum Genet. 2016; 135:327-43. [PMID: 26825853]
  34. Dalkara D, Goureau O, Marazova K, Sahel JA. Let There Be Light: Gene and Cell Therapy for Blindness. Hum Gene Ther. 2016; 27:134-47. [PMID: 26751519]
  35. Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, O’Donnell-Luria AH, Ware JS, Hill AJ, Cummings BB, Tukiainen T, Birnbaum DP, Kosmicki JA, Duncan LE, Estrada K, Zhao F, Zou J, Pierce-Hoffman E, Berghout J, Cooper DN, Deflaux N, DePristo M, Do R, Flannick J, Fromer M, Gauthier L, Goldstein J, Gupta N, Howrigan D, Kiezun A, Kurki MI, Moonshine AL, Natarajan P, Orozco L, Peloso GM, Poplin R, Rivas MA, Ruano-Rubio V, Rose SA, Ruderfer DM, Shakir K, Stenson PD, Stevens C, Thomas BP, Tiao G, Tusie-Luna MT, Weisburd B, Won HH, Yu D, Altshuler DM, Ardissino D, Boehnke M, Danesh J, Donnelly S, Elosua R, Florez JC, Gabriel SB, Getz G, Glatt SJ, Hultman CM, Kathiresan S, Laakso M, McCarroll S, McCarthy MI, McGovern D, McPherson R, Neale BM, Palotie A, Purcell SM, Saleheen D, Scharf JM, Sklar P, Sullivan PF, Tuomilehto J, Tsuang MT, Watkins HC, Wilson JG, Daly MJ, MacArthur DG, Exome Aggregation C. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016; 536:285-91. [PMID: 27535533]
  36. Tavtigian SV, Deffenbaugh AM, Yin L, Judkins T, Scholl T, Samollow PB, de Silva D, Zharkikh A, Thomas A. Comprehensive statistical study of 452 BRCA1 missense substitutions with classification of eight recurrent substitutions as neutral. J Med Genet. 2006; 43:295-305. [PMID: 16014699]
  37. Ng PC, Henikoff S. SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003; 31:3812-4. [PMID: 12824425]
  38. Schwarz JM, Rodelsperger C, Schuelke M, Seelow D. MutationTaster evaluates disease-causing potential of sequence alterations. Nat Methods. 2010; 7:575-6. [PMID: 20676075]
  39. Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR. A method and server for predicting damaging missense mutations. Nat Methods. 2010; 7:248-9. [PMID: 20354512]