Specialization of the photoreceptor transcriptome by Srrm3-dependent microexons is required for outer segment maintenance and vision.

Retinal photoreceptors have a distinct transcriptomic profile compared to other neuronal subtypes, likely reflecting their unique cellular morphology and function in the detection of light stimuli by way of the ciliary outer segment. We discovered a layer of this molecular specialization by revealing that the vertebrate retina expresses the largest number of tissue-enriched microexons of all tissue types. A subset of these microexons is included exclusively in photoreceptor transcripts, particularly in genes involved in cilia biogenesis and vesicle-mediated transport. This microexon program is regulated by Srrm3, a paralog of the neural microexon regulator Srrm4. Despite the fact that both proteins positively regulate retina microexons in vitro, only Srrm3 is highly expressed in mature photoreceptors. Its deletion in zebrafish results in widespread down-regulation of microexon inclusion from early developmental stages, followed by other transcriptomic alterations, severe photoreceptor defects, and blindness. These results shed light on the transcriptomic specialization and functionality of photoreceptors, uncovering unique cell type-specific roles for Srrm3 and microexons with implications for retinal diseases.

Light-responsive microRNA miR-211 targets Ezrin to modulate lysosomal biogenesis and retinal cell clearance.

Vertebrate vision relies on the daily phagocytosis and lysosomal degradation of photoreceptor outer segments (POS) within the retinal pigment epithelium (RPE). However, how these events are controlled by light is largely unknown. Here, we show that the light-responsive miR-211 controls lysosomal biogenesis at the beginning of light-dark transitions in the RPE by targeting Ezrin, a cytoskeleton-associated protein essential for the regulation of calcium homeostasis. miR-211-mediated down-regulation of Ezrin leads to Ca2+ influx resulting in the activation of calcineurin, which in turn activates TFEB, the master regulator of lysosomal biogenesis. Light-mediated induction of lysosomal biogenesis and function is impaired in the RPE from miR-211-/- mice that show severely compromised vision. Pharmacological restoration of lysosomal biogenesis through Ezrin inhibition rescued the miR-211-/- phenotype, pointing to a new therapeutic target to counteract retinal degeneration associated with lysosomal dysfunction.

Late-onset mucopolysaccharidosis type IIIA mimicking Usher syndrome.

Mucopolysaccharidosis type IIIA (MPS IIIA or Sanfilippo syndrome type A) is an autosomal recessive lysosomal storage disorder caused by pathogenic variants in the SGSH gene encoding N-sulfoglucosamine sulfohydrolase, an enzyme involved in the degradation of heparan sulfate. MPS IIIA is typically characterized by neurocognitive decline and hepatosplenomegaly with childhood onset. Here, we report on a 53-year-old male subject initially diagnosed with Usher syndrome for the concurrence of retinitis pigmentosa and sensorineural hearing loss. Clinical exome sequencing identified biallelic missense variants in SGSH, and biochemical assays showed complete deficiency of sulfamidase activity and increased urinary glycosaminoglycan excretion. Reverse phenotyping revealed left ventricle pseudo-hypertrophy, hepatosplenomegaly, bilateral deep white matter hyperintensities upon brain MRI, and decreased cortical metabolic activity by PET-CT. On neuropsychological testing, the proband presented only partial and isolated verbal memory deficits. This case illustrates the power of unbiased, comprehensive genetic testing for the diagnosis of challenging mild or atypical forms of MPS IIIA.

Novel and Recurrent Copy Number Variants in ABCA4-Associated Retinopathy.

ABCA4 is the most frequently mutated gene leading to inherited retinal disease (IRD) with over 2200 pathogenic variants reported to date. Of these, ~1% are copy number variants (CNVs) involving the deletion or duplication of genomic regions, typically >50 nucleotides in length. An in-depth assessment of the current literature based on the public database LOVD, regarding the presence of known CNVs and structural variants in ABCA4, and additional sequencing analysis of ABCA4 using single-molecule Molecular Inversion Probes (smMIPs) for 148 probands highlighted recurrent and novel CNVs associated with ABCA4-associated retinopathies. An analysis of the coverage depth in the sequencing data led to the identification of eleven deletions (six novel and five recurrent), three duplications (one novel and two recurrent) and one complex CNV. Of particular interest was the identification of a complex defect, i.e., a 15.3 kb duplicated segment encompassing exon 31 through intron 41 that was inserted at the junction of a downstream 2.7 kb deletion encompassing intron 44 through intron 47. In addition, we identified a 7.0 kb tandem duplication of intron 1 in three cases. The identification of CNVs in ABCA4 can provide patients and their families with a genetic diagnosis whilst expanding our understanding of the complexity of diseases caused by ABCA4 variants.

