Three patients of a consanguineous family from the south of Morocco were referred to our Department of Medical Genetics at the National of Health in Rabat for Corneal dystrophy. A five generational pedigree diagram was constructed after a thorough interview of the affected mother (III-2), the available three family members affected (IV-7, IV-8 and III-2) including two phenotypically unaffected (IV-1, IV-3). An autosomal dominant mode of inheritance (Fig. 1). Clinical features of patients are given below.
The patient 1 (IV-7) was a 34-year-old woman who complained of recurrent episodes of corneal pain since the age of 18, associated with a decrease in visual acuity. Her vision was 4/10 OD and 6/10 OS. Slit-lamp examinations revealed, in the two eyes, anterior epithelial and stromal corneal dystrophy, in the form of spaced microvacuoles by a heterogeneous thickening of the epithelium due to thickening of an abnormal subepithelial fibrous layer and poorly individualized anterior opacities, with fuzzy edges (Fig. 2). The rest of the examination of the anterior segment is normal, including a normal iris, a clear lens, and an intraocular pressure at 17mmHg in ODG. The Optical Coherence Tomography (OCT) scan shows an unevenness of the epithelial layers by a homogeneous confluent layer of hyper-reflected deposits with a serrated anterior border taking the sawtooth appearance, replacing the Bowman layer and reaching the anterior stroma (Fig. 3 (A , B, C, D)). It is thicker and becomes thinner on the periphery and disappears towards the limb. The pachymetry is 512 μm in OD and 523 μm in OS.
The patient 2 (IV-8) is a woman 44 years old; she presents episodes of recurrent keratitis with a decrease in visual acuity of progressive installation since the age of 20 years. Visual acuity in OD shows that she counts the fingers at 3m / OS at 2/10 (irremovable). Examination with the slit lamp shows in OD and OS anterior epithelial and stromal corneal dystrophy in the form of microvacuoles, especially in the periphery, with a heterogeneous thickening of the epithelium, more pronounced in the center giving a large central opacity with fuzzy edges at both eyes (Fig. 4 (A, B) ). Clear lens; ocular tone at 15 mmHg in OD 14 mmHg in OS, and normal fundus. The OCT of the cornea shows an unevenness of the anterior epithelial and stromal layers with thicker hyper-reflective deposits and a clear central opacity in both eyes (Fig. 5 (A, B, C, D)), As well as a pachymetry at 538 μm in OD, 543 μm in OS.
The patient 3 (III-2) is a woman 72 years old, she has bilateral osteoarthritis, and she reports episodes of recurrent corneal pain. As for her visual acuity, in far vision OD, she can barely see the movement of the fingers / OS: (she counts the fingers at 3m). The Slit lamp examination showed in the right eye, central and paracentric yellowish and gelatinous central and paracentric deposition associated with a corneal opacity deeper than the previous one and affecting the epithelial layers of the Bowman's membrane and the anterior stroma in the form of an epithelial fibrous layer, taking almost the entire corneal surface (Fig. 6). The rest of the examination is hampered by the very important dystrophy of the right eye, (i.e., ineclairable fundus). In the left eye, the slit lamp examination showed CD involving the epithelial layers of the Bowman's membrane and the anterior stroma in the form of a heterogeneous thickening more important at the center, associated with opacity with fuzzy boundaries and a corticonuclear cataract (Fig. 6). In the Fundus, the pupillary glow with flattened retina was noted. The OCT of the cornea shows in the right eye: the form of microvacuoles as thin patches (blue arrows) associated with a significant loss of epithelial cells with disorganization of the epithelial layers, the absence of Bowman's membrane, and an anterior intrastromal bubble (white arrows). The Pachymetry is at 586 μm (Fig. 7). In the Left eye, the OCT shows an irregularity of the corneal surface with disorganization of the epithelial layers, a discontinuous Bowman’s membrane, and anterior stromal reshaping especially in the center (Fig. 7), the process of corneal transplantation, in this patient, is ongoing.
After obtaining written consent from all participants involved in the study, we collected blood samples from family members. Blood samples were collected from five members of the family. The extraction of genomic DNA from whole peripheral blood by using QIAamp DNA Blood Mini Kit (Qiagen Valencia, CA), strictly following the manufacturer’s protocol. Whole-exome sequencing (WES) was performed for the three patients; 500ng of fragmented DNA by enzymatic fragmentation with Kapa Hyper Plus Kit (KapaBiosystems Inc. Wilmington, MA, USA) was amplified according to the manufacturer instructions and was subjected to enrichment with SeqCap EZ Human Exome v3.0 Roche Nimblegen (Roche, Basel, Switzerland).
The Illumina HiSeq 2500 system was used to sequence 64 enriched megabases in a fast-running double-ended mode (2x100bp). Use bcl2fastq v1.8.4 (Illumina) to convert the original data (bcl file) into a fastq file. The sequence was analyzed according to GATK best practice recommendations; BWA-MEM was used for mapping, and GATK (haplotype calling program) was used for variant calling. Adopt Variant Studio (Illumina) for annotation and filtering steps.
Candidate variants were selected using an autosomal dominant mode of inheritance according to the following criteria: i) heterozygous variants, ii) non-synonymous variants, iii) Variants predicted pathogenic or likely pathogenic, iv) Variants with Minor allele frequency (MAF) of <0.01 were selected from the 1000 Genome Project (http://www.1000genomes.org/) and ESP6500 exome project (https://evs.gs.washington.edu/EVS/) and v) segregation analysis. Candidate variants in favor of clinical manifestations and passing the criteria were validated by Sanger sequencing.
The established variants were cross-checked with the 1000 genomes database (http://www.1000genomes.org/), with the Exome Variant Server (http://evs.gs.washington.edu/EVS/), HGMD (http://www.biobase-international. com/product/hgmd) and with the «clinvar» database (http://www.ncbi.nlm.nih.gov/clinvar/).
To confirm the mutation detected by exome sequencing, standard PCR to amplify the exon 10 of TGFBI gene was carried by using the TGFBI_F:5’-GACCAGGCTAATTACCATTCTTG-3’ and TGFBI_R:5’-TGAGATATGTCCTGGAGCCC-3’ primer pair. Amplification products were electrophoresed on 1% agarose gel. Sanger sequencing was performed with dye terminator chemistry (ABI Prism BigDye v3.1) and run on automated sequencer Using 3130 Genetic Analyzer (Thermo Fisher Scientific). The results obtained were aligned with the reference genome (GRCh37/hg19) and then analyzed by DNA variant analysis software (Mutation Surveyor® software).
The transmission form of this family was consistent with an autosomal dominant inheritance (Fig. 1). Clinical examinations demonstrated bilaterally multiple superficial, epithelial and stromal anterior granular opacities; in different stages of severity among three patients of this family. The three patients with the same mutation who shared a mixed phenotype with a superficial form of Granular Corneal Dystrophy (GCD) type 1 and TBCD patterns.
Whole Exome Sequencing and Variant Validation: WES was performed on three family members (IV-7, IV-8 and III-2). The analysis of the data of the three patients shows 30887 variants in 19712 genes. After filtering, this number was reduced to a single allelic variant in the TGFBI gene Filtering of variants is illustrated in Fig. 8. Only a heterozygous mutation (c.1772C>A; p.Ser591Tyr) in the TGFBI gene, identified in these three patients, was proposed as the potential pathogenic mutation within this family. (Fig. 9 A)
Bioinformatics analysis using Polyphen-2 and SIFT suggested that this mutation was predicted to be probably damaging. A high degree of conservation of this amino acid was demonstrated by a score of 5.52 calculated by Genomic Evolutionary Rate Profiling (GERP) http://mendel.stanford.edu/SidowLab/downloads/gerp/index.html. Comparative amino acid sequence alignment of TGFBI across different species revealed that this mutation occurred at highly conserved positions (Fig. 9 B) (https://www.ensembl.org/Homo_sapiens/Tools/Blast). The mutation was confirmed by Sanger sequencing in affected sister and mother and two unaffected brothers.
Furthermore, the variant was not present in the GnomAD browser (accessed July 2020), nor in our in-house WES database of 100 unrelated individuals of Moroccan ethnicity that had been carried out for diseases other than corneal dystrophy.
The genomic and clinical data both supported a diagnosis of Thiel-Behnke corneal dystrophy in this family.