Snijders Blok-Campeau Syndrome: Description of 20 Additional Individuals with Variants in CHD3 and Literature Review.

Snijders Blok-Campeau syndrome (SNIBCPS, OMIM# 618205) is an extremely infrequent disease with only approximately 60 cases reported so far. SNIBCPS belongs to the group of neurodevelopmental disorders (NDDs). Clinical features of patients with SNIBCPS include global developmental delay, intellectual disability, speech and language difficulties and behavioral disorders like autism spectrum disorder. In addition, patients with SNIBCPS exhibit typical dysmorphic features including macrocephaly, hypertelorism, sparse eyebrows, broad forehead, prominent nose and pointed chin. The severity of the neurological effects as well as the presence of other features is variable among subjects. SNIBCPS is caused likely by pathogenic and pathogenic variants in CHD3 (Chromodomain Helicase DNA Binding Protein 3), which seems to be involved in chromatin remodeling by deacetylating histones. Here, we report 20 additional patients with clinical features compatible with SNIBCPS from 17 unrelated families with confirmed likely pathogenic/pathogenic variants in CHD3. Patients were analyzed by whole exome sequencing and segregation studies were performed by Sanger sequencing. Patients in this study showed different pathogenic variants affecting several functional domains of the protein. Additionally, none of the variants described here were reported in control population databases, and most computational predictors suggest that they are deleterious. The most common clinical features of the whole cohort of patients are global developmental delay (98%) and speech disorder/delay (92%). Other frequent features (51-74%) include intellectual disability, hypotonia, hypertelorism, abnormality of vision, macrocephaly and prominent forehead, among others. This study expands the number of individuals with confirmed SNIBCPS due to pathogenic or likely pathogenic variants in CHD3. Furthermore, we add evidence of the importance of the application of massive parallel sequencing for NDD patients for whom the clinical diagnosis might be challenging and where deep phenotyping is extremely useful to accurately manage and follow up the patients.

The neurodevelopmental and facial phenotype in individuals with a TRIP12 variant.

Haploinsufficiency of TRIP12 causes a neurodevelopmental disorder characterized by intellectual disability associated with epilepsy, autism spectrum disorder and dysmorphic features, also named Clark-Baraitser syndrome. Only a limited number of cases have been reported to date. We aimed to further delineate the TRIP12-associated phenotype and objectify characteristic facial traits through GestaltMatcher image analysis based on deep-learning algorithms in order to establish a TRIP12 gestalt. 38 individuals between 3 and 66 years (F = 20, M = 18) - 1 previously published and 37 novel individuals - were recruited through an ERN ITHACA call for collaboration. 35 TRIP12 variants were identified, including frameshift (n = 15) and nonsense (n = 6) variants, as well as missense (n = 5) and splice (n = 3) variants, intragenic deletions (n = 4) and two multigene deletions disrupting TRIP12. Though variable in severity, global developmental delay was noted in all individuals, with language deficit most pronounced. About half showed autistic features and susceptibility to obesity seemed inherent to this disorder. A more severe expression was noted in individuals with a missense variant. Facial analysis showed a clear gestalt including deep-set eyes with narrow palpebral fissures and fullness of the upper eyelids, downturned corners of the mouth and large, often low-set ears with prominent earlobes. We report the largest cohort to date of individuals with TRIP12 variants, further delineating the associated phenotype and introducing a facial gestalt. These findings will improve future counseling and patient guidance.

Biallelic intragenic tandem duplication of CPLANE1 in Joubert syndrome: A case report.

Joubert syndrome (JS) is a clinically and genetically heterogeneous genetic disorder. To date, 40 JS-causing genes have been reported and CPLANE1 is one of the most frequently mutated, with biallelic pathogenic missense and truncating variants explaining up to 14% of JS cases. We present a case of JS diagnosed after the identification of a novel biallelic intragenic duplication of exons 20-46 of CPLANE1. The quadruplication was identified by short-read sequencing and copy number variant analysis and confirmed in tandem by long PCR with the breakpoints defined by a nanopore-based long-read sequencing approach. Based on the genetic findings and the clinical presentation of the patient, a brain MRI was ordered, evidencing the molar tooth sign, which confirmed the diagnosis of JS in the patient. This is, to the best of our knowledge, the first report of an intragenic duplication in this gene as the potential molecular mechanism of JS.

Androgen receptor polyQ alleles and COVID-19 severity in men: A replication study.

BACKGROUND: Ample evidence indicates a sex-related difference in severity of COVID-19, with less favorable outcomes observed in men. Genetic factors have been proposed as candidates to explain this difference. The polyglutamine (polyQ) polymorphism in the androgen receptor gene has been recently described as a genetic biomarker of COVID-19 severity. OBJECTIVE: To test the association between the androgen receptor polyQ polymorphism and COVID-19 severity in a large cohort of COVID-19 male patients. MATERIALS AND METHODS: This study included 1136 male patients infected with SARS-CoV-2 as confirmed by positive PCR. Patients were retrospectively and prospectively enrolled from March to November 2020. Patients were classified according to their severity into three categories: oligosymptomatic, hospitalized and severe patients requiring ventilatory support. The number of CAG repeats (polyQ polymorphism) at the androgen receptor was obtained by PCR and patients were classified as either short (<23 repeats) or long (≥23 repeats) allele carriers. The association between polyQ alleles (short or long) and COVID-19 severity was assessed by Chi-squared (Chi2 ) and logistic regression analysis. RESULTS: The mean number of polyQ CAG repeats was 22 (±3). Patients were classified as oligosymptomatic (15.5%), hospitalized (63.2%), and severe patients (21.3%) requiring substantial respiratory support. PolyQ alleles distribution did not show significant differences between severity classes in our cohort (Chi2 test p > 0.05). Similar results were observed after adjusting by known risk factors such as age, comorbidities, and ethnicity (multivariate logistic regression analysis). DISCUSSION: Androgen sensitivity may be a critical factor in COVID-19 disease severity. However, we did not find an association between the polyQ polymorphism and the COVID-19 severity. Additional studies are needed to clarify the mechanism underlying the association between androgens and COVID-19 outcome. CONCLUSIONS: The results obtained in our study do not support the role of this polymorphism as biomarker of COVID-19 severity.

Age-dependent association of clonal hematopoiesis with COVID-19 mortality in patients over 60 years.

Clonal hematopoiesis, especially that of indeterminate potential (CHIP), has been associated with age-related diseases, such as those contributing to a more severe COVID-19. Four studies have attempted to associate CHIP with COVID-19 severity without conclusive findings. In the present work, we explore the association between CHIP and COVID-19 mortality. Genomic DNA extracted from peripheral blood of COVID-19 patients (n = 241 deceased, n = 239 survivors) was sequenced with the Myeloid Solutions™ panel of SOPHiA Genetics. The association between clonality and age and clonality and mortality was studied using logistic regression models adjusted for sex, ethnicity, and comorbidities. The association with mortality was performed with patients stratified into four groups of age according to the quartiles of the distribution: 60-74 years, 75-84 years, 85-91 years, and 92-101 years. Clonality was found in 38% of the cohort. The presence of CHIP variants, but not the number, significantly increased with age in the entire cohort of COVID-19 patients, as well as in the group of survivors (p < 0.001). When patients were stratified by age and the analysis adjusted, CHIP classified as pathogenic/likely pathogenic was significantly more represented in deceased patients compared with survivors in the group of 75-84 years (34.6% vs 13.7%, p = 0.020). We confirmed the well-established linear relationship between age and clonality in the cohort of COVID-19 patients and found a significant association between pathogenic/likely pathogenic CHIP and mortality in patients from 75 to 84 years that needs to be further validated.

Aggregated Genomic Data as Cohort-Specific Allelic Frequencies can Boost Variants and Genes Prioritization in Non-Solved Cases of Inherited Retinal Dystrophies.

The introduction of NGS in genetic diagnosis has increased the repertoire of variants and genes involved and the amount of genomic information produced. We built an allelic-frequency (AF) database for a heterogeneous cohort of genetic diseases to explore the aggregated genomic information and boost diagnosis in inherited retinal dystrophies (IRD). We retrospectively selected 5683 index-cases with clinical exome sequencing tests available, 1766 with IRD and the rest with diverse genetic diseases. We calculated a subcohort's IRD-specific AF and compared it with suitable pseudocontrols. For non-solved IRD cases, we prioritized variants with a significant increment of frequencies, with eight variants that may help to explain the phenotype, and 10/11 of uncertain significance that were reclassified as probably pathogenic according to ACMG. Moreover, we developed a method to highlight genes with more frequent pathogenic variants in IRD cases than in pseudocontrols weighted by the increment of benign variants in the same comparison. We identified 18 genes for further studies that provided new insights in five cases. This resource can also help one to calculate the carrier frequency in IRD genes. A cohort-specific AF database assists with variants and genes prioritization and operates as an engine that provides a new hypothesis in non-solved cases, augmenting the diagnosis rate.

The recurrent TCF4 missense variant p.(Arg389Cys) causes a neurodevelopmental disorder overlapping with but not typical for Pitt-Hopkins syndrome.

TCF4 haploinsufficiency by deletions, truncating variants or loss-of-function missense variants within the DNA-binding and protein interacting bHLH domain causes Pitt-Hopkins syndrome (PTHS). This neurodevelopmental disorder (NDD) is characterized by severe intellectual disability (ID), epilepsy, hyperbreathing and a typical facial gestalt. Only few aberrations of the N-terminus of TCF4 were associated with milder or atypical phenotypes. By personal communication and searching databases we assembled six cases with the novel, recurrent, de novo missense variant c.1165C > T, p.(Arg389Cys) in TCF4. This variant was identified by diagnostic exome or panel sequencing and is located upstream of the bHLH domain. All six individuals presented with moderate to severe ID with language impairment. Microcephaly occurred in two individuals, epilepsy only in one, and no breathing anomalies or myopia were reported. Facial gestalt showed some aspects of PTHS but was rather non-specific in most individuals. Interestingly, the variant is located within the AD2 activation domain next to a highly conserved coactivator-recruitment motif and might alter interaction with coactivator proteins independently from the bHLH domain. Our findings of a recurrent missense variant outside the bHLH domain in six individuals with an ID phenotype overlapping with but not typical for PTHS delineate a novel genotype-phenotype correlation for TCF4-related NDDs.

Presence of rare potential pathogenic variants in subjects under 65 years old with very severe or fatal COVID-19.

Rare variants affecting host defense against pathogens could be involved in COVID-19 severity and may help explain fatal outcomes in young and middle-aged patients. Our aim was to report the presence of rare genetic variants in certain genes, by using whole exome sequencing, in a selected group of COVID-19 patients under 65 years who required intubation or resulting in death (n = 44). To this end, different etiopathogenic mechanisms were explored using gene prioritization-based analysis in which genes involved in immune response, immunodeficiencies or blood coagulation were studied. We detected 44 different variants of interest, in 29 different patients (66%). Some of these variants were previously described as pathogenic and were located in genes mainly involved in immune response. A network analysis, including the 42 genes with candidate variants, showed three main components, consisting of 25 highly interconnected genes related to immune response and two additional networks composed by genes enriched in carbohydrate metabolism and in DNA metabolism and repair processes. In conclusion, we have detected candidate variants that may potentially influence COVID-19 outcome in our cohort of patients. Further studies are needed to confirm the ultimate role of the genetic variants described in the present study on COVID-19 severity.

Novel genes and sex differences in COVID-19 severity.

Here, we describe the results of a genome-wide study conducted in 11 939 coronavirus disease 2019 (COVID-19) positive cases with an extensive clinical information that were recruited from 34 hospitals across Spain (SCOURGE consortium). In sex-disaggregated genome-wide association studies for COVID-19 hospitalization, genome-wide significance (P < 5 × 10-8) was crossed for variants in 3p21.31 and 21q22.11 loci only among males (P = 1.3 × 10-22 and P = 8.1 × 10-12, respectively), and for variants in 9q21.32 near TLE1 only among females (P = 4.4 × 10-8). In a second phase, results were combined with an independent Spanish cohort (1598 COVID-19 cases and 1068 population controls), revealing in the overall analysis two novel risk loci in 9p13.3 and 19q13.12, with fine-mapping prioritized variants functionally associated with AQP3 (P = 2.7 × 10-8) and ARHGAP33 (P = 1.3 × 10-8), respectively. The meta-analysis of both phases with four European studies stratified by sex from the Host Genetics Initiative (HGI) confirmed the association of the 3p21.31 and 21q22.11 loci predominantly in males and replicated a recently reported variant in 11p13 (ELF5, P = 4.1 × 10-8). Six of the COVID-19 HGI discovered loci were replicated and an HGI-based genetic risk score predicted the severity strata in SCOURGE. We also found more SNP-heritability and larger heritability differences by age (<60 or ≥60 years) among males than among females. Parallel genome-wide screening of inbreeding depression in SCOURGE also showed an effect of homozygosity in COVID-19 hospitalization and severity and this effect was stronger among older males. In summary, new candidate genes for COVID-19 severity and evidence supporting genetic disparities among sexes are provided.

An evaluation of pipelines for DNA variant detection can guide a reanalysis protocol to increase the diagnostic ratio of genetic diseases.

Clinical exome (CE) sequencing has become a first-tier diagnostic test for hereditary diseases; however, its diagnostic rate is around 30-50%. In this study, we aimed to increase the diagnostic yield of CE using a custom reanalysis algorithm. Sequencing data were available for three cohorts using two commercial protocols applied as part of the diagnostic process. Using these cohorts, we compared the performance of general and clinically relevant variant calling and the efficacy of an in-house bioinformatic protocol (FJD-pipeline) in detecting causal variants as compared to commercial protocols. On the whole, the FJD-pipeline detected 99.74% of the causal variants identified by the commercial protocol in previously solved cases. In the unsolved cases, FJD-pipeline detects more INDELs and non-exonic variants, and is able to increase the diagnostic yield in 2.5% and 3.2% in the re-analysis of 78 cancer and 62 cardiovascular cases. These results were considered to design a reanalysis, filtering and prioritization algorithm that was tested by reassessing 68 inconclusive cases of monoallelic autosomal recessive retinal dystrophies increasing the diagnosis by 4.4%. In conclusion, a guided NGS reanalysis of unsolved cases increases the diagnostic yield in genetic disorders, making it a useful diagnostic tool in medical genetics.

High SARS-CoV-2 viral load is associated with a worse clinical outcome of COVID-19 disease.

COVID-19 severity and progression are determined by several host and virological factors that may influence the final outcome of SARS-CoV-2-infected patients. The objective of this work was to determine a possible association between viral load, obtained from nasopharyngeal swabs, and the severity of the infection in a cohort of 448 SARS-CoV-2-infected patients from a hospital in Madrid during the first outbreak of the pandemic in Spain. To perform this, we clinically classified patients as mild, moderate and severe COVID-19 according to a number of clinical parameters such as hospitalization requirement, need of oxygen therapy, admission to intensive care units and/or death. Also, Ct values were determined using SARS-CoV-2-specific oligonucleotides directed to ORF1ab. Here we report a statistically significant association between viral load and disease severity, a high viral load being associated with worse clinical prognosis, independently of several previously identified risk factors such as age, sex, hypertension, cardiovascular disease, diabetes, obesity and lung disease (asthma and chronic obstructive pulmonary disease). The data presented here reinforce viral load as a potential biomarker for predicting disease severity in SARS-CoV-2-infected patients. It is also an important parameter in viral evolution since it relates to the numbers and types of variant genomes present in a viral population, a potential determinant of disease progression.

Attention Deficit Hyperactivity and Autism Spectrum Disorders as the Core Symptoms of AUTS2 Syndrome: Description of Five New Patients and Update of the Frequency of Manifestations and Genotype-Phenotype Correlation.

Haploinsufficiency of AUTS2 has been associated with a syndromic form of neurodevelopmental delay characterized by intellectual disability, autistic features, and microcephaly, also known as AUTS2 syndrome. While the phenotype associated with large deletions and duplications of AUTS2 is well established, clinical features of patients harboring AUTS2 sequence variants have not been extensively described. In this study, we describe the phenotype of five new patients with AUTS2 pathogenic variants, three of them harboring loss-of-function sequence variants. The phenotype of the patients was characterized by attention deficit/hyperactivity disorder (ADHD) and autism spectrum disorder (ASD) or autistic features and mild global developmental delay (GDD) or intellectual disability (ID), all in 4/5 patients (80%), a frequency higher than previously reported for ADHD and autistic features. Microcephaly and short stature were found in 60% of the patients; and feeding difficulties, generalized hypotonia, and ptosis, were each found in 40%. We also provide the aggregated frequency of the 32 items included in the AUTS2 syndrome severity score (ASSS) in patients currently reported in the literature. The main characteristics of the syndrome are GDD/ID in 98% of patients, microcephaly in 65%, feeding difficulties in 62%, ADHD or hyperactivity in 54%, and autistic traits in 52%. Finally, using the location of 31 variants from the literature together with variants from the five patients, we found significantly higher ASSS values in patients with pathogenic variants affecting the 3' end of the gene, confirming the genotype-phenotype correlation initially described.

Schuurs-Hoeijmakers Syndrome (PACS1 Neurodevelopmental Disorder): Seven Novel Patients and a Review.

Schuurs-Hoeijmakers syndrome (SHMS) or PACS1 Neurodevelopmental disorder is a rare disorder characterized by intellectual disability, abnormal craniofacial features and congenital malformations. SHMS is an autosomal dominant hereditary disease caused by pathogenic variants in the PACS1 gene. PACS1 is a trans-Golgi-membrane traffic regulator that directs protein cargo and several viral envelope proteins. It is upregulated during human embryonic brain development and has low expression after birth. So far, only 54 patients with SHMS have been reported. In this work, we report on seven new identified SHMS individuals with the classical c.607C > T: p.Arg206Trp PACS1 pathogenic variant and review clinical and molecular aspects of all the patients reported in the literature, providing a summary of clinical findings grouped as very frequent (≥75% of patients), frequent (50-74%), infrequent (26-49%) and rare (less than ≤25%).

Comparison of the diagnostic yield of aCGH and genome-wide sequencing across different neurodevelopmental disorders.

Most consensus recommendations for the genetic diagnosis of neurodevelopmental disorders (NDDs) do not include the use of next generation sequencing (NGS) and are still based on chromosomal microarrays, such as comparative genomic hybridization array (aCGH). This study compares the diagnostic yield obtained by aCGH and clinical exome sequencing in NDD globally and its spectrum of disorders. To that end, 1412 patients clinically diagnosed with NDDs and studied with aCGH were classified into phenotype categories: global developmental delay/intellectual disability (GDD/ID); autism spectrum disorder (ASD); and other NDDs. These categories were further subclassified based on the most frequent accompanying signs and symptoms into isolated forms, forms with epilepsy; forms with micro/macrocephaly and syndromic forms. Two hundred and forty-five patients of the 1412 were subjected to clinical exome sequencing. Diagnostic yield of aCGH and clinical exome sequencing, expressed as the number of solved cases, was compared for each phenotype category and subcategory. Clinical exome sequencing was superior than aCGH for all cases except for isolated ASD, with no additional cases solved by NGS. Globally, clinical exome sequencing solved 20% of cases (versus 5.7% by aCGH) and the diagnostic yield was highest for all forms of GDD/ID and lowest for Other NDDs (7.1% versus 1.4% by aCGH) and ASD (6.1% versus 3% by aCGH). In the majority of cases, diagnostic yield was higher in the phenotype subcategories than in the mother category. These results suggest that NGS could be used as a first-tier test in the diagnostic algorithm of all NDDs followed by aCGH when necessary.

Sanger sequencing is no longer always necessary based on a single-center validation of 1109 NGS variants in 825 clinical exomes.

Despite the improved accuracy of next-generation sequencing (NGS), it is widely accepted that variants need to be validated using Sanger sequencing before reporting. Validation of all NGS variants considerably increases the turnaround time and costs of clinical diagnosis. We comprehensively assessed this need in 1109 variants from 825 clinical exomes, the largest sample set to date assessed using Illumina chemistry reported. With a concordance of 100%, we conclude that Sanger sequencing can be very useful as an internal quality control, but not so much as a verification method for high-quality single-nucleotide and small insertion/deletions variants. Laboratories might validate and establish their own thresholds before discontinuing Sanger confirmation studies. We also expand and validate 23 copy number variations detected by exome sequencing in 20 samples, observing a concordance of 95.65% (22/23).

IL-6-based mortality prediction model for COVID-19: Validation and update in multicenter and second wave cohorts.

BACKGROUND: Coronavirus disease 2019 (COVID-19) is a highly variable condition. Validated tools to assist in the early detection of patients at high risk of mortality can help guide medical decisions. OBJECTIVE: We sought to validate externally, as well as in patients from the second pandemic wave in Europe, our previously developed mortality prediction model for hospitalized COVID-19 patients. METHODS: Three validation cohorts were generated: 2 external with 185 and 730 patients from the first wave and 1 internal with 119 patients from the second wave. The probability of death was calculated for all subjects using our prediction model, which includes peripheral blood oxygen saturation/fraction of inspired oxygen ratio, neutrophil-to-lymphocyte ratio, lactate dehydrogenase, IL-6, and age. Discrimination and calibration were evaluated in the validation cohorts. The prediction model was updated by reestimating individual risk factor effects in the overall cohort (N = 1477). RESULTS: The mortality prediction model showed good performance in the external validation cohorts 1 and 2, and in the second wave validation cohort 3 (area under the receiver-operating characteristic curve, 0.94, 0.86, and 0.86, respectively), with excellent calibration (calibration slope, 0.86, 0.94, and 0.79; intercept, 0.05, 0.03, and 0.10, respectively). The updated model accurately predicted mortality in the overall cohort (area under the receiver-operating characteristic curve, 0.91), which included patients from both the first and second COVID-19 waves. The updated model was also useful to predict fatal outcome in patients without respiratory distress at the time of evaluation. CONCLUSIONS: This is the first COVID-19 mortality prediction model validated in patients from the first and second pandemic waves. The COR+12 online calculator is freely available to facilitate its implementation (https://utrero-rico.shinyapps.io/COR12_Score/).

Prioritizing variants of uncertain significance for reclassification using a rule-based algorithm in inherited retinal dystrophies.

Inherited retinal dystrophies (IRD) are a highly heterogeneous group of rare diseases with a molecular diagnostic rate of >50%. Reclassification of variants of uncertain significance (VUS) poses a challenge for IRD diagnosis. We collected 668 IRD cases analyzed by our geneticists using two different clinical exome-sequencing tests. We identified 114 unsolved cases pending reclassification of 125 VUS and studied their genomic, functional, and laboratory-specific features, comparing them to pathogenic and likely pathogenic variants from the same cohort (N = 390). While the clinical exome used did not show differences in diagnostic rate, the more IRD-experienced geneticist reported more VUS (p = 4.07e-04). Significantly fewer VUS were reported in recessive cases (p = 2.14e-04) compared to other inheritance patterns, and of all the genes analyzed, ABCA4 and IMPG2 had the lowest and highest VUS frequencies, respectively (p = 3.89e-04, p = 6.93e-03). Moreover, few frameshift and stop-gain variants were found to be informed VUS (p = 6.73e-08 and p = 2.93e-06). Last, we applied five pathogenicity predictors and found there is a significant proof of deleteriousness when all score for pathogenicity in missense variants. Altogether, these results provided input for a set of rules that correctly reclassified ~70% of VUS as pathogenic in validation datasets. Disease- and setting-specific features influence VUS reporting. Comparison with pathogenic and likely pathogenic variants can prioritize VUS more likely to be reclassified as causal.

Genetic landscape of 6089 inherited retinal dystrophies affected cases in Spain and their therapeutic and extended epidemiological implications.

Inherited retinal diseases (IRDs), defined by dysfunction or progressive loss of photoreceptors, are disorders characterized by elevated heterogeneity, both at the clinical and genetic levels. Our main goal was to address the genetic landscape of IRD in the largest cohort of Spanish patients reported to date. A retrospective hospital-based cross-sectional study was carried out on 6089 IRD affected individuals (from 4403 unrelated families), referred for genetic testing from all the Spanish autonomous communities. Clinical, demographic and familiar data were collected from each patient, including family pedigree, age of appearance of visual symptoms, presence of any systemic findings and geographical origin. Genetic studies were performed to the 3951 families with available DNA using different molecular techniques. Overall, 53.2% (2100/3951) of the studied families were genetically characterized, and 1549 different likely causative variants in 142 genes were identified. The most common phenotype encountered is retinitis pigmentosa (RP) (55.6% of families, 2447/4403). The most recurrently mutated genes were PRPH2, ABCA4 and RS1 in autosomal dominant (AD), autosomal recessive (AR) and X-linked (XL) NON-RP cases, respectively; RHO, USH2A and RPGR in AD, AR and XL for non-syndromic RP; and USH2A and MYO7A in syndromic IRD. Pathogenic variants c.3386G > T (p.Arg1129Leu) in ABCA4 and c.2276G > T (p.Cys759Phe) in USH2A were the most frequent variants identified. Our study provides the general landscape for IRD in Spain, reporting the largest cohort ever presented. Our results have important implications for genetic diagnosis, counselling and new therapeutic strategies to both the Spanish population and other related populations.