Next Article in Journal
Multisystemic Manifestations in Rare Diseases: The Experience of Dyskeratosis Congenita
Next Article in Special Issue
Clinical Features and Genetic Findings of Autosomal Recessive Bestrophinopathy
Previous Article in Journal
Comparative Study of Classification Algorithms for Various DNA Microarray Data
Previous Article in Special Issue
Maternal Uniparental Isodisomy of Chromosome 4 and 8 in Patients with Retinal Dystrophy: SRD5A3-Congenital Disorders of Glycosylation and RP1-Related Retinitis Pigmentosa
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Survey of Multigenic Protein-Altering Variant Frequency in Familial Exudative Vitreo-Retinopathy (FEVR) Patients by Targeted Sequencing of Seven FEVR-Linked Genes

1
Eye Research Institute, Rochester, MI 48309, USA
2
Oakland University William Beaumont School of Medicine, Rochester, MI 48309, USA
3
Associated Retinal Consultants LLC, Royal Oak, MI 48073, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Genes 2022, 13(3), 495; https://doi.org/10.3390/genes13030495
Submission received: 16 June 2021 / Revised: 4 March 2022 / Accepted: 9 March 2022 / Published: 11 March 2022
(This article belongs to the Special Issue Study of Inherited Retinal Diseases)

Abstract