Definition of the transcriptional units of inherited retinal disease genes by meta-analysis of human retinal transcriptome data.

BACKGROUND: Inherited retinal diseases (IRD) are genetically heterogeneous disorders that cause the dysfunction or loss of photoreceptor cells and ultimately lead to blindness. To date, next-generation sequencing procedures fail to detect pathogenic sequence variants in coding regions of known IRD disease genes in about 30-40% of patients. One of the possible explanations for this missing heritability is the presence of yet unidentified transcripts of known IRD genes. Here, we aimed to define the transcript composition of IRD genes in the human retina by a meta-analysis of publicly available RNA-seq datasets using an ad-hoc designed pipeline. RESULTS: We analysed 218 IRD genes and identified 5,054 transcripts, 3,367 of which were not previously reported. We assessed their putative expression levels and focused our attention on 435 transcripts predicted to account for at least 5% of the expression of the corresponding gene. We looked at the possible impact of the newly identified transcripts at the protein level and experimentally validated a subset of them. CONCLUSIONS: This study provides an unprecedented, detailed overview of the complexity of the human retinal transcriptome that can be instrumental in contributing to the resolution of some cases of missing heritability in IRD patients.

Novel variants in genes related to vesicle-mediated-transport modify Parkinson's disease risk.

OBJECTIVES: VPS35 and VPS13 have been associated with Parkinson's disease (PD), and their shared phenotype in yeast when reduced in function is abnormal vacuolar transport. We aim to test if additional potentially deleterious variants in other genes that share this phenotype can modify the risk for PD. METHODS: 77 VPS and VPS-related genes were analyzed using whole-genome-sequencing data from 202 PD patients of Ashkenazi Jewish (AJ) ancestry. Filtering was done based on quality and functionality scores. Ten variants in nine genes were further genotyped in 1200 consecutively recruited unrelated AJ-PD patients, and allele frequencies and odds ratio calculated compared to gnomAD-AJ-non-neuro database, in un-stratified (n = 1200) and stratified manner (LRRK2-G2019S-PD patients (n = 145), GBA-PD patients (n = 235), and non-carriers of these mutations (NC, n = 787)). RESULTS: Five variants in PIK3C3, VPS11, AP1G2, HGS and VPS13D were significantly associated with PD-risk. PIK3C3-R768W showed a significant association in an un-stratified (all PDs) analysis, as well as in stratified (LRRK2, GBA, and NC) analyses (Odds ratios = 2.71, 5.32, 3.26. and 2.19 with p = 0.0015, 0.002, 0.0287, and 0.0447, respectively). AP1G2-R563W was significantly associated in LRRK2-carriers (OR = 3.69, p = 0.006) while VPS13D-D2932N was significantly associated in GBA-carriers (OR = 5.45, p = 0.0027). VPS11-C846G and HGS-S243Y were significantly associated in NC (OR = 2.48 and 2.06, with p = 0.022 and 0.0163, respectively). CONCLUSIONS: Variants in genes involved in vesicle-mediated protein transport and recycling pathways, including autophagy and mitophagy, may differentially modify PD-risk in LRRK2-carriers, GBA carriers, or NC. Specifically, PIK3C3-R768W is a PD-risk allele, with the highest effect size in LRRK2-G2019S carriers. These results suggest oligogenic effect that may depends on the genetic background of the patient. An unbiased burden of mutations approach in these genes should be evaluated in additional PD and control groups. The mechanisms by which these novel variants interact and increase PD-risk should be researched in depth for better tailoring therapeutic intervention for PD prevention or slowing disease progression.

Patients with DeSanto-Shinawi syndrome: Further extension of phenotype from Italy.

