DOI: https://doi.org/10.21203/rs.3.rs-16667/v2
Background: Usher syndrome is a king of phenotypic and genetic heterogeneous disease. Our purpose was to identify the gene mutation in a Chinese family with Usher syndrome type 2 and describe the clinical features.
Case presentation: A 23-year-old man complained of 10-year nyctalopia and a 3-year decline in visual acuity of both eyes accompanied by congenital dysaudia. To clarify the diagnosis, the clinical symptoms were observed and analysed in combination with comprehensive ophthalmologic examinations as well as genetic analysis (Targeted exome sequencing, TES). Typical clinical presentation of Usher syndrome on fundus was found including a wax yellow-like disc, bone-spicule formations and retinal vessel stenosis. Optical coherence tomography (OCT) and optical coherence tomography angiography (OCTA) showed the loss of ellipsoid zone and the reduction in paracaval vessel density of both eyes. Genetic analysis identified a novel homozygote of c.8483_8486del (p.Ser2828*) in USH2A. The mutation of premature translation termination causes the deletion of 19 fibronectin type 3 domains(FN3), transmembrane region (TM) and PDZ-binding motif domain, which plays an important role in protein binding. Combining the clinical manifestation and genetic result, the patient was diagnosed with Usher syndrome type 2.
Conclusions: We found a novel mutation of c.8483_8486del in USH2A gene through TES techniques. The result broaden the spectrum of mutations in Usher syndrome type 2 and suggest that the combination of clinical information and TES molecular diagnosis could help USH patients obtain a better diagnosis.
Usher syndrome (USH) is an autosomal recessive disease that is characterized by retinal pigmentosa (RP), sensorineural hearing impairment and vestibule dysfunction. The prevalence is approximately 3.2 to 6.2 per 100,000 [1-4]. It is clinical and genetic heterogeneity. To date, 19 genes and loci are related to USH and atypical USH (RetNet [https://sph.uth.edu/retnet]; October 2019). Among them, 16 genes were identified as causative genes. USH could be divided into three types according to the age of onset, the severity of visual and hearing impairment and vestibule dysfunction. However, due to the point that genetic manifestation is currently not well understood, the rate of missed diagnosis (4%) is high in Asia, especially in China[5].
In patients with USH1, deafness occurs early, and their abnormal visual function is easily ignored. In patients with USH2 and USH3, visual function and hearing abnormalities are gradually progressive. Accurate clinical and molecular diagnoses are the basis of prognosis, treatment selection and eugenics. TES provides us with a new opportunity of revealing the genetic defects in usher syndrome patients[6]. Here, we screened 381 inherited retinal disease (IRD) related genes in a USH2 family and identified a novel mutation of c.8483_8486del (p.Ser2828*) in the USH2A gene.
Case presentation
A 23-year-old man visited our clinic and suffered from deafness from childhood with occasional dizziness, nyctalopia for 10 years and visual acuity declined for nearly 3 years of both eyes. Previously, he was diagnosed with "sensorineural deafness" by an otolarynologist. The patient and his family members gave informed consent for the study, which was approved by the Ethics Committee of Tianjin Medical University Eye Hospital (Tianjin, China). Then, peripheral venous blood samples were collected for NGS or Sanger sequencing. For clinical diagnosis, we performed a comprehensive ocular examination that contained best-corrected visual acuity (BCVA), slit-lamp examination, visual field tests, OCT, OCTA, ultra-wide field fundus photography, and fundus autofluorescence (FAF).
DNA Library Preparation
Genomic DNA was extracted from peripheral blood leukocytes of the patient and his family members using the DNA Extraction Kit (TIANGEN, Beijing, China) following the manufacturer’s instructions. The DNA was quantified with a Nanodrop 2000 (Thermal Fisher Scientific, DE). A minimum of 3 µg of DNA was used for the indexed Illumina libraries according to the manufacturer’s protocol (My Genostics, Inc., Beijing, China). DNA fragments with sizes ranging from 350 bp to 450 bp and those including the adaptor sequences were selected for the DNA libraries.
Targeted Gene Capture and Sequencing
Next, 381 known genes associated with IRD, including USH (Table 1), were selected by a gene capture strategy using the GenCapCustom Enrichment Kit (My Genostics Inc., Beijing, China) following the manufacturer’s protocol. The biotinylated capture probes were designed to tile all of the exons with non-repeated regions. Sequencing was performed on an Illumina HiSeq 2000 sequencer (Illumina, San Diego, CA, USA) for paired-end reads of 150 bp.
Bioinformatics Analysis
Following sequencing, raw image files were processed using Bcl2Fastq software (Bcl2Fastq, Illumina, Inc.) for base calling and raw data generation. Low-quality variations (score≥20) were filtered out. The clean reads were then aligned to the reference human genome using Short Oligonucleotide Analysis Package (SOAP) aligner software (SOAP2.21; soap.genomics.org.cn/soapsnp.html) (hg19). After removing polymerase chain reaction (PCR) duplicates using the Picard program, single nucleotide polymorphisms (SNPs) were determined using the SOAP SNP program, and the deletions and insertions (InDels) were detected using Genome Analysis Toolkit software 3.7. Subsequently, we annotated the identified SNPs and InDels with the Exome-assistant program (http://122.228.158.106/exomeassistant)and viewed the short read alignment using MagicViewer to confirm the candidate SNPs and InDels. Non-synonymous variants were evaluated for pathogenicity using Sorting Intolerant From Tolerant [SIFT; (http://sift.jcvi.org/)] and PolyPhen (http://genetics.bwh.harvard.edu/pph2/). Protein Analysis Through Evolutionary Relationships (PANTHER; www.pantherdb.org) and Pathogenic Mutation Prediction (Pmut; http://mmb.pcb.ub.es/PMut/) were also used.
Expanded Validation and Protein Function Prediction
Genomic DNA of the proband was subjected to TES. Filtered candidate variants identified by anIllumina HiSeq 2000sequencer were confirmed by Sanger sequencing. The coding exons containing the detected mutations were amplified using Ex Tag DNA polymerase (Takara, Dalian). The purified PCR samples were sequenced using the ABI PRISM 3730 genetic analyser (Applied Biosystems; Thermo Fisher Scientific, Inc.), and then sequence traces were analysed with the Mutation Surveyor (Softgenetics, PA). The mutation in the family members was confirmed by the same procedure. Multiple sequence alignments were performed using ClustalW2 with the default setting (http://www.ebi.ac.uk/Tools/clustalw2/).Protein structures were determined by SMART (http://smart.emblheidelberg.de). The 3D structure of the protein variation caused by gene mutation was analysed using Protein Data Bank(PDB) and the homology modelling software Swiss-Model. We collected all genomic DNA samples upon informed consent.
Clinical Findings
A 23-year-old man presented a 10-year history of deafness and poor night vision. His best corrected visual acuity (BCVA) was 0.6/1.0 (R/L). His parents had a consanguineous marriage, and his grandparents had passed away, but they were healthy according to their past medical history and eye conditions. In addition, his parents and his sister were unaffected (Fig. 1). To better clarify the proband's condition, ophthalmologic investigations were performed. Slit-lamp examination showed that the anterior segment of the eyes was normal. The fundus was typical of RP: the appearance of the wax yellow-like disc, the large amount of osteoblast-like pigmentation, and tapering of the retinal vessels, of which the arteries were obvious (Fig. 2a). An abnormal parafoveal ring of increased autofluorescence of ultra-wide-angle images was observed, and there was a ring-like hypoautofluorescence region around the macular and optic disc on FAF imaging (Fig. 2b). We also found decreased retinal thickness and absence of the ellipsoid zone in macular (Fig. 2c). Macular OCTA revealed an enlarging foveal avascular area (FAZ) in the superficial capillary plexus and deep capillary plexus, while macular vascular flow density was also decreased (Fig. 3). Examination with Octopus perimeter device put up a tubular visual field in both eyes (Fig. 4).
Genetic and Molecular Analysis
DNA extracted from the peripheral blood was subjected to TES. Genetic tests showed that the patient had a novel mutation (c.8483_8486del) in the USH2A gene. Moreover, DNA samples extracted from proband’s sister and parents were used for Sanger sequencing. The results confirmed there was a pedigree genetic co-segregation in this family (Fig. 5a). A model structure for USH2A was generated from homology modelling (Fig. 5b). The mutation resulted in premature translation termination, and the stop-gain variant was predicted to remove 2375 amino acids from the encoded proteins, which would result in truncation of the α-βhydrolase domain (Fig. 5c). This may change the overall function of the folded state of the protein (Fig. 5d).
We reported the case of a 23-year-old patient who presented a series of typical clinical features with a novel homozygous mutation, p.Ser2828* (rs1199684717) in USH2A, a gene responsible for USH2 (OMIM:276901). The frequency of the mutation is 0.000004 in the Genome Aggregation Database and it’s found at heterozygous state in one European non Finnish individual. Mutations in USH2A are associated with USH2, it is responsible for almost 50% of USH cases[7].
USH2A codes two alternatively spliced isoforms of usherin. Short ~170 kDa isoform a, consisting of 21 exons, is regarded as an extracellular protein. Full-length ~580 kDa isoform b is a complex transmembrane protein composed of three regions: a large extracellular region consisting of an N-terminal signal peptide, laminin G-like domain (LamGL), laminin domain N-terminal (LamNT), laminin-type EGF-like modules (EGF-Lam), fibronectin type III (FN3) repeats, laminin G domains (LamG); a transmembrane region(TM); and a cytoplasmic C-terminal domain containing a PDZ-binding motif[8, 9]. Usherin is distributed in the periciliary membrane complex and synapse in photoreceptors. All USH1 and USH2 proteins are organized as protein networks by the scaffold proteins harmonin (USH1C), whirlin (USH2D) and SANS (USH1G). Usherin(USH2A) and VLGR1b(USH2C) are part of the links that are intracellularly attached to the scaffold proteins. On the other hand, during the differentiation of the hair bundle, both USH1 and USH2 proteins contribute to the formation of side links located at the tip and the base of the stereocilia, respectively. They exist in multiprotein complexes that work together as molecular networks to anchor them to the stereocilia actin filaments[10-14].
The homozygous mutation (p.Ser2828*) in USH2A made the premature termination translation, as a result, 19 FN3, TM and PDZ-binding motif domains were deleted. FN3 plays a key role in cell adhesion, cell morphology, thrombosis, cell migration, and embryonic differentiation and pathophysiologic processes such as angiogenesis and vascular remodeling[15]. TM exists at the base of the differentiating stereocilia, and it makes up the mechanosensitive hair bundles receptive to sound. PDZ-binding motif domains provide the anchoring of interstereocilia lateral links to the F-acin core of stereocilia[16]. In this regard, we suppose that the absence these domains corresponding to the incompleteness of usherin, which might probably have in turn affected the process of stereocilia differentiation and maturation, resulting in a milder stereocilia dysmorphic phenotype. Several positions are found associated exclusively with pathogenic of FN3 in usherin[17], which support our hypothesis. However, the pathway needs to be confirmed by molecular experiments in the future.
Whole-genome sequencing (WGS), whole-exome sequencing (WES) and TES is three major methodologies for molecular diagnosis of IRD. WGS is useful for detecting copy number variations and structural variations[18]. WES is especially useful for identifying novel IRD related genes. TES is an accurate, rapid and cost-effective approach for screening of multiple genes[19], but still have some major limitations, such as detecting variants in low-depth regions and copy number variations[18, 20]. Because of the higher cost, both of the WGS and WES is less widely used than TES. TES is suitable for molecular diagnosis of USH. As the great diversity of various types of pathogenic genes and the frequent occurrence of new mutations, array-based diagnosis often could not accurately reflect the pathogenicity. USH pathogenic genes have many subtypes and numerous exons. At present, more than 400 coding exons have been commented.[21] Therefore, a higher diagnosis rate can be obtained by using sequence-based diagnosis method.
Here, we report a novel homozygous mutation, c.8483_8486del, in the USH2A gene through TES techniques. The mutation truncated the translation of the USH2A gene, and 19 FN3, TM and PDZ-binding motif domain were lost, which influenced the function of stereocilia. We broadened the spectrum of mutations in the disease and provided a new locus for gene therapy of the disease. The combination of TES molecular diagnosis and clinical information can help USH patients obtain more accurate diagnoses.
Ethics approval and consent to participate
This study adhered to the tenets of the Declaration of Helsinki and was approved by institution review board of Tianjin Medical University Eye Hospital(Tianjin, China), Number:2019KY-02. The
written informed consent to participate was obtained from the patients.
Consent to publication
Written informed consent for publication of their clinical images and genetic results in this manuscript was obtained from the patients according to Ethic Committee Regulations.
Availability of data and materials
All data are fully available without restriction.
Competing interests
The authors declare that they have no competing interests.
Fundings
This study was supported by Natural Science Foundation of Tianjin City (Award number: 17JCYBJC27200; Recipient: Zhiqing Li) and Tianjin Clinical Key Subject Construction Project (Award number: TJLCZDXKQ022; Recipient: Dongjun Xing) in China.
Author Contributions
XL and ZL conceived the idea and take responsibility for the integrity of the data. DX and HZ collected the samples, performed data analyses and wrote the manuscript. Both of authors also contributed equally to this work. RY, LW and LH performed the experiments. All authors have read and approved the final manuscript.
Acknowledgements
We thank all the family members for their participation in this study.
USH: Usher syndrome; OCT: Optical coherence tomography; OCTA: Optical coherence tomography angiography; IRD: inherited retinal disease; FN3: Fibronectin type 3 domains; TES: Targeted exome sequencing; RP: retinal pigmentosa; BCVA: best-corrected visual acuity; FAF: fundus autofluorescence; SOAP: Short Oligonucleotide Analysis Package; SNPs: single nucleotide polymorphisms; PCR: polymerase chain reaction; PDB: Protein Data Bank; FAZ: foveal avascular area; LamGL: laminin G-like domain; LamNT: laminin domain N-terminal; EGF-Lam: laminin-type EGF-like modules; LamG: laminin G domains
Table 1 List of 381 genes chosen for targeted gene capture
ABCA4 |
ADAMTSL4 |
ALMS1 |
ARMS2 |
BBS12 |
BLOC1S6 |
C5AR2 |
CAPN5 |
CEP290 |
CFI |
CLRN1 |
COL11A1 |
CRX |
CYP4V2 |
DTNBP1 |
ERCC8 |
FSCN2 |
GNB3 |
GUCA1B |
HK1 |
HTRA1 |
IMPG2 |
KIAA1549 |
LRIT3 |
MAPKAPK3 |
MKS1 |
NDP |
NRL |
OPN1MW2 |
PCDH15 |
PDE6H |
PEX16 |
PEX7 |
PNPLA6 |
PRPF31 |
RABGGTA |
RD3 |
RIMS1 |
RPGR |
SEMA4A |
SLC25A15 |
SPATA7 |
TENM3 |
TMEM231 |
TRPM1 |
TUBGCP4 |
USH2A |
YAP1 |
ABCB6 |
ADGRA3 |
ANO5 |
ASRGL1 |
BBS2 |
BMP4 |
C5orf42 |
CC2D2A |
CEP41 |
CHM |
CLUAP1 |
COL11A2 |
CSPP1 |
DHCR7 |
EFEMP1 |
EXOSC2 |
FXN |
GNPTG |
GUCY2D |
HMCN1 |
IDH3B |
INPP5E |
KIF11 |
LRP5 |
MC1R |
MLPH |
NEK2 |
NYX |
OPN1SW |
PCYT1A |
PDZD7 |
PEX19 |
PGK1 |
POC1B |
PRPF4 |
RABGGTB |
RDH11 |
RLBP1 |
RPGRIP1 |
SGCD |
SLC26A4 |
SPP2 |
TIMP3 |
TMEM237 |
TSPAN12 |
TUBGCP6 |
VCAN |
ZNF408 |
ABCC6 |
ADGRV1 |
AP3B1 |
ATF6 |
BBS4 |
C10orf11 |
C8orf37 |
CCDC28B |
CEP78 |
CIB2 |
CNGA1 |
COL18A1 |
CST3 |
DHDDS |
ELOVL4 |
EYA1 |
FZD4 |
GPR143 |
GUSB |
HMX1 |
IDUA |
INVS |
KIF7 |
LRPAP1 |
MCOLN1 |
MPDZ |
NMNAT1 |
OAT |
OR2W3 |
PDCD2 |
PEX1 |
PEX2 |
PGR |
POMGNT1 |
PRPF6 |
RAX |
RDH12 |
ROM1 |
RPGRIP1L |
SHH |
SLC38A8 |
STRA6 |
TINF2 |
TMEM67 |
TTC21B |
TULP1 |
VHL |
ZNF423 |
ABHD12 |
ADIPOR1 |
APOE |
ATOH7 |
BBS5 |
C1QTNF5 |
C9 |
CDH23 |
CERKL |
CLDN19 |
CNGA3 |
COL2A1 |
CTC1 |
DHX38 |
ERCC2 |
EYS |
GDF3 |
GPR179 |
HARS |
HPS1 |
IFT140 |
IQCB1 |
KIZ |
LYST |
MERTK |
MTHFR |
NPHP1 |
OCA2 |
OTX2 |
PDE6A |
PEX10 |
PEX26 |
PHYH |
PPT1 |
PRPF8 |
RAX2 |
RDH5 |
RP1 |
RS1 |
SHOX |
SLC39A5 |
TBK1 |
TLR3 |
TMEM98 |
TTC8 |
TYR |
VPS13B |
ZNF513 |
ACBD5 |
AGBL5 |
ARL13B |
ATXN7 |
BBS7 |
C2 |
CA4 |
CDH3 |
CFB |
CLN3 |
CNGB1 |
COL4A1 |
CTNNA1 |
DMD |
ERCC3 |
FAM161A |
GDF6 |
GRK1 |
HEXA |
HPS3 |
IFT172 |
ITGA2B |
KLHL7 |
LZTFL1 |
MFRP |
MTTP |
NPHP3 |
OFD1 |
P3H2 |
PDE6B |
PEX11B |
PEX3 |
PITPNM3 |
PRCD |
PRPH2 |
RB1 |
RGR |
RP1L1 |
SAG |
SIX5 |
SLC45A2 |
TCTN1 |
TLR4 |
TOPORS |
TTLL5 |
TYRP1 |
VSX2 |
ZNF644 |
ACO2 |
AHI1 |
ARL2BP |
BBIP1 |
BBS9 |
C21orf2 |
CABP4 |
CDHR1 |
CFH |
CLN5 |
CNGB3 |
COL9A1 |
CTSD |
DNAJC5 |
ERCC4 |
FBLN5 |
GMPPB |
GRM6 |
HEXB |
HPS4 |
IFT27 |
ITGB3 |
LAMA1 |
MAK |
MFSD8 |
MVK |
NPHP4 |
OPA3 |
PANK2 |
PDE6C |
PEX12 |
PEX5 |
PLA2G5 |
PRDM13 |
PRSS56 |
RBP3 |
RGS9 |
RP2 |
SALL2 |
SIX6 |
SLC7A14 |
TCTN2 |
TMEM126A |
TPP1 |
TTPA |
UNC119 |
WDPCP |
ADAM9 |
AIPL1 |
ARL3 |
BBS1 |
BEST1 |
C2orf71 |
CACNA1F |
CEP164 |
CFHR1 |
CLN6 |
CNNM4 |
COL9A2 |
CTSF |
DRAM2 |
ERCC5 |
FBN2 |
GNAT1 |
GRN |
HFE |
HPS5 |
IMPDH1 |
ITM2B |
LCA5 |
MAN2B1 |
MITF |
MYO5A |
NR2E1 |
OPN1LW |
PAX2 |
PDE6D |
PEX13 |
PEX5L |
PLEKHA1 |
PROM1 |
RAB27A |
RBP4 |
RGS9BP |
RP9 |
SCO2 |
SLC24A1 |
SNRNP200 |
TCTN3 |
TMEM138 |
TRIM32 |
TTR |
USH1C |
WDR19 |
ADAMTS18 |
ALDH1A3 |
ARL6 |
BBS10 |
BLOC1S3 |
C3 |
CACNA2D4 |
CEP250 |
CFHR3 |
CLN8 |
CNOT9 |
CRB1 |
CX3CR1 |
DTHD1 |
ERCC6 |
FLVCR1 |
GNAT2 |
GUCA1A |
HGSNAT |
HPS6 |
IMPG1 |
KCNV2 |
LRAT |
MANBA |
MKKS |
MYO7A |
NR2E3 |
OPN1MW |
PAX6 |
PDE6G |
PEX14 |
PEX6 |
PLK4 |
PRPF3 |
RAB28 |
RCBTB1 |
RHO |
RPE65 |
SDCCAG8 |
SLC24A5 |
SOD1 |
TEAD1 |
TMEM216 |
TRNT1 |
TUB |
USH1G |
WHRN |
- |
- |
- |