- Primary research
- Open Access
A novel LRAT mutation affecting splicing in a family with early onset retinitis pigmentosa
- Yabin Chen†1,
- Li Huang†2,
- Xiaodong Jiao1,
- Sheikh Riazuddin3, 4, 5,
- S. Amer Riazuddin†6 and
- J. Fielding Hetmancik†1Email author
© The Author(s). 2018
Received: 24 January 2018
Accepted: 18 June 2018
Published: 4 July 2018
Background and purpose
Retinitis pigmentosa is an important cause of severe visual dysfunction. This study reports a novel splicing mutation in the lecithin retinol acyltransferase (LRAT) gene associated with early onset retinitis pigmentosa and characterizes the effects of this mutation on mRNA splicing and structure.
Genome-wide linkage analysis followed by dideoxy sequencing of the linked candidate gene LRAT was performed in a consanguineous Pakistani family with autosomal recessive retinitis pigmentosa. In silico prediction and minigene assays were used to investigate the effects of the presumptive splicing mutation.
ARRP in this family was linked to chromosome 4q31.21-q32.1 with a maximum LOD score of 5.40. A novel homozygous intronic mutation (NM_004744.4: c.541-15T>G) was detected in LRAT. In silico tools predicted that the AG-creating mutation would activate an intronic cryptic acceptor site, but cloning fragments of wild-type and mutant sequences of LRAT into Exontrap Cloning Vector pET01 and Expression Cloning Vector pCMV-(DYKD4K)-C showed that the primary effect of the sequence change was to weaken the nearby authentic acceptor site and cause exon skipping, with only a small fraction of transcripts utilizing the acceptor site producing the reference transcript.
The c.541-15T>G mutation in LRAT results in aberrant splicing and is therefore predicted to be causal for the early onset retinitis pigmentosa in this family. In addition, this work suggests that minigenes adapted to the specific gene and exon may need to be designed for variants in the first and last exon and intron to mimic the authentic splicing mechanism in vivo.
Retinitis pigmentosa (RP, [MIM 268000]) is a clinically and genetically heterogeneous disorder affecting approximately 1 in 4000 individuals worldwide . Clinically, patients initially exhibit night blindness followed by progressive loss of peripheral visual fields and eventually complete loss of central vision. Typical fundus changes include bone spicule-like pigmentation in the mid-peripheral retina, waxy pallor of the optic discs, and attenuation of retinal blood vessels. Since RP initially affects the rod photoreceptors, followed by the degeneration of cone photoreceptors, patients often have severely diminished or extinguished rod response in electroretinography (ERG) recordings even in early stages, while the cone response is relatively preserved initially but becomes undetectable as the disease progresses . The genetic inheritance patterns of RP include autosomal-dominant (about 30–40% of cases), autosomal-recessive (50–60%), and X-linked (5–15%) [3, 4] inheritance. More than 82 causative genes have been identified for RP so far, of which 58 genes have been identified in families with autosomal recessive RP (arRP; RetNet).
Lecithin retinol acyltransferase (LRAT) is a retinyl ester synthase, catalyzing the formation of fatty acid retinyl esters, a crucial step in the retinoid cycle. In the eye, it is specifically expressed in the retinal pigmented epithelium (RPE), so that LRAT dysfunction leads to diminished visual chromophores and eventual retinal degeneration . Mutations in the LRAT gene can cause Leber congenital amaurosis, juvenile retinitis pigmentosa, and early-onset severe retinitis pigmentosa with autosomal recessive inheritance (MIM 604863). Leber congenital amaurosis is the most severe retinal dystrophic disease. Patients usually present in the first decade of life with severe visual impairment and pendular nystagmus .
It is estimated that about 10% of disease-causing mutations affect splicing . Most of the splice site mutations affect the invariant GT or AG dinucleotides in the 3′ and 5′ splice sites, respectively [8, 9]. However, mutations involving other positions of the 5′ or 3′ splice sites can also impair splicing and typically lead to exon skipping, activation of a cryptic splice site, or intron retention . This study demonstrates that single base pair substitution in the intron even 15 bases upstream of 3′ splice site can disrupt spliceosomal recognition of the splice site.
Materials and methods
This study was approved by the Institutional Review Boards (IRB) of the National Centre of Excellence in Molecular Biology, Lahore, Pakistan, and the CNS IRB at the National Institutes of Health. Participating individuals or their guardians gave written informed consent consistent with the tenets of the Declaration of Helsinki before the study. Family 61254 is a consanguineous Pakistani family with non-syndromic RP. All participants underwent a thorough family, ophthalmic, and medical history, and selected individuals were assessed by best-corrected visual acuity, slit-lamp biomicroscopy, fundus photography, and electroretinography (ERG). Blood samples were collected from potentially informative family members, and genomic DNA was extracted from leukocytes according to standard protocols .
Genome-wide linkage analysis
We designed linkage mapping panels to complete a genome-wide scan for family 61254. Primer sequences and PCR conditions are shown in Additional file 1: Table S1. Based on the initial results of the linkage analysis, markers with logarithm of the odds (LOD) scores greater than 2 were selected for further fine mapping by microsatellite markers closely spaced 1–2 cM apart.
Two-point linkage analyses were performed with alleles of family 61254 obtained through the genome-wide scan and fine mapping using the FASTLINK modification of the MLINK program in the LINKAGE program package [12, 13] using the pedegree file shown in Additional file 2. Maximum LOD scores were calculated using ILINK. Autosomal recessive RP was analyzed as a fully penetrant trait with an affected allele frequency of 0.0001. Haplotypes were drawn by the Cyrillic 2.1 program (Cyrillic Software, Wallingford, Oxfordshire, UK) and confirmed by inspection.
Primer pairs for the coding exons and 100 bp of flanking intronic regions of LRAT were designed using the Primer 3 program (http://primer3.ut.ee/) and are shown in Additional file 3: Table S3. Amplifications were completed in 10 μl reactions containing 40 ng of genomic DNA, 0.5 μl of 5 μM of each primer, 1 μl GeneAmp 10×PCR Gold buffer (Applied Biosystems), 0.6 μl of 25 mM MgCl2, 0.8 μl of 2.5 mM dNTP mixture, and 1 μl of 5 U/μl Taq DNA polymerase. PCR reactions consisted of a denaturation step at 95 °C for 5 min followed by a two-step touchdown procedure. The first 15 cycles consisted of denaturation at 95 °C for 30 s, followed by a touchdown annealing step for 30 s in which the annealing temperature was initially set at 64 °C and decreased by 0.5 °C per cycle, and an elongation step at 72 °C for 30 s. The second step of 20 cycles consisted of denaturation at 95 °C for 30 s followed by annealing at 57 °C for 30 s and elongation at 72 °C for 30 s. These were followed by a final elongation at 72 °C for 10 min. PCR products were purified using the AMPure XP system (Beckman coulter Biomek NX, Brea, CA). The PCR primers for each exon were used for bidirectional sequencing using the BigDye Terminator Ready Reaction mix (Applied Biosystems), according to the manufacturer’s instructions. The sequencing products were purified using the Agencourt CleanSEQ system (Beckman Coulter Biomeck NX). Sequencing was performed on an ABI PRISM 3130 Automated sequencer (Applied Biosystems, Foster City, CA). Sequencing results were analyzed using Mutation Surveyor (SoftGenetics, State College, PA) and Seqman software Lasergene 12.2.0 (DNASTAR, Madison, WI).
In silico splice site prediction
To evaluate the potential pathogenicity of the splice site variation, several in silico splice site prediction programs were used, including Human Splice Finder (http://www.umd.be/HSF3/), SpliceView (http://bioinfo.itb.cnr.it/oriel/splice-view.html), Berkeley Drosophila Genome Project Neural Network (BDGP: http://www.fruitfly.org/), and NetGene2 (http://www.cbs.dtu.dk/services/NetGene2/). Default settings are applied for all analyses.
Generation of vector constructs
Minigene assays were performed to confirm the predicted pathogenic effects of the LRAT mutation. Since mutation c.541-15T>G (NM_004744.4) was in intron 2, a 5228-bp genomic DNA fragment that included exon 2, intron 2, and part of exon 3 of LRAT was amplified using primers LRAT-gDNA-F and LRAT-gDNA-R. Because limited minigene amplicon length can be included in the minigene assay, and exon 3 is 4136-bp long with 153-bp coding sequence and 3983-bp 3′ UTR sequence, only part of exon 3 (402-bp) containing the coding region and several predicted downstream cryptic splice donor and acceptor sites were included in the 5228-bp fragment used for the minigene assay. The wildtype and mutant DNA fragments were amplified from genomic DNA of unaffected sibling 6125441 and patient 6125408 carrying the wild-type and homozygous mutant alleles, respectively. PCR products were cloned into a TA vector (pCR®II TOPO) and then subcloned into Exontrap cloning Vector pET01 (MoBiTec GmbH, Göttingen, Germany) via BamHI and NotI restriction sites to generate the wildtype WT-LRAT-pET01 and mutant Mut-LRAT-pET01 constructs.
Because exon 3 is the last exon of LRAT, to mimic authentic splicing mechanism in vivo, we also used a Q5 Site-Directed Mutagenesis Kit (New England Biolabs, Ipswich, MA) to abolish the splice acceptor site of 3′ exon on pET01 by altering the final acceptor AG to AT, creating the modified empty construct pET01-2, modified wildtype construct WT-2-LRAT-pET01, and modified mutant construct Mut-2-LRAT-pET01. In addition, the same 5228-bp genomic DNA fragment was also amplified by PCR with primers containing the appropriate restriction enzyme sites (LRAT-gDNA-F-SalI and LRAT-gDNA-R-KpnI) and then cloned into the expression vector pCMV-(DYKD4K)-C (Clontech Laboratories, Inc. CA, USA), which has no vector exons before and behind the inserted PCR products. SalI and KpnI restriction sites were used to insert the PCR products and construct wildtype construct WT-LRAT-pCMV and mutant construct Mut-LRAT-pCMV.
Two-point LOD scores of markers around 4q31.21-q32.1
Cell culture and transfection
Human embryonic kidney 293T cells and human retinal pigment epithelial cells (ARPE-19) were maintained in Dulbecco’s modified Eagle medium supplemented with 10% fetal bovine serum in 5% CO2 at 37 °C. The Exontrap cloning vectors pET01, WT-LRAT-pET01, Mut-LRAT-pET01, pET01-mut, WT-2-LRAT-pET01, and Mut-2-LRAT-pET01 were transfected into 293T cells using PolyJet In Vitro DNA Transfection Reagent (SignaGen Laboratories, MD). The expression vectors pCMV-(DYKD4K)-C, WT-LRAT-pCMV, and Mut-LRAT-pCMV were transfected into 293T cells and ARPE-19 cells, respectively, using PolyJet In Vitro DNA Transfection Reagent. For nonsense-mediated decay inhibition studies, non-transfected and transfected ARPE19 cells were treated with 0, 5 μM, and 10 μM PTC124 separately (Enzo Life Sciences, Inc., NY) for 24 h. Cells were harvested 24 h after transfection.
RNA isolation, cDNA synthesis, PCR, and sequencing
Total cellular RNA was extracted using Trizol reagent (Life Technologies), and cDNA was synthesized using a reverse transcriptase kit (Invitrogen, Carlsbad, CA) with oligo (dT)20 primers. Primers used for PCR amplification for cDNA are listed in Additional file 2: Table S2. The amplified products were separated by electrophoresis on 10% Novex TBE gels (Invitrogen) and stained with ethidium bromide. Visualized bands were extracted from the gels, and the same primers were used for sequencing to validate the splice effects in the in vitro assays.
During a genome-wide scan for family 61254, a LOD score greater than 2.0 was found only for markers D4S1575, D4S413, D14S261, and D17S921. Except for a single chromosome 4 locus, which was narrowed from D4S429 to D4S424, other regions were excluded by fine mapping markers at 1–2 cM intervals and haplotype inspection. Among the markers selected for fine mapping of the chromosome 4 locus, D4S2998 obtained a maximum LOD score of 5.40 at θ = 0 (Table 1). Thus, two-point linkage mapping in family 61254 identified a linked region of 14.09 cM (18.89 Mb) on chromosome 4q31.21-q32.1 flanked by markers D4S424 and D4S2982 (Table 1, Fig. 1a).
Because there is no exon following exon 3 in human genomic DNA of LRAT, we abolished the splice acceptor site of the pre-proinsulin 3′ exon on pET01 to reduce its effect on mRNA splicing of the inserted sequence of LRAT. Cells transfected with WT-2-LRAT-pET01 still produced mRNA with correct splicing at the 5′ end of exon 3, and cells transfected with Mut-2-LRAT-pET01 still produced mRNA skipping the exon 3. However, because the authentic acceptor site of the pre-proinsulin 3′ exon was abolished, the transcription utilized a cryptic acceptor site within the 3′ exon (Fig. 5).
We have identified a c.541-15T>G mutation in LRAT that results in aberrant splicing as the cause of autosomal recessive early-onset retinitis pigmentosa in a Pakistani family. The mammalian 3′ splice site consensus consists of three critical parts: YAG (Y = pyrimidine) at positions − 1 to − 3 relative to the first nucleotide of exon, a polypyrimidine tract (PPT) starting at position − 5 and extending 10 or more nucleotides into the intron, and the branch point sequence (BPS; mammalian consensus YNCURAY, where Y = pyrimidine, R = purine, and N = any nucleotide) located upstream of the PPT, usually 11–40 nucleotides from the YAG [14, 15]. The PPT is a crucial recognition element for branch site definition prior to splicing, and the factor responsible for this early PPT recognition is U2AF,65 a subunit of the U2 snRNP auxiliary factor (U2AF) . Since the c.541-15T>G mutation replaces a pyrimidine with a purine at the 5′ end of the PPT, it might possibly affect the interaction of the U2AF65 with the PPT, thus weakening activity of the authentic splice acceptor site. Both bioinformatic tools and minigene assays showed that the novel c.541-15T>G mutation in the LRAT gene causes aberrant splicing and is therefore causal for early onset retinitis pigmentosa, further broadening the mutation spectrum in the LRAT gene.
Two in silico prediction tools, Human Splice Finder and SpliceView, predicted creation of a new splice acceptor site by the mutation, consistent with the AG/GT rule. Minigene assays using both vector pET01 and pCMV-(DYKD4K)-C based constructs showed that the resulting mutant cDNA skipped exon 3 of LRAT. The wild-type transcript was absent in 293T cells transfected with the pET01 construct containing the mutant fragment of LRAT, and identified at reduced levels in 293T cells transfected with mutated pCMV-(DYKD4K)-C construct compared with cells transfected with the wildtype pCMV construct. This suggests that the mutation does not completely inactivate the authentic splice site but rather suppresses it so that some WT LRAT RNA might be present, consistent with the concept of leaky splicing . This RNA level might not be sufficient to maintain adequate LRAT function in the retina so that patients homozygous for this mutation might manifest a RP phenotype. Although the minigene assays are not in perfect accordance with the in silico predictions about splicing patterns, all indicated that the mutation would result in aberrant splicing.
Minigene assays provide a rapid way to evaluate the impact of splice site mutations, especially when isolation and characterization of mRNA from the original source is technically difficult or ethically impossible. LRAT is expressed in RPE cells and is important for the retinoid cycle . Since it is not expressed in peripheral blood cells and retinal biopsy is too invasive for patients, minigene assays provided an optimal approach to studying the in vitro effect of LRAT spice site mutations. However, the system also does not always reflect the exact splicing patterns of in vivo studies. While it allows prediction of whether a mutation is going to influence splicing, several studies have found occasional differences in splice patterns between minigene and RNA analyses in patient samples [17–19]. This might well relate to differences between the vector sequence surrounding the cloned insert and the genomic DNA in vivo. In an attempt to address this, we abolished the splice acceptor site of 3′ exon on pET01 and also used the pCMV-(DYKD4K)-C vector to investigate further the splicing effect of the mutation, aiming to weaken the influence of vector sequence on the cloned insert. The results confirmed that the mutation could cause splicing impairment compared with the wildtype sequence.
In conclusion, bioinformatics tools may be useful for predicting the pathogenicity of variants, but minigene assays are more reliable in the evaluation of the potential effects of splice-site mutations. Adapted minigenes may need to be designed for variants in the first and last exon and introns to mimic the authentic splicing mechanism in vivo. Despite the limitations of both in silico predictions and minigene assays, it is likely that this mutation c.541-15T>G is pathogenic, when all the information including linkage and co-segregation in the pedigree excluding other candidates, the known role of LRAT in retinal degeneration, in silico analyses predicting splicing impairment, and two minigene assays are consistent with weakening of the authentic splice acceptor site resulting in a truncated LRAT protein and arRP.
The authors are thankful to all the family members for their participation in this study.
This work was supported by NEI project EYE000272.
Availability of data and materials
All relevant data are included in the manuscript and Additional files.
YC carried out the experiments, assisted with the analysis, and helped write the manuscript. LH carried out the experiments, assisted with the analysis, and helped write the manuscript. XDJ carried out the experiments, assisted with the analysis, and helped write the manuscript. SR helped plan the project. SAR helped plan the project, conceive of the experiments, analyze the data, and write the manuscript. JFH helped plan the project, conceive of the experiments, analyze the data, and write the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
This study was approved by the Institutional Review Boards (IRB) of the National Centre of Excellence in Molecular Biology, Lahore, Pakistan, and the CNS IRB at the National Institutes of Health. Participating individuals or their guardians gave written informed consent to participate and to publish consistent with the tenets of the Declaration of Helsinki before the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Hartong DT, Berson EL, Dryja TP. Retinitis pigmentosa. Lancet. 2006;368:1795–809.View ArticlePubMedGoogle Scholar
- Bird AC. Retinal photoreceptor dystrophies LI. Edward Jackson memorial lecture. AmJOphthalmol. 1995;119:543–62.Google Scholar
- Bunker CH, Berson EL, Bromley WC, Hayes RP, Roderick TH. Prevalence of retinitis pigmentosa in Maine. AmJOphthalmol. 1984;97:357–65.Google Scholar
- Rivolta C, Sharon D, DeAngelis MM, Dryja TP. Retinitis pigmentosa and allied diseases: numerous diseases, genes, and inheritance patterns. HumMolGenet. 2002;11:1219–27.Google Scholar
- Sears AE, Palczewski K. Lecithin:retinol acyltransferase: a key enzyme involved in the retinoid (visual) cycle. Biochemistry. 2016;55:3082–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Koenekoop RK. An overview of Leber congenital amaurosis: a model to understand human retinal development. SurvOphthalmol. 2004;49:379–98.Google Scholar
- Lopez-Bigas N, Audit B, Ouzounis C, Parra G, Guigo R. Are splicing mutations the most frequent cause of hereditary disease? FEBS Lett. 2005;579:1900–3.View ArticlePubMedGoogle Scholar
- Krawczak M, Thomas NS, Hundrieser B, Mort M, Wittig M, Hampe J, Cooper DN. Single base-pair substitutions in exon-intron junctions of human genes: nature, distribution, and consequences for mRNA splicing. Hum Mutat. 2007;28:150–8.View ArticlePubMedGoogle Scholar
- Krawczak M, Reiss J, Cooper DN. The mutational spectrum of single base-pair substitutions in mRNA splice junctions of human genes: causes and consequences. Hum Genet. 1992;90:41–54.View ArticlePubMedGoogle Scholar
- Ward AJ, Cooper TA. The pathobiology of splicing. J Pathol. 2010;220:152–63.PubMedPubMed CentralGoogle Scholar
- Smith RJH, Holcomb JD, Daiger SP, Caskey CT, Pelias MZ, Alford BR, Fontenot DD, Hejtmancik JF. Exclusion of Usher syndrome gene from much of chromosome 4. CytogenetCell Genet. 1989;50:102–6.View ArticlePubMedGoogle Scholar
- Lathrop GM, Lalouel JM. Easy calculations of lod scores and genetic risks on small computers. AmJ Hum Genet. 1984;36:460–5.Google Scholar
- Schaffer AA, Gupta SK, Shriram K, Cottingham RW. Avoiding recomputation in genetic linkage analysis. Hum Hered. 1994;44:225–37.View ArticlePubMedGoogle Scholar
- Moore MJ. Intron recognition comes of AGe. Nat Struct Biol. 2000;7:14–6.View ArticlePubMedGoogle Scholar
- Peled-Zehavi H, Berglund JA, Rosbash M, Frankel AD. Recognition of RNA branch point sequences by the KH domain of splicing factor 1 (mammalian branch point binding protein) in a splicing factor complex. Mol Cell Biol. 2001;21:5232–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Caminsky N, Mucaki EJ, Rogan PK. Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis. F1000Res. 2014;3:282.PubMedPubMed CentralGoogle Scholar
- Schneider B, Koppius A, Sedlmeier R. Use of an exon-trapping vector for the evaluation of splice-site mutations. Mamm Genome. 2007;18:670–6.View ArticlePubMedGoogle Scholar
- Bonnet C, Krieger S, Vezain M, Rousselin A, Tournier I, Martins A, Berthet P, Chevrier A, Dugast C, Layet V, et al. Screening BRCA1 and BRCA2 unclassified variants for splicing mutations using reverse transcription PCR on patient RNA and an ex vivo assay based on a splicing reporter minigene. J Med Genet. 2008;45:438–46.View ArticlePubMedGoogle Scholar
- Acedo A, Sanz DJ, Duran M, Infante M, Perez-Cabornero L, Miner C, Velasco EA. Comprehensive splicing functional analysis of DNA variants of the BRCA2 gene by hybrid minigenes. Breast Cancer Res. 2012;14:R87.View ArticlePubMedPubMed CentralGoogle Scholar