Here we describe three patients with neurodevelopmental disorders characterized by mild-to-moderate intellectual disability, mildly dysmorphic features, and hirsutism, all of which carry de novo sequence variants in the WW domain-containing adaptor of the coiled-coil (WAC) gene; two of these-c.167delA, p.(Asn56I1efs*136) and c.1746G>C, p.(Gln582His)-are novel pathogenic variants, and the third-c.1837C>T, p(Arg613*)-has been previously described. Diseases associated with WAC include DeSanto-Shinawi syndrome; to date, de novo heterozygous constitutional pathogenic WAC variants have caused a syndromic form of intellectual disability and mild dysmorphic features in 33 patients, yet potential associations with other clinical manifestations, such as oligomenorrhea and hyperandrogenism, remain unknown, because the phenotypic spectrum of the condition has not yet been delineated. The patient bearing the novel c.167delA WAC gene variant presented a normal psychomotor development, oligomenorrhea, hyperandrogenism, and hirsutism, and hirsutism was also observed in the patient with the c.1746G>C WAC gene variant. Hypertrichosis and hirsutism have been described in nine DeSanto-Shinawi patients, only in 17 of the 33 aforementioned patients thus far reported this aspect, and no hormonal-pattern data are available. In conclusion, we note that the pathogenic c.167delA WAC variant may be associated with a mild phenotype; and in addition to the neurodevelopmental problems nearly all DeSanto-Shinawi patients experience (i.e., intellectual disability and/or developmental delay), we recommend the addition of mild dysmorphic features, hirsutism, and hypertrichosis to this clinical presentation.

An Alu-mediated duplication in NMNAT1, involved in NAD biosynthesis, causes a novel syndrome, SHILCA, affecting multiple tissues and organs.

We investigated the genetic origin of the phenotype displayed by three children from two unrelated Italian families, presenting with a previously unrecognized autosomal recessive disorder that included a severe form of spondylo-epiphyseal dysplasia, sensorineural hearing loss, intellectual disability and Leber congenital amaurosis (SHILCA), as well as some brain anomalies that were visible at the MRI. Autozygome-based analysis showed that these children shared a 4.76 Mb region of homozygosity on chromosome 1, with an identical haplotype. Nonetheless, whole-exome sequencing failed to identify any shared rare coding variants, in this region or elsewhere. We then determined the transcriptome of patients' fibroblasts by RNA sequencing, followed by additional whole-genome sequencing experiments. Gene expression analysis revealed a 4-fold downregulation of the gene NMNAT1, residing indeed in the shared autozygous interval. Short- and long-read whole-genome sequencing highlighted a duplication involving 2 out of the 5 exons of NMNAT1 main isoform (NM_022787.3), leading to the production of aberrant mRNAs. Pathogenic variants in NMNAT1 have been previously shown to cause non-syndromic Leber congenital amaurosis (LCA). However, no patient with null biallelic mutations has ever been described, and murine Nmnat1 knockouts show embryonic lethality, indicating that complete absence of NMNAT1 activity is probably not compatible with life. The rearrangement found in our cases, presumably causing a strong but not complete reduction of enzymatic activity, may therefore result in an intermediate syndromic phenotype with respect to LCA and lethality.

Pathogenic variants in IMPG1 cause autosomal dominant and autosomal recessive retinitis pigmentosa.

BACKGROUND: Inherited retinal disorders are a clinically and genetically heterogeneous group of conditions and a major cause of visual impairment. Common disease subtypes include vitelliform macular dystrophy (VMD) and retinitis pigmentosa (RP). Despite the identification of over 90 genes associated with RP, conventional genetic testing fails to detect a molecular diagnosis in about one third of patients with RP. METHODS: Exome sequencing was carried out for identifying the disease-causing gene in a family with autosomal dominant RP. Gene panel testing and exome sequencing were performed in 596 RP and VMD families to identified additional IMPG1 variants. In vivo analysis in the medaka fish system by knockdown assays was performed to screen IMPG1 possible pathogenic role. RESULTS: Exome sequencing of a family with RP revealed a splice variant in IMPG1. Subsequently, the same variant was identified in individuals from two families with either RP or VMD. A retrospective study of patients with RP or VMD revealed eight additional families with different missense or nonsense variants in IMPG1. In addition, the clinical diagnosis of the IMPG1 retinopathy-associated variant, originally described as benign concentric annular macular dystrophy, was also revised to RP with early macular involvement. Using morpholino-mediated ablation of Impg1 and its paralog Impg2 in medaka fish, we confirmed a phenotype consistent with that observed in the families, including a decreased length of rod and cone photoreceptor outer segments. CONCLUSION: This study discusses a previously unreported association between monoallelic or biallelic IMPG1 variants and RP. Notably, similar observations have been reported for IMPG2.

Inherited Retinal Diseases Due to RPE65 Variants: From Genetic Diagnostic Management to Therapy.

Inherited retinal diseases (IRDs) are a heterogeneous group of conditions that include retinitis pigmentosa (RP) and Leber congenital amaurosis (LCA) and early-onset severe retinal dystrophy (EO[S]RD), which differ in severity and age of onset. IRDs are caused by mutations in >250 genes. Variants in the RPE65 gene account for 0.6-6% of RP and 3-16% of LCA/EORD cases. Voretigene neparvovec is a gene therapy approved for the treatment of patients with an autosomal recessive retinal dystrophy due to confirmed biallelic RPE65 variants (RPE65-IRDs). Therefore, the accurate molecular diagnosis of RPE65-IRDs is crucial to identify 'actionable' genotypes-i.e., genotypes that may benefit from the treatment-and is an integral part of patient management. To date, hundreds of RPE65 variants have been identified, some of which are classified as pathogenic or likely pathogenic, while the significance of others is yet to be established. In this review, we provide an overview of the genetic diagnostic workup needed to select patients that could be eligible for voretigene neparvovec treatment. Careful clinical characterization of patients by multidisciplinary teams of experts, combined with the availability of next-generation sequencing approaches, can accelerate patients' access to available therapeutic options.

Clinical and Molecular Characterization of Achromatopsia Patients: A Longitudinal Study.

Achromatopsia (ACHM) is a rare genetic disorder of infantile onset affecting cone photoreceptors. To determine the extent of progressive retinal changes in achromatopsia, we performed a detailed longitudinal phenotyping and genetic characterization of an Italian cohort comprising 21 ACHM patients (17 unrelated families). Molecular genetic testing identified biallelic pathogenic mutations in known ACHM genes, including four novel variants. At baseline, the patients presented a reduced best corrected visual acuity (BCVA), reduced macular sensitivity (MS), normal dark-adapted electroretinogram (ERG) responses and undetectable or severely reduced light-adapted ERG. The longitudinal analysis of 16 patients (mean follow-up: 5.4 ± 1.0 years) showed a significant decline of BCVA (0.012 logMAR/year) and MS (-0.16 dB/year). Light-adapted and flicker ERG responses decreased below noise level in three and two patients, respectively. Only two patients (12.5%) progressed to a worst OCT grading during the follow-up. Our findings corroborate the notion that ACHM is a progressive disease in terms of BCVA, MS and ERG responses, and affects slowly the structural integrity of the retina. These observations can serve towards the development of guidelines for patient selection and intervention timing in forthcoming gene replacement therapies.

Resolving the dark matter of ABCA4 for 1054 Stargardt disease probands through integrated genomics and transcriptomics.

PURPOSE: Missing heritability in human diseases represents a major challenge, and this is particularly true for ABCA4-associated Stargardt disease (STGD1). We aimed to elucidate the genomic and transcriptomic variation in 1054 unsolved STGD and STGD-like probands. METHODS: Sequencing of the complete 128-kb ABCA4 gene was performed using single-molecule molecular inversion probes (smMIPs), based on a semiautomated and cost-effective method. Structural variants (SVs) were identified using relative read coverage analyses and putative splice defects were studied using in vitro assays. RESULTS: In 448 biallelic probands 14 known and 13 novel deep-intronic variants were found, resulting in pseudoexon (PE) insertions or exon elongations in 105 alleles. Intriguingly, intron 13 variants c.1938-621G>A and c.1938-514G>A resulted in dual PE insertions consisting of the same upstream, but different downstream PEs. The intron 44 variant c.6148-84A>T resulted in two PE insertions and flanking exon deletions. Eleven distinct large deletions were found, two of which contained small inverted segments. Uniparental isodisomy of chromosome 1 was identified in one proband. CONCLUSION: Deep sequencing of ABCA4 and midigene-based splice assays allowed the identification of SVs and causal deep-intronic variants in 25% of biallelic STGD1 cases, which represents a model study that can be applied to other inherited diseases.

The Pervasive Role of the miR-181 Family in Development, Neurodegeneration, and Cancer.

MicroRNAs (miRNAs) are small noncoding RNAs playing a fundamental role in the regulation of gene expression. Evidence accumulating in the past decades indicate that they are capable of simultaneously modulating diverse signaling pathways involved in a variety of pathophysiological processes. In the present review, we provide a comprehensive overview of the function of a highly conserved group of miRNAs, the miR-181 family, both in physiological as well as in pathological conditions. We summarize a large body of studies highlighting a role for this miRNA family in the regulation of key biological processes such as embryonic development, cell proliferation, apoptosis, autophagy, mitochondrial function, and immune response. Importantly, members of this family have been involved in many pathological processes underlying the most common neurodegenerative disorders as well as different solid tumors and hematological malignancies. The relevance of this miRNA family in the pathogenesis of these disorders and their possible influence on the severity of their manifestations will be discussed. A better understanding of the miR-181 family in pathological conditions may open new therapeutic avenues for devasting disorders such as neurodegenerative diseases and cancer.

MIB2 variants altering NOTCH signalling result in left ventricle hypertrabeculation/non-compaction and are associated with Ménétrier-like gastropathy.

We performed whole exome sequencing in individuals from a family with autosomal dominant gastropathy resembling Ménétrier disease, a premalignant gastric disorder with epithelial hyperplasia and enhanced EGFR signalling. Ménétrier disease is believed to be an acquired disorder, but its aetiology is unknown. In affected members, we found a missense p.V742G variant in MIB2, a gene regulating NOTCH signalling that has not been previously linked to human diseases. The variant segregated with the disease in the pedigree, affected a highly conserved amino acid residue, and was predicted to be deleterious although it was found with a low frequency in control individuals. The purified protein carrying the p.V742G variant showed reduced ubiquitination activity in vitro and white blood cells from affected individuals exhibited significant reductions of HES1 and NOTCH3 expression reflecting alteration of NOTCH signalling. Because mutations of MIB1, the homolog of MIB2, have been found in patients with left ventricle non-compaction (LVNC), we investigated members of our family with Ménétrier-like disease for this cardiac abnormality. Asymptomatic left ventricular hypertrabeculation, the mildest end of the LVNC spectrum, was detected in two members carrying the MIB2 variant. Finally, we identified an additional MIB2 variant (p.V984L) affecting protein stability in an unrelated isolated case with LVNC. Expression of both MIB2 variants affected NOTCH signalling, proliferation and apoptosis in primary rat cardiomyocytes.In conclusion, we report the first example of left ventricular hypertrabeculation/LVNC with germline MIB2 variants resulting in altered NOTCH signalling that might be associated with a gastropathy clinically overlapping with Ménétrier disease.

Non-coding RNAs in retinal development and function.

Accumulating evidence on the role of non-protein-coding RNA sequences in the regulation of gene expression is greatly expanding our understanding of the flow of genetic information within biological systems. The interplay between protein-coding and non-coding RNAs (ncRNAs) is essential for tissue development, homeostasis, and function. NcRNAs can be divided in short ncRNAs, whose main subtype is represented by microRNAs, and long ncRNAs, which constitute a more heterogeneous class. The retina is a light-sensitive tissue consisting of highly interconnected cell types and is the primary target of many genetic diseases. Among these, the genetically heterogeneous group of inherited retinal diseases (IRDs) represents the most frequent monogenic cause of visual impairment that can ultimately lead to blindness. Here, we provide an overview on the role of ncRNAs in retinal development and function with an emphasis on microRNAs and on different types of long ncRNAs. We also review how sequence variations in ncRNAs can play a pathogenic role in IRDs as well as in multifactorial ocular disorders. These data indicate that a comprehensive study of the contribution of ncRNAs to the mutation repertoire associated with retinal disease can shed light on previously unknown pathophysiological mechanisms and open new therapeutic avenues. We conclude that a more comprehensive dissection of the pathogenic role of non-coding RNAs in retinal function and disease will not only improve our diagnostic ability, but will allow the development of novel targeted therapies for ocular disease.

microRNAs as biomarkers in Pompe disease.

PURPOSE: We studied microRNAs as potential biomarkers for Pompe disease. METHODS: We analyzed microRNA expression by small RNA-seq in tissues from the disease murine model at two different ages (3 and 9 months), and in plasma from Pompe patients. RESULTS: In the mouse model we found 211 microRNAs that were differentially expressed in gastrocnemii and 66 in heart, with a different pattern of expression at different ages. In a preliminary analysis in plasma from six patients 55 microRNAs were differentially expressed. Sixteen of these microRNAs were common to those dysregulated in mouse tissues. These microRNAs are known to modulate the expression of genes involved in relevant pathways for Pompe disease pathophysiology (autophagy, muscle regeneration, muscle atrophy). One of these microRNAs, miR-133a, was selected for further quantitative real-time polymerase chain reaction analysis in plasma samples from 52 patients, obtained from seven Italian and Dutch biobanks. miR-133a levels were significantly higher in Pompe disease patients than in controls and correlated with phenotype severity, with higher levels in infantile compared with late-onset patients. In three infantile patients miR-133a decreased after start of enzyme replacement therapy and evidence of clinical improvement. CONCLUSION: Circulating microRNAs may represent additional biomarkers of Pompe disease severity and of response to therapy.

Correction: Biallelic sequence and structural variants in RAX2 are a novel cause for autosomal recessive inherited retinal disease.

The original version of this Article contained an incorrect version of Fig. 3, which included two variants initially shown in black text in Fig. 3a that the authors removed from the final manuscript. The correct version of Fig. 3 without the two variants now appears in the PDF and HTML versions of the Article.

Biallelic sequence and structural variants in RAX2 are a novel cause for autosomal recessive inherited retinal disease.

PURPOSE: RAX2 encodes a homeobox-containing transcription factor, in which four monoallelic pathogenic variants have been described in autosomal dominant cone-dominated retinal disease. METHODS: Exome sequencing in a European cohort with inherited retinal disease (IRD) (n = 2086) was combined with protein structure modeling of RAX2 missense variants, bioinformatics analysis of deletion breakpoints, haplotyping of RAX2 variant c.335dup, and clinical assessment of biallelic RAX2-positive cases and carrier family members. RESULTS: Biallelic RAX2 sequence and structural variants were found in five unrelated European index cases, displaying nonsyndromic autosomal recessive retinitis pigmentosa (ARRP) with an age of onset ranging from childhood to the mid-40s (average mid-30s). Protein structure modeling points to loss of function of the novel recessive missense variants and to a dominant-negative effect of the reported dominant RAX2 alleles. Structural variants were fine-mapped to disentangle their underlying mechanisms. Haplotyping of c.335dup in two cases suggests a common ancestry. CONCLUSION: This study supports a role for RAX2 as a novel disease gene for recessive IRD, broadening the mutation spectrum from sequence to structural variants and revealing a founder effect. The identification of biallelic RAX2 pathogenic variants in five unrelated families shows that RAX2 loss of function may be a nonnegligible cause of IRD in unsolved ARRP cases.

Clinical and Genetic Analysis of a European Cohort with Pericentral Retinitis Pigmentosa.

Retinitis pigmentosa (RP) is a clinically heterogenous disease that comprises a wide range of phenotypic and genetic subtypes. Pericentral RP is an atypical form of RP characterized by bone-spicule pigmentation and/or atrophy confined in the near mid-periphery of the retina. In contrast to classic RP, the far periphery is better preserved in pericentral RP. The aim of this study was to perform the first detailed clinical and genetic analysis of a cohort of European subjects with pericentral RP to determine the phenotypic features and the genetic bases of the disease. A total of 54 subjects from 48 independent families with pericentral RP, non-syndromic and syndromic, were evaluated through a full ophthalmological examination and underwent clinical exome or retinopathy gene panel sequencing. Disease-causative variants were identified in 22 of the 35 families (63%) in 10 different genes, four of which are also responsible for syndromic RP. Thirteen of the 34 likely pathogenic variants were novel. Intra-familiar variability was also observed. The current study confirms the mild phenotype of pericentral RP and extends the spectrum of genes associated with this condition.

High-resolution analysis of the human retina miRNome reveals isomiR variations and novel microRNAs.

MicroRNAs play a fundamental role in retinal development and function. To characterise the miRNome of the human retina, we carried out deep sequencing analysis on sixteen individuals. We established the catalogue of retina-expressed miRNAs, determined their relative abundance and found that a small number of miRNAs accounts for almost 90% of the retina miRNome. We discovered more than 3000 miRNA variants (isomiRs), encompassing a wide range of sequence variations, which include seed modifications that are predicted to have an impact on miRNA action. We demonstrated that a seed-modifying isomiR of the retina-enriched miR-124-3p was endowed with different targeting properties with respect to the corresponding canonical form. Moreover, we identified 51 putative novel, retina-specific miRNAs and experimentally validated the expression for nine of them. Finally, a parallel analysis of the human Retinal Pigment Epithelium (RPE)/choroid, two tissues that are known to be crucial for retina homeostasis, yielded notably distinct miRNA enrichment patterns compared to the retina. The generated data are accessible through an ad hoc database. This study is the first to reveal the complexity of the human retina miRNome at nucleotide resolution and constitutes a unique resource to assess the contribution of miRNAs to the pathophysiology of the human retina.