:
While Inherited Retinal Diseases (IRDs) are typically considered rare diseases, Familial Exudative Vitreo-Retinopathy (FEVR) and Norrie Disease (ND) are more rare than retinitis pigmentosa. We wanted to determine if multigenic protein-altering variants are common in FEVR subjects within a set of FEVR-related genes. The potential occurrence of protein-altering variants in two different genes has been documented in a very small percentage of patients, but potential multigenic contributions to FEVR remain unclear. Genes involved in these orphan pediatric retinal diseases are not universally included in available IRD targeted-sequencing panels, and cost is also a factor limiting multigenic-sequence-based testing for these rare conditions. To provide an accurate solution at lower cost, we developed a targeted-sequencing protocol that includes seven genes involved in Familial Exudative Vitreo-Retinopathy (FEVR) and Norrie disease. Seventy-six DNA samples from persons refered to clinic with possible FEVR and some close relatives were sequenced using a novel Oakland-ERI orphan pediatric retinal disease panel (version 2) providing 900 times average read coverage. The seven genes involved in FEVR/ND were: NDP (ChrX), CTNNB1 (Chr3); TSPAN12 (Chr7); KIF11 (Chr10), FZD4 (Chr11), LRP5 (Chr11), ZNF408 (Chr11). A total of 33 variants were found that alter protein sequence, with the following relative distribution: LRP5 13/33 (40%), FZD4 9/33 (27%), ZNF408 6/33 (18%), (KIF11 3/33 (9%), NDP 1/33 (3%), CTNNB1 1/33 (3%). Most protein-altering variants, 85%, were found in three genes: FZD4, LRP5, and ZNF408. Four previously known pathogenic variants were detected in five families and two unrelated individuals. Two novel, likely pathogenic variants were detected in one family (FZD4: Cys450ter), and a likely pathogenic frame shift termination variant was detected in one unrelated individual (LRP5: Ala919CysfsTer67). The average number of genes with protein-altering variants was greater in subjects with confirmed FEVR (1.46, n = 30) compared to subjects confirmed unaffected by FEVR (0.95, n = 20), (p = 0.009). Thirty-four percent of persons sequenced had digenic and trigenic protein-altering variants within this set of FEVR genes, which was much greater than expected in the general population (3.6%), as derived from GnomAD data. While the potential contributions to FEVR are not known for most of the variants in a multigenic context, the high multigenic frequency suggests that potential multigenic contributions to FEVR severity warrant future investigation. The targeted-sequencing format developed will support such exploration by reducing the testing cost to $250 (US) for seven genes and facilitating greater access to genetic testing for families with this very rare inherited retinal disease.

1. Introduction

FEVR (Familial Exudative Vitreo-retinopathy) and Norrie disease are inherited disorders that affect the development of the neural retinal vasculature, resulting in avascular regions in the peripheral retina [1,2,3]. FEVR was first described in 1968 and can result in blindness from retinal traction, folding, detachments, neovasuclarization, and vitreous hemorrhage [4]. FEVR phenotypically can show large variations in severity, even within a single family, and Shastry and Trese (2004) were the first to suggest that digenic alleles could exist in some FEVR families as a factor contributing to the variable penetrance [5]. More recently, Li et al. [6] reported evidence that digenic protein-altering variants occurred in 2.7% of their study cohort (13/487) using a four-gene panel and suggested that increasing the number of FEVR-linked genes sequenced might increase this percentage [6]. For this study, our goal was to increase the number of genes sequenced to seven genes. Currently, variants assocated with FEVR impact at least seven proteins (genes): Norrin Cystine-Knot Growth Factor NDP (NDP), Catenin β-1 (CTNNB1), Frizzled-4 (FZD4), Kinesin Family Member-11 (KIF11), LDL-Receptor-Related Protein-5 (LRP5), Tetraspanin-12 (TSPAN12), and Zinc-Finger Protein-408 (ZNF408). Due to the rarity of some of these conditions, DNA testing is still difficult to access, expensive, not available in most countries, and not supported by medical insurance in the United States.
We first pre-tested an amplicon targeted-sequencing panel with six FEVR-related genes and one for retinoschisis, and we confirmed that it detected known variants from 15 subjects that we had analyzed in years past using traditional PCR and Sanger sequencing. After that preliminary feasiblity test, we settled on using the Illumina AmpliSeq targeted-sequencing process to develop our current panel (version 2) for eight genes involved in three very rare pediatric retinal conditions: FEVR, Norrie disease, and retinoschisis. The seven FEVR/Norrie disease-related genes listed above are included. Using the AmpliSeq workflow for sequencing on the iSeq-100 instrument, the cost of testing was reduced to only $250 per sample ($32 US per gene), providing a practical research-sequencing service for families with these rare IRDs. The goal of this initial study was to ask if a significant percentage of FEVR subjects may have multigenic protein-altering variants in these seven FEVR-relevant genes. We then analyzed a novel group of 76 persons for DNA sequence variants that alter the amino acid sequence of the seven proteins of interest.
All the genes of interest are important to the biological function of retinal vascular endothelial cells. Most of the gene products are members of the Norrin-based signaling pathway specific to endothelial cells of the neural retinal microvasculature. Retinal endothelial cells are reliant on Norrin with their unique combination of receptor proteins that make them sensitive and responsive to Norrin. This combination consists of the co-expression of the FZD-4, TSPAN-12, and LRP-5 proteins (see Figure 1). While different Wnt family members can bind to FZD-4 in other cell types, the combination with LRP-5 and TSPAN-12 is required for the unique and strong response to Norrin [7]. In addition to the five proteins shown in Figure 1, two other proteins that are abundantly expressed and particularly required in retinal endothelial cells are KIF-11 and ZNF 408. ZNF 408 is a zinc-finger transcription factor that is essential for the adequate expression of genes required for normal endothelial cell function, and protein-altering variants of this gene are also reported to cause FEVR [8]. KIF-11 is the predominant Kinesin protein family member expressed in retinal endothelial cells. KIF-11 variants are also associated with several medical conditions including microcephaly, lymphedema, and chorioretinal dysplasia (MLCRD) as well as chorioretinal dysplasia, microcephaly, and mental retardation (CDMMR) [9,10]. Patients with KIF-11 mutations can have one or more of these conditions in combination.
Due to the rarity and complexity of these diseases, diagnosing FEVR can be very challenging. In addition, their polygenic nature makes identifying the cause difficult because of the small population of patients and the limited amount of research available [12]. Understanding the potential role of multiple protein-altering variants in more than one FEVR-relevant gene in the same patient faces two challenges. First, obtaining enough persons with these very rare conditions to consent for analysis is difficult. Second, an economical methodology is required to accurately sequence not just one, two, or four genes but a multigene panel that includes as many FEVR-related genes as possible. The first challenge was solved by the fact that the Associated Retinal Consultants practice (Royal Oak, MI, USA) is a tertiary referral center that sees FEVR patients from all over the United States and the world. The second challenge remaining to be solved was addressed here by the development of a custom AmpliSeq targeted-sequencing panel.
Before we can start to address the question of potential multigenic involvement in FEVR, we must first answer this question: what percentage of FEVR subjects have multiple protein-altering variants in this pool of seven FEVR-related genes? That question is the subject of this report. If multigenic protein-altering variants are not common in this pool of the most frequently involved FEVR-related genes, then it is unlikely we would need to be highly concerned with the potential for multigenic contributions in most FEVR subjects. On the other hand, if a significant proportion of FEVR subjects have multigenic protein-altering variants, then the multiple combinations would warrant some consideration in the analysis of families with this condition.

2. Materials and Methods

2.1. IRB Protocol and Test Group Characteristics

All of the human DNA sequencing work used for this research was accomplished under the approval of the Oakland University Institutional Research Committee through an IRB management agreement for the Creation of an Ophthalmology BioBank (OU-IRB approval #1314454-2) with the Western IRB (WIRB). Subjects providing DNA consented under WIRB-approved protocols for the DNA Eye-bank at Associated Retinal Consultants (ARC, Royal Oak, MI, USA) and for DNA-sequencing.
For this analysis, 76 persons’ samples were processed for targeted DNA-sequencing and analysis. Subjects included persons referred to the ARC clinic for FEVR evaluation and some close relatives of FEVR or suspected FEVR patients. Many patients travelled from other US states or other countries, so an examination of all first-degree relatives by the ARC clinic was often not possible. After ARC clinical evaluation, 35/76 subjects were considered to have a confirmed FEVR phenotype and were graded for FEVR clinical staging as described previously by the ARC group [13,14]. Additionally, 20/76 subjects were clinically confirmed to be unaffected by FEVR, 5/76 subjects were considered to have potential mixed FEVR and ROP (FROP) (Retinopathy of Prematurity), 1/76 was clinically classified as just ROP, and 15/76 subjects were of unknown FEVR status. The latter subjects were mostly relatives of FEVR patients who provided blood samples to the ARC Eye Biobank but who could not be clinically evaluated for FEVR.

2.2. Sample Storage and Handling

Whole-blood samples were collected at the clinical site and frozen at −70 °C. Frozen aliquots (100–200 µL) of blood samples from consenting persons were made in 1.5 mL sterile microfuge tubes and kept frozen until use. All blood samples obtained from subjects were de-identified and analyzed with a randomly assigned research ID number prior to transportation to the Oakland University Eye Research Institute (ERI) with lists of sample IDs. The research side of the project (ERI) did not have access to any key for decoding sample identification. Family relationships between de-identified samples and their gender were known for analysis purposes only. Samples were stored at −70 °C in the ERI until extraction of genomic DNA.

2.3. Extraction of Genomic DNA

Genomic DNA was extracted from 100–200 µL of frozen whole blood with the Thermo-Fisher Pure-Link Genomic DNA Extraction Kit (ThermoFisher, Waltham, MA, USA). DNA concentrations were measured with a Qubit Fluorimeter, from 1 µL of genomic DNA extracted, with the Qubit dsDNA HS Assay Kit (ThermoFisher, Waltham, MA, USA).

2.4. AmpliSeq Targeted Library Design In Silico

The targeted AmpiSeq panel was designed in silico using Illumina’s DesignStudio (https://designstudio.illumina.com/, accessed on 12 January 2022). Parameters included eight target genes, with seven of the genes related to FEVR/ND and another very rare pediatric disease gene (RS1) for retinoschisis: NDP (ChrX), RS1 (ChrX); CTNNB1 (Chr3); TSPAN12 (Chr7); KIF11 (Chr10), FZD4 (Chr11), LRP5 (Chr11), ZNF408 (Chr11). Options were selected for an average 250 base pair amplicon size, targeting exons with amplicon overlap, and at least 25 base pairs of 5-prime and 3-prime adjacent intron sequence. Designs of different pool number were compared in silico, and 3 pools were selected from predictions of >99.95% base coverage. The final Amplicon targeted panel numbered 360 PCR-primers to generate 180 amplicons, covering just over 32,000 base pairs in a sequence for 83 exons.

2.5. AmpliSeq Targeted Library Preparation

Illumina’s AmpliSeq Library preparation protocol was followed as per the manufacturer’s instructions. (Illumina, Sand Diego, CA, USA) Other Illumina reagents and kits were used per the Illumina workflow. The AmpliSeq Library Plus kit was used for 16 samples per kit using a three-pool amplicon set (beginning with three PCRs before combining into one PCR product per patient). The AmpliSeq CD Indexes Set B for Illumina (96) was used in the workflow to provide unique bar-coded end indexes for each patient library. These provided for a different index for each strand to enable separate sequencing for each paired amplicon strand.
Final targeted amplicon libraries were quality-checked using an Agilent Bioanalyzer and the Agilent DNA 1000 kit. After a few groups of libraries were found to always be of sufficient quality for sequencing, libraries could also simply be quantified using the Qubit DNA assay to establish library DNA concentration. In total, 16 to 48 libraries were diluted to equal concentrations according to Illumina’s protocol, and a 5% non-human PhiX DNA spike-in standard (Illumina) was included for monitoring of sequencing error rate and quality by the iSeq-100 system. Final total pooled library concentration was 50 pM, divided equally among samples (up to 50, of which 20 µL was loaded into a single-use iSeq-100 Reagent V2 cartridge and flow cell). Typical run times started in the afternoon, continuing overnight, and took 18–19 h from instrument start to the delivery of final variant call file data for all samples.

2.6. Sequencing on the iSeq-100 System

For sequencing, up to 48 patient sequencing libraries were pooled for a single instrument run. The on-board software workflow provided by the iSeq-100 included: Illumina’s DNA Amplicon Workflow (3.24.1.8+master), BWA-MEM Whole-Genome (Aligner, 0.7.9a-isis-1.0.2), Pisces Variant Caller (5.2.9.23), Illumina Annotation Engine (2.0.11-0-g7fb24a09), Bam Metrics (v.0.0.22), and SAMtools (0.1.19-isis-1.0.3). The annotation source option was selected as RefSeq, and the dataset was version 91.26.44. Options for sequencing included paired end sequencing format, with a minimum threshold read depth of 10 for final base calling. Reference genome selection was Homo Sapiens NCBI GRCh38 with decoys.

2.7. Bioinformatics Data and Analysis

The data from sequencing the eight genes involved in Norrie disease, FEVR, and retinoschisis typically provided about 20–30 DNA sequence variants per person using the reference human genome sequence (GRCh38 - hg38). Variant impacts on protein sequences and allele frequency data were determined from the ClinVar database (https://www.ncbi.nlm.nih.gov/clinvar/), the Genome Aggregation Database, gnomAD https://gnomad.broadinstitute.org/, and the Variant Effect Predictor (VEP) database at ENSEMBL (https://uswest.ensembl.org/info/docs/tools/vep/index.html). Databases were last accessed on 12 January 2022. Bar graphs and statistical testing were completed using software tools in R, managed in the RStudio Desktop open-source environment (Software available from https://www.rstudio.com).

3. Results

3.1. Sequencing Data Quality and Performance

Table 1 contains the sequencing quality metrics obtained in a multiple-sample run using the designed AmpliSeq targeted library. Overall, 95.5% of base reads were >Q30 quality (error < 1/1000). The percentage of on-target bases passing filter (PF) was 92.2%. Average sequencing depth of coverage was 978.

3.2. Total Variant Detection Results

Numerous variants were detected from the 76 persons tested. Table 2 contains the counts for SNVs (single-nucleotide variants), insertions, and deletions according to variant category and location. Average numbers of variants per sample were: 19.5 SNVs, 1.7 inserts, and 1.5 deletions. A total of 937 SNVs were detected, of which 3 SNVs were not found in the dbSNP database. Of a total of 81 insertions, 3 were not found in dbSNP, and of 71 deletions, 10 were not present in dbSNP.

3.3. Protein-Altering Variants in Cohort of 76 FEVR Subjects and Relatives

Persons included in this round of sequencing were refered to Associated Retinal Consultants (ARC), and age, diagnosis, and any family relationships at presentation to the ARC clinic are listed in Table 3. The sample set included related persons from 18 families grouped first in the table, followed by individual patients without relatives. (Note: there is no Family 6.) FEVR stages (base-grade) are provided in Table 3 if stages were confirmed by examination at the ARC clinic. The base stages were: (1) avascular periphery, (2) retinal neovascularization present, (3) extramacular retinal detachment, (4) macula-involving retinal detachment, and (5) total retinal detachment. The final distribution of subjects grouped into five categories is shown in Figure 2: FEVR, U (Unaffected by FEVR), ROP (Retinopathy of Prematurity), ROP with possible FEVR complication (FROP), and NA (unknown, Not Available).
Of the entire cohort sequenced, 50 of the subjects were present with a close relative, representing 18 family groupings. Other subjects in the cohort were included as unrelated individuals who were referred to the ARC clinic for potential diagnosis of FEVR. Results for the family groupings are summarized below. Note: there was no Family 6 designation in this study.
Thirty-three (33) protein-altering variants were detected in the entire group of sequenced patients. These variants with their HGVS descriptions are listed in Table 4, as well as any known fractional allele frequencies (AF). Of the 33 variants, 85% were found in just three of the seven genes: FZD4, LRP5, and ZNF408. Figure 3 contains a chart of the relative variant-to-gene distributions: LRP5 13/33 (40%), FZD4 9/33 (27%), ZNF408 6/33 (18%), KIF11 3/33 (9%), NDP 1/33 (3%), CTNNB1 1/33 (3%). The impacts on protein sequences were: 27 single amino substitutions (81%), 3 amino acid deletions (9%), 1 stop-gained alteration (3%), 1 amino acid insertion duplication (3%), and 1 splice donor site loss (3%).
Two novel, likely pathogenic variants were detected in one family and one unrelated individual based upon their similarity to known pathogenic variants that truncate significant portions of FZD4 and LRP5. FZD4: Cys450ter (Family 11) and LRP5: Ala919CysfsTer67 (Subject 59). Four previously known pathogenic protein-altering variants were found in five families and two unrelated individuals. They included: FZD4: Met105Val (Family 7, 9, 11), LRP5: c.4488+2T>G (splice donor) (Family 10, 14), LRP5: Cys913LeufsTer73 (Subject 58), and NDP: His42Arg (Subject 56) (Table 3).
Family 1: Affected twins both had LRP5: Val667Met and LRP5: Ala1330Val with high FEVR grades. The LRP5: Ala1330Val Allele Frequency (AF) is 13.4%, currently categorized in ClinVar as likely benign. LRP5: Val667Met is also considered benign (AF = 13.4%) Interestingly, Proband-1 has more severe disease (grade 5) than their twin (grade 4) and the LRP5: Ala1330Val homozygous variant, while their twin is heterozygous for both variants. The visual status of the mother, maternal aunt, and maternal uncle were not known.
Family 2: Both the proband son and the mother had lower FEVR grades, but they had different protein-altering variants. The proband’s variant CTTNB1: Val273Met has a very low Allele Frequency (AF = 0.0021%), was not found in ClinVar and would be of uncertain significance at this time. The mother had digenic variants that are listed as benign in ClinVar, LRP5: Leu20dup and ZNF408: Arg337Pro.
Family 3: The LRP5: Ala1330Val homozygous variant was combined with the LRP5: Leu20dup benign variant in proband-3, an adult son with a FEVR grade of 4. The retinal status of the proband’s father was unknown, but the father was heterozygous for he LRP5: Ala1330Val variant.
Family 4: Sequencing was only possible for the the parents of an affected, not-yet-sequenced infant, who were themselves unaffected by FEVR. Only the father had a combination of two currently likely benign and benign variants: LRP: Met1086Val and ZNF408: Arg337Pro.
Family 5: The proband infant daughter (FEVR grade 3) shared the same two individually benign variants in LRP5 and ZNF408 as her unaffected mother. The father had a FEVR base-grade 1 and the same ZNF408 variant but two completely different likely benign and benign variants of LRP5.
Family 7: A known pathogenic variant FZD4: Met105Val was present in the proband teenage daughter with a FEVR base-grade 2. Clinical diagnosis of the father, mother, and one half-sibling was not available. The father had the same FZD4: Met105Val pathogenic variant, while the mother had no detected protein changes. The related sister did not have the pathogenic variant, only one benign variant (ZNF408: Val194_Val197del) inherited from the father.
Family 8: The proband teenage daughter diagnosed with FEVR grade 1 and two diagnosed unaffected pediatric siblings were sequenced. The proband had two variants. One, KIF11: Glu129Ala, of low allele frequency < 0.0001, not found in ClinVar, we considered to be of uncertain significance. The second variant present was homozygous, LRP5: Ala1330Val, a variant categorized as benign in ClinVar. The LRP5: Ala1330Val variant was present and heterozygous in an unaffected brother. An unaffected sister had the relatively common and benign variant ZNF408: Val194_Val197del (AF = 13.4%).
Family 9: The proband pediatric son had the known pathogenic variant FZD4: Met105Val (AF = 0.0043%), correlating with a severe FEVR base-grade of 5. Sequencing was completed for the unaffected father, who was found to have no protein changes from the panel. Presumably the pathogenic variant was inherited maternally.
Family 10: The proband was a young adult female, 18 years of age at presentation to our clinical group, with a FEVR base-grade of 5. A splice donor variant in LRP5: c.4488+2T>G considered to be pathogenic was detected in the proband and in the mother and a maternal uncle. Both the mother and uncle presented with a FEVR grade of 1. This variant was not inherited by the uncle’s son, who is unaffected and was found to have no protein changes from the panel. The mother had one additional variant, LRP5: Leu20dup (AF = 10.1%), considered benign. The proband also had a second variant, FZD4: Pro11Gln (AF = 0.0021%), not in ClinVar that we currently consider of uncertain significance.
Family 11: The proband pediatric daughter had a FEVR grade of 5 at presentation to the clinic and was found to have two known pathogenic variants, FZD4: Cys450ter and FZD4: Met105Val. The first variant FZD4: Cys450ter was presumably inherited from the mother, who had a FEVR grade of 1. The father was not available for sequencing. The mother had a second variant, the benign and common LRP5: Leu20dup, which was also detected in the proband’s maternal grandfather, who was unaffected. The proband’s maternal grandfather also had a second common benign variant, ZNF408: Val194_Val197del (AF = 13.1%).
Family 12: The proband, pediatric daughter had a base FEVR grade of 5 and two LRP5 variants considered benign in ClinVar: LRP: Val667Met and LRP: Ala1330Val. The retinal status of the parents was not known. Both variants appeared to be inherited from the father, who also had a third common variant, ZNF408: Val194_Val197del. The mother had no protein changes from the panel.
Family 13: The proband pediatric son with a FEVR grade of 4 was found to have the same two benign LRP5 variants as in Family 12. LRP: Val667Met and LRP: Ala1330Val. The proband’s mother with a FEVR grade of 1 had no protein-altering variants from the panel. The same two variants were shared by an unaffected sister and were inherited from the unaffected father who also had the common variant ZNF408: Val194_Val197del. The maternal grandmother was unaffected and had no protein changes from the test panel, while the maternal grandfather with a FEVR base grade of 1 had four common benign variants, LRP5: Val667Met, LRP5: Ala1330Val, FZD4: Pro33Ser, and FZD4:Pro168Ser.
Family 14: The proband, pediatric son with FEVR grade 4 was found to have the pathogenic LRP5: c.4488+2T>G splice donor variant. The mother with a FEVR grade of 1 did not have this variant but was found have a common variant, ZNF408: Val194_Val197del (AF = 13%), and a deletion, LRP5: Leu16_Leu20del (AF = 0.0007%), which was not in ClinVar and we currently categorize as of uncertain significance.
Family 15: Currently, two twin boys were sequenced with base FEVR grades of 1 and 2 at the time of presentation. Both had the same variant ZNF408: Glu230Gly, listed in ClinVar as likely benign (AF = 0.037%).
Family 16: The proband was a male pediatric case referred with possible FEVR but subsequently considered unaffected by FEVR. One sister without a clinical work up was also sequenced. Both had two common benign variants, LRP5: Ala1330Val and ZNF408: Val194_Val197del.
Family 17: The proband was an infant with uncertain FEVR, whose sample was not available for sequencing. Both parents of unknown retinal status were sequenced and found to have a common benign deletion variant, ZNF408: Val194_Val197del (AF = 13%).
Family 18: The proband was an infant with severe bilateral FEVR with retinal detachment whose sample was not available for sequencing. The mother is diagnosed as unaffected by FEVR, and no missense variants were detected upon sequencing. The father’s retinal status is unknown, but DNA sequencing detected two common benign variants, LRP5: Val667Met and LRP5: Ala1330Val, and a more rare variant, FZD4: Arg127Cys (AF = 0.0046%), that was listed as likely benign in ClinVar.

3.4. Multiple Variants and Multigenic Variants

Data from all the subjects sequenced revealed the existence of multiple protein variants in single genes, as well as variants in multiple FEVR-related genes in the same individuals. The latter included both digenic and trigenic variant combinations, as summarized in Table 5. Considering just the subset of patients with an ARC-confirmed diagnosis of FEVR or Unaffected (U), the distribution of the number of genes with protein-altering variants is shown in Figure 4.
The average number of genes with protein-altering variants in the group of confirmed FEVR diagnosis was greater than the unaffected group of subjects, averaging 1.46 genes for FEVR and 0.95 genes for Unaffected (p = 0.009, t-test). Of the total subjects sequenced in this cohort, 26/76 were di- or tri-genic for protein-altering variants within the group of 7- FEVR-linked genes sequenced. Using data from the GnomAD database, derived from over 14,000 genome sequences, the occurrence of the same combined frequencies in the general population was estimated to total 3.6%. This was substantially lower than the frequency in our cohort, 34.2% (Chi-squared, p = 3.0 × 10−17).

4. Discussion

4.1. FEVR Phenotype and Genes

While numerous genes are involved in the FEVR phenotype, there is no simple correlation between disease severity and specific genes. This is most likely because the impact on the neural retina results from an avascular peripheral retina that results in hypoxia and elevated vascular endothelial growth factor. Any genetic variant that inhibits development of the neural retinal vasculature will generate the situation of peripheral hypoxia. Similar to ROP (retinopathy of prematurity), the smaller the final zone of vascularized retina, the more severe the depth of hypoxia and the greater the potential for a more severe FEVR phenotype [1].
The Norrin-Wnt signaling pathway is unique to the neural retinal endothelium, and indeed, active Norrin signaling is essential for normal growth and proliferation of the developing retinal vasculature [7,15]. Likewise, a loss of KIF11′s kinesin activity can impair growth and proliferation of the retinal vasculature and again results in peripheral avascular neural retina [9]. KIF11 FEVR patients tend to have other significant syndromes (MLCRD, CDMMR), none of which were noted in the 76 subjects presented for analysis here. ZNF408 is a multiple zinc-finger DNA-binding protein, and point mutations that impair ZNF408′s DNA-binding activity result in diminished expression of genes required for vascular development [16]. Thus, disrupting the normal function of Norrin-Wnt signaling, KIF11 kinesin action, or ZNF408′s gene activation activity can result in mild to severe FEVR phenotypes.
An AmpliSeq targeted-sequencing panel including seven FEVR-related genes revealed protein-altering variants in most of the subjects that were sequenced for this initial study. The cohort included persons diagnosed with FEVR, close relatives of some proband FEVR patients, and individual persons referred to the ARC clinic for potential FEVR. Of the 33 protein-altering variants detected, 85% were found in three of the genes: FZD4, LRP5, and ZNF408. ZNF408 ranked third in the number of protein-changing variants we found in this study, which was also found in one small study from China [6]. Another study from China of a larger cohort of several hundred FEVR patients found that TSPAN12 variants ranked third after FZD4 and LRP5 [17]. FZD4 and LRP5 variants were the two most frequent in our cohort, which seems to be a distribution shared with previously published studies [6,17,18,19]. Together, these studies indicate some regional variation occurs, but in general, the two most commonly variant genes were FZD4 and LRP5. Most of the variants, 20/33, had low allele frequencies, <0.01%, and are candidate variants to be analyzed further for their potential contribution to FEVR in this group of subjects.

4.2. Frequency of Multigenic Protein-Altering Variants

Our targeted panel covered seven FEVR-related genes with over 900 times average read coverage and >95.5% Q-30 quality reads, which provided us with high confidence in sensitivity for variant detection and for determination of zygosity. These sequencing results have permitted us to answer the question posed in the introduction: what percentage of FEVR subjects may have multiple protein-altering variants in this pool of seven FEVR-related genes? In our cohort, 34% (26/76) of subjects had multiple protein-altering variants in two or three different FEVR-related genes. This was greater than the 3.6% expected from the general population, based on GnomAD genome data. Selections of patients and families based on FEVR phenotypes greatly increased the relative presence of multigenic protein-altering variants in the FEVR-related gene set. While this association is not sufficient proof that multigenic variants contribute to the FEVR phenotype, the high percentage of multigenic variants indicates that we cannot absolutely rule out such potential contributions due to scarcity.
Li et al. [6] reported detection of potential digenic variants in just under 3% of their study cohort when sequencing four FEVR-related genes (FZD4, LRP5, TSPAN12, ZNF408) [6]. A relatively higher incidence of multigenic variants in our study may be due to several factors, including the inclusion of more FEVR-related genes, sequencing to a greater depth for sensitivity, and the fact that ZNF408 variants were relatively highly represented in our cohort. ZNF408 variants occurred in 18 subjects in digenic and trigenic combinations.
Comparing only the confirmed FEVR-positive to confirmed FEVR-negative (unaffected) individuals in our cohort, we found that the average number of FEVR-related genes per person with protein-altering variations was higher in the confirmed FEVR subset. While this also was not sufficient proof to conclude that multigenic contributions are important in FEVR incidence, this relative occurrence would be expected in that scenario. Some of our specific findings have also provided us with guidance for continued follow-up of specific families, where we can extend pedigrees with both sequencing and clinical diagnosis.

4.3. Known and Novel Pathogenic Variants

Pathogenic variants were uncovered in several families (Table 3 and Table 4). FZD4: Met105Val was found in Family 7, 9, and 11. A new likely pathogenic variant, FZD4: Cys450ter, was also found in Family 11. The splice donor variant LRP5: c.4488+2T>G was detected in Family 10 and 14. The other pathogenic variants detected in individual subjects included LRP5: Cys913LeufsTer73 and a new similar frame shift termination variant LRP5: Ala919CysfsTer67 (likely pathogenic). Finally, one known pathogenic variant affecting Norrin itself is NDP: His42Arg.
Family 1 leads us to consider an additional question going forward. Can heterozygous variants currently considered benign on their own contribute to FEVR occurrence when homozygous or when two different variants impact the same gene such as LRP5? That situation may be present here, where both twins had two variants in LRP5: Val 667Met and Ala1330Val but the more severely affected twin inherited the LRP5: Ala1330Val variant homozygous. Interestingly, this homozygous variant LRP5: Ala1330Val also occurred in the Family 3 proband, again in combination with a different LRP5 benign variant. Again, LRP5: Ala1330Val appeared homozygous in the Family 8 proband with a novel KIF11 variant of unknown significance, while the LRP5: Ala1330Val variant was present and heterozygous in an unaffected sibling.
Family 10 was identified for future analysis of a potential digenic contribution to FEVR severity. Both the mother and uncle of proband10 had a base FEVR grade of 1 and shared a pathogenic splice variant in LRP5: c.4488+2T>G. The proband was more severely affected with FEVR (stage 5) and had a second variant of uncertain significance in FZD4: Pro11Gln (AF = 0.0021%). The uncle’s son, without any of these variants, was unaffected by FEVR.
Another example of features that can complicate our understanding of FEVR severity within the same family is the occurrence of completely different variants, with some present in a multigenic manner. Family 14 was identified for future follow-up analysis as an example of this situation. The proband son had a FEVR base grade of 4 and the LRP5 pathogenic splice donor variant, LRP5: c.4488+2T>G, while the proband’s mother has a lower FEVR grade of 1 but did not have this variant. Instead, the mother was digenic for the benign ZNF408: Val194_Val197del (AF = 13%) and LRP5: Leu16_Leu20del (AF = 0.0007%). The latter was not in ClinVar, and we currently categorize this variant as of uncertain significance.

4.4. Sequencing Approach

Targeted sequencing, Whole-Exome Sequencing (WES), and Whole-Genome Sequencing (WGS) have all been employed for genetic analysis of Inherited Retinal Diseases (IRDs). All three have their advantages and disadvantages. A large IRD targeted panel applied to a disease-appropriate cohort of 85 children demonstrated 80% identification of the likely disease-causing variant [20]. Governed by the concept that “perfection is the enemy of progress” (Winston Churchill, speech, 11 October 1952), we designed a customized seven-gene targeted panel, to balance several parameters including the patient consenting process, FEVR-gene sequencing coverage, SNV sensitivity, and cost.
Targeted IRD disease panels exist that include some FEVR-related genes, but they are missing one or more key FEVR-linked genes. While WES can survey most of the known exome, some IRD gene comparison studies have shown that targeted sequencing of IRD genes is generally superior to WES, as the latter may have poor capture efficiency for some exome regions [21,22]. WGS sequencing is theoretically the best method to survey all genes without a predetermined bias, but the ability to improve accuracy may require sequencing to a depth of coverage that greatly reduces the number of patients that can be tested with the available resources. While genome-wide variant data can contain a potentially novel disease-linked variant, a logical functional analysis may not be available to identify which variants from tens of thousands are relevant to the disease process.
Another consideration when adopting the WES and WGS approach concerns the potential detection of non-ocular genetic disease variants. In the United States, the American College of Medical Genetics and Genomics (ACMG) recommends that additional protocols should be included for reporting a pathological variant from the current ACMG list of medically actionable genes [23]. Currently, medical actionable gene reporting is not policy in most countries. Full patient education for ethical consent is even more important, and the choice to participate is more complicated for families considering the potential discovery of serious non-ocular health conditions using genome-wide analysis.

5. Conclusions

From this initial study, we can now plan future analysis of individual probands and additional relatives to determine if any of the family-specific multigenic variant combinations might influence phenotype penetrance in FEVR. The initial question posed by this study was answered: what is the incidence of multigenic protein-altering variants among known FEVR-linked genes in our FEVR cohort? In this study, the answer was 34%, far exceeding the percentage expected in the general population. This high proportion of multigenic protein-altering variants in a specific FEVR-linked gene panel establishes that there is a theoretical potential for multigenic variants to contribute to FEVR’s variable phenotypic presentation. While we cannot point to absolute proof of multigenic contributions to the FEVR phenotype from this initial study, this data set forms the basis for further investigation of this possibility.

Author Contributions

Conceptualization, W.D., K.A.D. and K.P.M.; Data curation, A.P.C., W.D., K.M., K.A.D. and K.P.M.; Formal analysis, A.P.C., W.D., M.S., A.S., D.J., L.J., K.K., M.D., K.S., K.A.D. and K.P.M.; Funding acquisition, K.P.M.; Investigation, W.D., M.T.T., A.C., K.A.D. and K.P.M.; Methodology, W.D., M.S., N.H., K.A.D. and K.P.M.; Project administration, W.D., K.P.M., K.A.D. and K.M.; Resources, J.A.F., A.E.G., K.P.M., M.T.T., A.C. and K.M.; Supervision, W.D., K.P.M. and K.A.D.; Validation, K.P.M., W.D., K.M. and K.A.D.; Visualization, A.P.C. and K.P.M.; Writing—original draft, A.P.C., W.D. and K.P.M.; Writing—review and editing, A.P.C., W.D., M.S., N.H., J.A.F., A.E.G., M.T.T., A.C., K.A.D. and K.P.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by awards to K.P.M. from the Pediatric Retinal Research Foundation, Novi, MI 48375, USA, (2019, 2020, 2021) and The Carls Foundation, Bloomfield Hills, MI 48304, USA (2021).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, through an IRB management agreement for the Creation of an Ophthalmology BioBank (OU-IRB approval # 1314454-2) with the Western IRB (WIRB). Subjects providing DNA consented under WIRB-approved protocols for the DNA Eye-bank at Associated Retinal Consultants (Royal Oak, Michigan) and for DNA-sequencing.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

For more detalied information on the targeted DNA-sequenicng methodology and process used for this study, as well as the integration of this sequencing service into biomedical science and medical student education, please contact Dr. Kenneth P. Mitton ([email protected]).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Drenser, K.A.; Dailey, W.; Vinekar, A.; Dalal, K.; Capone, A.; Trese, M.T. Clinical Presentation and Genetic Correlation of Patients With Mutations Affecting the FZD4 Gene. Arch. Ophthalmol. 2009, 127, 1649–1654. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Tauqeer, Z.; Yonekawa, Y. Familial Exudative Vitreoretinopathy: Pathophysiology, Diagnosis, and Management. Asia Pacific J. Ophthalmol. 2018, 7, 176–182. [Google Scholar] [CrossRef]
  3. Liu, J.; Zhu, J.; Yang, J.; Zhang, X.; Zhang, Q.; Zhao, P. Prenatal diagnosis of familial exudative vitreoretinopathy and Norrie disease. Mol. Genet. Genom. Med. 2018, 7, e00503. [Google Scholar] [CrossRef] [PubMed]
  4. Criswick, V.G.; Schepens, C.L. Familial exudative vitreoretinopathy. Am. J. Ophthalmol. 1969, 68, 578–594. [Google Scholar] [CrossRef]
  5. Shastry, B.S.; Trese, M.T. Cosegregation of two unlinked mutant alleles in some cases of autosomal dominant familial exudative vitreoretinopathy. Eur. J. Hum. Genet. 2003, 12, 79–82. [Google Scholar] [CrossRef]
  6. Li, Y.; Peng, J.; Li, J.; Zhang, Q.; Li, J.; Zhang, X.; Fei, P.; She, K.; Zhao, P. The characteristics of digenic familial exudative vitreoretinopathy. Graefes Arch. Clin. Exp. Ophthalmol. 2018, 256, 2149–2156. [Google Scholar] [CrossRef]
  7. Wang, Y.; Cho, C.; Williams, J.; Smallwood, P.M.; Zhang, C.; Junge, H.J.; Nathans, J. Interplay of the Norrin and Wnt7a/Wnt7b signaling systems in blood–brain barrier and blood–retina barrier development and maintenance. Proc. Natl. Acad. Sci. USA 2018, 115, E11827–E11836. [Google Scholar] [CrossRef] [Green Version]
  8. Collin, R.W.J.; Nikopoulos, K.; Dona, M.; Gilissen, C.; Hoischen, A.; Boonstra, F.N.; Poulter, J.A.; Kondo, H.; Berger, W.; Toomes, C.; et al. ZNF408 is mutated in familial exudative vitreoretinopathy and is crucial for the development of zebrafish retinal vasculature. Proc. Natl. Acad. Sci. USA 2013, 110, 9856–9861. [Google Scholar] [CrossRef] [Green Version]
  9. Robitaille, J.M.; Gillett, R.M.; LeBlanc, M.A.; Gaston, D.; Nightingale, M.; Mackley, M.P.; Parkash, S.; Hathaway, J.; Thomas, A.; Ells, A.; et al. Phenotypic Overlap Between Familial Exudative Vitreoretinopathy and Microcephaly, Lymphedema, and Chorioretinal Dysplasia Caused byKIF11Mutations. JAMA Ophthalmol. 2014, 132, 1393–1399. [Google Scholar] [CrossRef] [Green Version]
  10. Hu, H.; Xiao, X.; Li, S.; Jia, X.; Guo, X.; Zhang, Q. KIF11mutations are a common cause of autosomal dominant familial exudative vitreoretinopathy. Br. J. Ophthalmol. 2015, 100, 278–283. [Google Scholar] [CrossRef] [Green Version]
  11. Li, V.; Ng, S.S.; Boersema, P.J.; Low, T.Y.; Karthaus, W.R.; Gerlach, J.; Mohammed, S.; Heck, A.; Maurice, M.; Mahmoudi, T.; et al. Wnt Signaling through Inhibition of β-Catenin Degradation in an Intact Axin1 Complex. Cell 2012, 149, 1245–1256. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Tian, T.; Zhang, X.; Zhang, Q.; Zhao, P. Variable reduction in Norrin signaling activity caused by novel mutations in FZD4 identified in patients with familial exudative vitreoretinopathy. Mol. Vis. 2019, 25, 60–69. [Google Scholar] [PubMed]
  13. Pendergast, S.D.; Trese, M.T. Familial exudative vitreoretinopathy. Results of surgical management. Ophthalmology 1998, 105, 1015–1023. [Google Scholar] [CrossRef]
  14. Ranchod, T.M.; Ho, L.Y.; Drenser, K.A.; Capone, A.; Trese, M.T. Clinical Presentation of Familial Exudative Vitreoretinopathy. Ophthalmology 2011, 118, 2070–2075. [Google Scholar] [CrossRef]
  15. Ohlmann, A.; Scholz, M.; Goldwich, A.; Chauhan, B.K.; Hudl, K.; Ohlmann, A.V.; Zrenner, E.; Berger, W.; Cvekl, A.; Seeliger, M.W.; et al. Ectopic Norrin Induces Growth of Ocular Capillaries and Restores Normal Retinal Angiogenesis in Norrie Disease Mutant Mice. J. Neurosci. 2005, 25, 1701–1710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Karjosukarso, D.W.; Van Gestel, S.H.C.; Qu, J.; Kouwenhoven, E.N.; Duijkers, L.; Garanto, A.; Zhou, H.; Collin, R.W.J. An FEVR-associated mutation in ZNF408 alters the expression of genes involved in the development of vasculature. Hum. Mol. Genet. 2018, 27, 3519–3527. [Google Scholar] [CrossRef]
  17. Li, J.-K.; Li, Y.; Zhang, X.; Chen, C.-L.; Rao, Y.-Q.; Fei, P.; Zhang, Q.; Zhao, P.; Li, J. Spectrum of Variants in 389 Chinese Probands With Familial Exudative Vitreoretinopathy. Investig. Opthalmol. Vis. Sci. 2018, 59, 5368–5381. [Google Scholar] [CrossRef] [Green Version]
  18. Xiao, H.; Tong, Y.; Zhu, Y.; Peng, M. Familial Exudative Vitreoretinopathy-Related Disease-Causing Genes and Norrin/ β -Catenin Signal Pathway: Structure, Function, and Mutation Spectrums. J. Ophthalmol. 2019, 2019, 5782536. [Google Scholar] [CrossRef]
  19. Seo, S.H.; Yu, Y.S.; Park, S.W.; Kim, J.H.; Kim, H.K.; Cho, S.I.; Park, H.; Lee, S.J.; Seong, M.-W.; Park, S.S.; et al. Molecular Characterization of FZD4, LRP5, and TSPAN12 in Familial Exudative Vitreoretinopathy. Investig. Opthalmol. Vis. Sci. 2015, 56, 5143. [Google Scholar] [CrossRef] [Green Version]
  20. Taylor, R.L.; Parry, N.R.; Barton, S.J.; Campbell, C.; Delaney, C.M.; Ellingford, J.; Hall, G.; Hardcastle, C.; Morarji, J.; Nichol, E.J.; et al. Panel-Based Clinical Genetic Testing in 85 Children with Inherited Retinal Disease. Ophthalmology 2017, 124, 985–991. [Google Scholar] [CrossRef]
  21. LaDuca, H.; Farwell, K.D.; Vuong, H.; Lu, H.-M.; Mu, W.; Shahmirzadi, L.; Tang, S.; Chen, J.; Bhide, S.; Chao, E.C. Exome sequencing covers >98% of mutations identified on targeted next generation sequencing panels. PLoS ONE 2017, 12, e0170843. [Google Scholar] [CrossRef] [PubMed]
  22. Consugar, M.B.; Navarro-Gomez, D.; Place, E.M.; Bujakowska, K.M.; Sousa, M.E.; Fonseca-Kelly, Z.D.; Taub, D.G.; Janessian, M.; Wang, D.Y.; Au, E.D.; et al. Panel-based genetic diagnostic testing for inherited eye diseases is highly accurate and reproducible, and more sensitive for variant detection, than exome sequencing. Genet. Med. 2014, 17, 253–261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Miller, D.T.; Lee, K.; Chung, W.K.; Gordon, A.S.; Herman, G.E.; Klein, T.E.; Stewart, D.R.; Amendola, L.M.; Adelman, K.; Bale, S.J.; et al. ACMG SF v3.0 list for reporting of secondary findings in clinical exome and genome sequencing: A policy statement of the American College of Medical Genetics and Genomics (ACMG). Genet. Med. 2021, 23, 1381–1390. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Norrin-based Wnt-signaling in retinal endothelial cells. Lack of Norrin-based Wnt-signaling results in underdevelopment of the microvascular beds and failure to develop a patent high-barrier vascular endothelium in the neural retina and in the cochlea. A Norrin dimer can bind to its central cognate receptor Frizzled-4 (FZD-4) as well as the co-receptors LRP-5 and TSPAN-12. Activation of the receptor complex is increased by the synergy of multiple protein–protein interactions (dashed arrows), inhibiting the ubiquitination of phosphorylated β-Catenin, which then accumulates on the destruction complex [11]. When Norrin is available, this effect causes an increase in cytoplasmic β-Catenin concentration, which enters the nucleus to modulate target gene expression through interactions with the TCF/LEF family of transcription factors as well as other transcription factors that remain unknown.
Figure 1. Norrin-based Wnt-signaling in retinal endothelial cells. Lack of Norrin-based Wnt-signaling results in underdevelopment of the microvascular beds and failure to develop a patent high-barrier vascular endothelium in the neural retina and in the cochlea. A Norrin dimer can bind to its central cognate receptor Frizzled-4 (FZD-4) as well as the co-receptors LRP-5 and TSPAN-12. Activation of the receptor complex is increased by the synergy of multiple protein–protein interactions (dashed arrows), inhibiting the ubiquitination of phosphorylated β-Catenin, which then accumulates on the destruction complex [11]. When Norrin is available, this effect causes an increase in cytoplasmic β-Catenin concentration, which enters the nucleus to modulate target gene expression through interactions with the TCF/LEF family of transcription factors as well as other transcription factors that remain unknown.
Genes 13 00495 g001
Figure 2. Final group distribution of 76 sequenced subjects. The largest grouping contained patients confirmed in the ARC clinic with FEVR (35). ARC confirmed FEVR-Unaffected subjects (U) numbered 20. Five pediatric subjects were categorized as ROP with potential FEVR complication (FROP), while one pediatric subject was ROP. Fifteen subjects, mostly relatives of FEVR patients, were not available to ARC for retinal diagnosis (NA).
Figure 2. Final group distribution of 76 sequenced subjects. The largest grouping contained patients confirmed in the ARC clinic with FEVR (35). ARC confirmed FEVR-Unaffected subjects (U) numbered 20. Five pediatric subjects were categorized as ROP with potential FEVR complication (FROP), while one pediatric subject was ROP. Fifteen subjects, mostly relatives of FEVR patients, were not available to ARC for retinal diagnosis (NA).
Genes 13 00495 g002
Figure 3. Gene distribution of protein sequence variants detected in the FEVR cohort. Of the 33 protein-altering variants listed in Table 4, the distribution by gene was as follows: LRP5 13/33 (40%), FZD4 9/33 (27%), ZNF408 6/33 (18%), (KIF11 3/33 (9%), NDP 1/33 (3%), CTNNB1 1/33 (3%).
Figure 3. Gene distribution of protein sequence variants detected in the FEVR cohort. Of the 33 protein-altering variants listed in Table 4, the distribution by gene was as follows: LRP5 13/33 (40%), FZD4 9/33 (27%), ZNF408 6/33 (18%), (KIF11 3/33 (9%), NDP 1/33 (3%), CTNNB1 1/33 (3%).
Genes 13 00495 g003
Figure 4. Distribution of the number of affected genes with protein-altering variants in affected and unaffected FEVR subjects. Numbers of genes are shown for 35 FEVR subjects and 20 Unaffected (U) subjects as confirmed by ARC diagnosis.
Figure 4. Distribution of the number of affected genes with protein-altering variants in affected and unaffected FEVR subjects. Numbers of genes are shown for 35 FEVR subjects and 20 Unaffected (U) subjects as confirmed by ARC diagnosis.
Genes 13 00495 g004
Table 1. Sequencing Parameters. General sequencing metrics as obtained with 48 pooled samples on a single iSeq-100 run.
Table 1. Sequencing Parameters. General sequencing metrics as obtained with 48 pooled samples on a single iSeq-100 run.
ParameterValue
Number of amplicon regions180.0
Total length of target regions32,731.0
Percent on-target aligned reads98.3
Percent on-target PF reads92.9
Percent aligned read94.6
Percent on-target aligned bases98.4
Percent on-target PF bases92.2
Percent Q30 bases95.5
Percent aligned bases93.7
Amplicon mean coverage978.0
Table 2. Variant types. Variants detected, by type, from the sequencing of 48 pooled sequencing library samples with the 8-gene panel. Metrics for SNVs (single nucleotide variants), insertions, and deletions are shown as total counts and the average count per sample.
Table 2. Variant types. Variants detected, by type, from the sequencing of 48 pooled sequencing library samples with the 8-gene panel. Metrics for SNVs (single nucleotide variants), insertions, and deletions are shown as total counts and the average count per sample.
Variant TypeTotal CountAverage per Sample
SNVs total93719.5
SNVs in genes92719.3
SNVs in exons2685.6
SNVs in coding regions2204.6
SNVs in UTR regions481.0
SNVs in splice site regions551.1
Stop-gained SNVs20.04
Stop-lost SNVs00.0
Non-synonymous SNVs440.9
Synonymous SNVs1743.6
SNVs not in dbSNP30.1
Insertions total811.7
Insertions in genes280.6
Insertions in exons100.2
Insertions in coding regions80.2
Insertions in UTR regions20.04
Insertions in splice site regions00.0
Stop-gained insertions00.0
Stop-lost insertions00.0
Frame shift insertions20.04
Non-synonymous insertions60.1
Insertions not in dbSNP30.1
Deletions total711.5
Deletions in genes260.5
Deletions in exons190.4
Deletions in coding regions150.3
Deletions in UTR regions40.1
Deletions in splice site regions00.0
Stop-gained deletions00.0
Stop-lost deletions00.0
Frame shift deletions00.0
Non-synonymous deletions150.3
Deletions not in dbSNP100.2
Table 3. Protein-altering variants from targeted sequencing of 76 subjects. Patients were referred to the clinic with FEVR diagnosis or to evaluate potential FEVR status. Some close relatives of FEVR subjects were also sequenced. FEVR grade at clinical presentation is shown for confirmed FEVR group subjects. Other group notations are: U (Unaffected, confirmed not-FEVR), ROP (Retinopathy of Prematurity), FROP (FEVR/ROP, potential ROP overlapping with FEVR but not confirmed FEVR), and Not Available (NA). Known FEVR-causing pathogenic variants, likely pathogenic variants, and their corresponding Allele Frequencies (AF) are shown in bold-face type. Allele frequency estimates from the GnomAD data base are provided as well as the estimated frequency of combinations. Subjects with family members referred to clinic are grouped at the top of the table by family designation. A summary listing of the coding variants in HGVS format is provided in Table 4. (Note: there was no Family #6 in this FEVR-based study).
Table 3. Protein-altering variants from targeted sequencing of 76 subjects. Patients were referred to the clinic with FEVR diagnosis or to evaluate potential FEVR status. Some close relatives of FEVR subjects were also sequenced. FEVR grade at clinical presentation is shown for confirmed FEVR group subjects. Other group notations are: U (Unaffected, confirmed not-FEVR), ROP (Retinopathy of Prematurity), FROP (FEVR/ROP, potential ROP overlapping with FEVR but not confirmed FEVR), and Not Available (NA). Known FEVR-causing pathogenic variants, likely pathogenic variants, and their corresponding Allele Frequencies (AF) are shown in bold-face type. Allele frequency estimates from the GnomAD data base are provided as well as the estimated frequency of combinations. Subjects with family members referred to clinic are grouped at the top of the table by family designation. A summary listing of the coding variants in HGVS format is provided in Table 4. (Note: there was no Family #6 in this FEVR-based study).
GroupFEVR GradeFamilyGenderAgeRelationshipCoding-VariationsAF 1AF 2AF 3AF 4Combined AF
1NANA1M43maternal uncle of 4 and 5LRP5:Val 667Met LRP5:Ala1330Val0.0340.134NANA0.004556
2NANA1F44maternal aunt of 4 and 5LRP5:Val 667Met LRP5:Ala1330Val0.0340.134NANA0.004556
3NANA1F41mother of 4 and 5LRP5:Val 667Met LRP5:Ala1330Val0.0340.134NANA0.004556
4FEVR51M1pediatric, son, proband 1LRP5:Val 667Met LRP5:Ala1330Val (homo)0.0340.016NANA0.000544
5FEVR41M1pediatric, twin of 4LRP5:Val 667Met LRP5:Ala1330Val0.0340.134NANA0.004556
6FEVR12F42mother of 7LRP5: Leu20dup ZNF408:Arg337Pro0.10140.0086NANA0.000872
7FEVR22M1pediatric, son, proband 2CTTNB1:Val273Met0.000021NANANA0.000021
8NANA3M67father of 9LRP5:Ala1330Val0.134NANANA0.134
9FEVR43M29son, proband 3LRP5:Leu20dup LRP5:Ala1330Val (homo)0.10140.016NANA0.001622
10U04F28mother of pediatric not collected yetLRP5: Met1086Val ZNF408:Arg337Pro0.0020770.0086NANA0.000018
11U04M30father of pediatric not collected yetNo protein changesNANANANANA
12U05F42mother of 14LRP5:Leu20dup ZNF408:Val194_Val197del0.10140.1307NANA0.013253
13FEVR15M54father of 14LRP5:Val667Met LRP5:Ala1330Val ZNF408:Val194_Val197del0.0340.1340.1307NA0.000595
14FEVR35F3pediatric, daughter, proband-5LRP5:Leu20dup ZNF408:Val194_Val197del0.10140.1307NANA0.013253
15NANA7M45father of 17 and 18FZD4:Met105Val LRP5:Ala1330Val ZNF408:Val194_Val197del0.0000430.1340.1307NA0.000001
16NANA7F43mother of 17No protein changesNANANANANA
17FEVR27F16pediatric, daughter, proband-7FZD4:Met105Val LRP5:Ala1330Val0.0000430.134NANA0.000006
18NANA7F19sibling of 17ZNF408:Val194_Val197del0.1307NANANA0.1307
19FEVR18F14pediatric, daughter, proband-8KIF11:Glu129Ala LRP5:Ala1330Val (homo)00.016NANA0
20U08F11pediatric, sibling of 19ZNF408:Val194_Val197del0.1307NANANA0.1307
21U08M13pediatric, sibling of 19LRP5:Ala1330Val0.134NANANA0.134
22U09M31father of 23No protein changesNANANANANA
23FEVR59M1pediatric, son, proband-9FZD4:Met105Val0.000043NANANA0.000043
24FEVR110M57maternal uncle of 26LRP5:c.4488+2T>G (splice donor)0NANANA0
25FEVR110F45mother of 26LRP5:c.4488+2T>G (splice donor) LRP5:Leu20dup00.1014NANA0
26FEVR510F18daughter, proband-10LRP5:c.4488+2T>G (splice donor) FZD4:Pro11Gln00.000021NANA0
27U010M17son of maternal uncle 24No protein changesNANANANANA
28U011M68maternal grandfather of 30LRP5:Leu20dup ZNF408:Val194_Val197del0.10140.1307NANA0.013253
29FEVR111F30mother of 30FZD4:Cys450ter LRP5:Leu20dup00.1014NANA0
30FEVR511F3pediatric, daughter, proband-11FZD4:Cys450ter FZD4:Met105Val00.000043NANA0
31NANA12M28father of 33LRP5:Val667Met LRP5:Ala1330Val ZNF408:Val194_Val197del0.0340.1340.1307NA0.000595
32NANA12F26mother of 33No protein changesNANANANANA
33FEVR512F1pediatric, daughter, proband-12LRP5:Val667Met LRP5:Ala1330Val0.0340.134 NA0.004556
34NANA13F74grandmother of 38No protein changesNANANANANA
35FEVR113M73grandfather of 38LRP5:Val667Met LRP5:Ala1330Val FZD4:Pro33Ser FZD4:Pro168Ser0.0340.1340.012360.0140770.000001
36U013M49father of 39LRP5:Val667Met LRP5:Ala1330Val ZNF408:Val194_Val197del0.0340.1340.1307NA0.000595
37FEVR113F47mother of 38No protein changesNANANANANA
38FEVR413M11pediatric, son, proband-13LRP5:Val667Met LRP5:Ala1330Val0.0340.134NANA0.004556
39U013F12pediatric, daughter, sibling of 38LRP5:Val667Met LRP5:Ala1330Val0.0340.134NANA0.004556
40FEVR114F46mother of 41LRP5:Leu16_Leu20del ZNF408:Val194_Val197del0.0000070.1307NANA0.000001
41FEVR414M6pediatric, son, proband-14LRP5:c.4488+2T>G (splice donor)0NANANA0
42FEVR215M3pediatric, son, twin of 43ZNF408:Glu230Gly0.000371NANANA0.000371
43FEVR115M3pediatric, son, twin of 42ZNF408:Glu230Gly0.000371NANANA0.000371
44NANA16F8pediatric, daughter, sibling of 47LRP5:Ala1330Val ZNF408:Val194_Val197del0.1340.1307NANA0.017514
45U016M11sibling of 46LRP5: Ala1330Val ZNF408:Val194_Val197del0.1340.1307NANA0.017514
46NANA17MNAfather of child (FEVR not sequenced)ZNF408:Val194_Val197del0.1307NANANA0.1307
47NANA17F37mother of child (FEVR not sequenced)ZNF408:Val194_Val197del0.1307NANANA0.1307
48NANA18M41father of child (FEVR not sequenced)FZD4:Arg127Cys LRP5:Val667Met LRP5:Ala1330Val0.0000460.0340.134NA0
49U018F38mother of child (FEVR not sequenced)No protein changesNANANANA0
50NANA19F60mother of child (FEVR not sequenced)LRP5:Pro6Thr ZNF408:Val194_Val197del0.0029320.1307NANA0.000383
51FEVR3NAF2pediatric, no relationsLRP5:Val667Met LRP5:Ala1330Val LRP5:Pro848Leu LRP5:Thr852Met0.0340.13400.0000130
52FEVR1NAF90adult, no relationsKIF11:His526Gln ZNF408:Val194_Val197del0.0025320.1307NANA0.000331
53FEVR2NAM11pediatric, no relationsLRP5:Ala1330Val ZNF408:Val194_Val197del(homo)0.1340.0159NANA0.002131
54FEVR4NAF6pediatric, no relationsLRP5:Leu20dup ZNF408:Val194_Val197del0.10140.1307NANA0.013253
55FEVR1NAF54adult, no relationsLRP5:Ala1330Val0.134NANANA0.134
56FEVR4NAM9pediatric, no relationsNDP:His42Arg0.000005NANANA0.000005
57FEVR2NAM1pediatric, no relationsLRP5:Val667Met LRP5:Ala1330Val ZNF408:Val194_Val197del0.0340.1340.1307NA0.000595
58FEVR5NAF1pediatric, no relationsLRP5:Cys913LeufsTer730.000004NANANA0.000004
59FEVR5NAF1pediatric, no relationsLRP5:Ala919CysfsTer670NANANA0
60FEVR1NAM14pediatric, no relationsLRP5:Leu20dup0.1014NANANA0.1014
61FEVR5NAM1pediatric, no relationsLRP5:Pro1522Leu FZD4:Ala408del0.0002910NANA0
62FEVR3NAF34adult, no relationsLRP5:Val667Met LRP5:Ala1330Val ZNF408:Val194_Val197del FZD4 Gly161Arg0.0340.1340.130700
63FROP2NAM2pediatric, no relationsNo protein changesNANANANANA
64FROPNANAF1pediatric, no relationsZNF408:Val194_Val197del0.1307NANANA0.1307
65FROP5NAF1pediatric, no relationsLRP5:Gln89Arg0.008262NANANA0.008262
66FROP5NAF1pediatric, no relationsNo protein changesNANANANANA
67FROP4NAF39adult, no relationsLRP5:Ala1330Val0.134NANANA0.134
68U0NAM1pediatric, no relationsLRP5:Leu20dup0.1014NANANA0.1014
69U0NAM7pediatric, no relationsZNF408:Val194_Val197del0.1307NANANA0.1307
70ROPNANAF1pediatric, no relationsZNF408:Val194_Val197del FZD4:Trp139Ser0.13070NANA0
71U0NAF9adult, no relationsLRP5:Leu20dup0.1014NANANA0.1014
72U0NAF15pediatric, no relationsNo protein changesNANANANANA
73U0NAM14pediatric, no relationsKIF11:Pro642Ala0.000293NANANA0.000293
74U0NAF69adult, no relationsLRP5:Leu20dup0.1014NANANA0.1014
75U0NAM43adult, no relationsNo protein changesNANANANANA
76U0NAF6pediatric, no relationsZNF408:Leu67Val ZNF408:Ala372Thr ZNF408:Pro647Gln0.003680.0005770NA0
Table 4. Variants altering the primary amino acid sequence. For six of the seven FEVR genes sequenced (NDP, CTNNB1, KIF11, FZD4, LRP5, ZNF408), thirty-three variants were detected that alter the primary amino acid sequence. Only synonymous variants were found in TSPAN12. Seven new variants are noted, including a likely pathogenic Cys450Ter stop-gained variant in FZD4. Allele Frequency (AF %) was derived from GnomDB. Accession numbers in dbSNP and ClinVar are shown where available. Consequence categories were from ClinVar, except those listed in parentheses, which are the current determinations of the authors.
Table 4. Variants altering the primary amino acid sequence. For six of the seven FEVR genes sequenced (NDP, CTNNB1, KIF11, FZD4, LRP5, ZNF408), thirty-three variants were detected that alter the primary amino acid sequence. Only synonymous variants were found in TSPAN12. Seven new variants are noted, including a likely pathogenic Cys450Ter stop-gained variant in FZD4. Allele Frequency (AF %) was derived from GnomDB. Accession numbers in dbSNP and ClinVar are shown where available. Consequence categories were from ClinVar, except those listed in parentheses, which are the current determinations of the authors.
No.GeneNucleotideProteinMolecular ChangeConsequence AF %AF (Homo) %dbSNPClinVar
Accession
1CTNNB1NM_001098209.2:c.817G>ANP_001091679.1:p.Val273Metmissense variant(Uncertain significance)0.0021%0.0%rs1183899293Not in ClinVar
2FZD4NM_012193.4:c.32C>ANP_036325.2:p.Pro11Glnmissense variant(Uncertain significance)0.0021%0.0%rs766393047Not in ClinVar
3FZD4NM_012193.4:c.97C>TNP_036325.2:p.Pro33Sermissense variantBenign1.236%0.89%rs61735304RCV000387944.3
4FZD4NM_012193.4:c.313A>GNP_036325.2:p.Met105Valmissense variantPathogenic0.0043%0.0%rs80358284RCV000210241.1
5FZD4NM_012193.4:c.379C>T NP_036325.2:p.Arg127Cysmissense variantLikely benign0.0046%0.0%rs376854255RCV001111581.1
6 NewFZD4NM_012193.4:c.416G>CNP_036325.2:p.Trp139Sermissense variant(Uncertain significance)NANANANot in ClinVar
7 NewFZD4NM_012193.4:c.481G>CNP_036325.2:p.Gly161Argmissense variant(Uncertain significance)NANANANot in ClinVar
8FZD4NM_012193.4:c.502C>TNP_036325.2:p.Pro168Sermissense variantBenign1.4077%0.01379%rs61735303RCV000368489.2
9 NewFZD4NM_012193.4:c.1221_1223delTCGNP_036325.2:p.Ala408delin-frame deletion(Uncertain significance)NANANANot in ClinVar
10 NewFZD4NM_012193.4:c.1350T>ANP_036325.2:p.Cys450Terstop-gained(Likely pathogenic)NANANANot in ClinVar
11KIF11NM_004523.4:c.386A>CNP_004514.2:p.Glu129Alamissense variant(Uncertain significance)0.00004 NArs779558239Not in ClinVar
12KIF11NM_004523.4:c.1924C>GNP_004514.2:p.Pro642Alamissense variantBenign0.02930%0.0%rs79865214RCV000915477.2
13KIF11NM_004523.4:c.1578C>ANP_004514.2:p.His526Glnmissense variantBenign0.25320%0.0013%rs112145870RCV000967853.2
14LRP5NM_002335.4:c.16C>ANP_002326.2:p.Pro6Thrmissense variantBenign0.2932%0.0035%rs771718186RCV000592927.1
15LRP5NM_002335.4:c.34_36CTG [4]NP_002326.2:p.Leu16_Leu20deldeletionUncertain significance0.0006928%0.0%rs72555376Not in ClinVar
16LRP5NM_002335.4:c.58_60dupCTGNP_002326.2:p.Leu20dupinsert duplicationBenign10.14%0.57%rs564221347VCV000193231.2
17LRP5NM_002335.4:c.266A>GNP_002326.2:p.Gln89Argmissense variantBenign0.83%NArs41494349RCV000175719.2
18LRP5NM_002335.4:c.1999G>ANP_002326.2:p.Val667Metmissense variantLikely benign3.40%0.092%rs4988321RCV000250939.1
19 NewLRP5NM_002335.4:c.2543C>TNP_002326.2:p.Pro848Leumissense variantNANANANANot in ClinVar
20LRP5NM_002335.4:c.2555C>TNP_002326.2:p.Thr852Metmissense variantNA0.0013140%0.0%rs1398692057Not in ClinVar
21LRP5NM_002335.4:c.2737dupNP_002326.2:p.Cys913fsTer73frameshift variantPathogenic0.0004000%0.0%rs886043590RCV000761295.1
22 NewLRP5NM_002335.4:c.2754dupNP_002326.2:p.Ala919CysfsTer67frameshift variant(Likely pathogenic)NANANANot in ClinVar
23LRP5NM_002335.4:c.3256A>GNP_002326.2:p.Met1086Valmissense variantLikely benign0.20770%0.0019730%rs145774832RCV000592263.1
24LRP5NM_002335.4:c.3989C>TNP_002326.2:p.Ala1330Valmissense variantBenign13.40%1.60%rs3736228RCV000242123.2
25LPR5NM_002335.4:c.4565C>TNP_002326.2:p.Pro1522Leumissense variantUncertain significance0.02910%0.0%rs200624778RCV000724481.4
26LRP5NM_002335.4:c.4488+2T>Gintron 21 donor site lostsplice donorPathogenicNANArs80358322RCV000006667.3
27NDPNM_000266.4:c.125A>GNP_000257.1:p.His42Argmissense variantPathogenic0.0005%NArs104894874RCV000011437.5
28ZNF408NM_024741.3:c.199C>GNP_079017.1:p.Leu67Valmissense variantBenign0.3680%0.0019710%rs35652367RCV000969067.2
29ZNF408NM_024741.3:c.581_592delNP_079017.1:p.Val194_Val197deldeletionBenign13.1%1.59%rs148055528RCV000248955.1
30ZNF408NM_024741.3:c.689A>GNP_079017.1:p.Glu230Glymissense variantLikely benign0.0371%0.0006570%rs147850078RCV000974525.2
31ZNF408NM_024741.3:c.1010G>CNP_079017.1:p.Arg337Promissense variantBenign0.8609%0.0151100%rs36017347RCV000959911.2
32ZNF408NM_024741.3:c.1114G>ANP_079017.1:p.Ala372Thrmissense variantUncertain significance0.0577%0.0%rs141624151RCV001053827.1
33 NewZNF408NM_024741.3:c.1940_1941delinsAANP_079017.1:p.Pro647Glnmissense variant(Uncertain significance)NANANANot in ClinVar
Table 5. Single and multigenic variants combinations. Protein-altering variants were found in multiple combinations in many patients. Multiple variants in single genes are indicated within parentheses. Totals are for the numbers of unique, different combinations detected. Genes are listed for the digenic and trigenic variant combinations.
Table 5. Single and multigenic variants combinations. Protein-altering variants were found in multiple combinations in many patients. Multiple variants in single genes are indicated within parentheses. Totals are for the numbers of unique, different combinations detected. Genes are listed for the digenic and trigenic variant combinations.
Variant NatureGene/sTotal
monogenicNDP1
monogenicLRP516
monogenicLRP5 (2)7
monogenicLRP5 (3)1
monogenicFZD42
monogenicFZD4 (2)1
monogenicCTNNB11
monogenicZNF4081
monogenicZNF408 (3)1
monogenicKIFll1
digenicLRP5, ZNF40814
digenicLRP5, FZD46
digenicLRP5, KIF112
digenicKIF11, ZNF4081
digenicFZD4, ZNF4081
trigenicLRP5, ZNF408, FZD42
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cicerone, A.P.; Dailey, W.; Sun, M.; Santos, A.; Jeong, D.; Jones, L.; Koustas, K.; Drekh, M.; Schmitz, K.; Haque, N.; et al. A Survey of Multigenic Protein-Altering Variant Frequency in Familial Exudative Vitreo-Retinopathy (FEVR) Patients by Targeted Sequencing of Seven FEVR-Linked Genes. Genes 2022, 13, 495. https://doi.org/10.3390/genes13030495

AMA Style

Cicerone AP, Dailey W, Sun M, Santos A, Jeong D, Jones L, Koustas K, Drekh M, Schmitz K, Haque N, et al. A Survey of Multigenic Protein-Altering Variant Frequency in Familial Exudative Vitreo-Retinopathy (FEVR) Patients by Targeted Sequencing of Seven FEVR-Linked Genes. Genes. 2022; 13(3):495. https://doi.org/10.3390/genes13030495

Chicago/Turabian Style

Cicerone, Amanda Petrelli, Wendy Dailey, Michael Sun, Andrew Santos, Daeun Jeong, Lance Jones, Konstaninos Koustas, Mary Drekh, Keaton Schmitz, Naomi Haque, and et al. 2022. "A Survey of Multigenic Protein-Altering Variant Frequency in Familial Exudative Vitreo-Retinopathy (FEVR) Patients by Targeted Sequencing of Seven FEVR-Linked Genes" Genes 13, no. 3: 495. https://doi.org/10.3390/genes13030495

APA Style

Cicerone, A. P., Dailey, W., Sun, M., Santos, A., Jeong, D., Jones, L., Koustas, K., Drekh, M., Schmitz, K., Haque, N., Felisky, J. A., Guzman, A. E., Mellert, K., Trese, M. T., Capone, A., Drenser, K. A., & Mitton, K. P. (2022). A Survey of Multigenic Protein-Altering Variant Frequency in Familial Exudative Vitreo-Retinopathy (FEVR) Patients by Targeted Sequencing of Seven FEVR-Linked Genes. Genes, 13(3), 495. https://doi.org/10.3390/genes13030495

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop