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.

Background: The most frequent mutation of Friedreich ataxia (FRDA) is the abnormal expansion of a GAA repeat located within the first intron of FXN gene. It is known that the length of GAA is directly correlated with disease severity. The effect of mutation is a severe reduction of mRNA. Recently, a link among aberrant CpG methylation, chromatin organisation and GAA repeat was proposed. Methods: In this study, using pyrosequencing technology, we have performed a quantitative analysis of the methylation status of five CpG sites located within the region upstream of GAA repeat, in 67 FRDA patients. Results: We confirm previous observation about differences in the methylation degree between FRDA individuals and controls. We showed a direct correlation between CpG methylation and triplet expansion size. Significant differences were found for each CpG tested (ANOVA p<0.001). These differences were largest for CpG1 and CpG2: 84.45% and 76.80%, respectively, in FRDA patients compared to 19.65% and 23.34% in the controls. Most importantly, we found a strong inverse correlation between CpG2 methylation degree and age of onset (Spearman’s ρ = −0.550, p<0.001). Conclusion: Because epigenetic changes may cause or contribute to gene silencing, our data may have relevance for the therapeutic approach to FRDA. Since the analysis can be performed in peripheral blood leucocytes (PBL), evaluation of the methylation status of specific CpG sites in FRDA patients could be a convenient biomarker.

The human retina is a specialized tissue involved in light stimulus transduction. Despite its unique biology, an accurate reference transcriptome is still missing. Here, we performed gene expression analysis (RNA-seq) of 50 retinal samples from non-visually impaired post-mortem donors. We identified novel transcripts with high confidence (Observed Transcriptome (ObsT)) and quantified the expression level of known transcripts (Reference Transcriptome (RefT)). The ObsT included 77 623 transcripts (23 960 genes) covering 137 Mb (35 Mb new transcribed genome). Most of the transcripts (92%) were multi-exonic: 81% with known isoforms, 16% with new isoforms and 3% belonging to new genes. The RefT included 13 792 genes across 94 521 known transcripts. Mitochondrial genes were among the most highly expressed, accounting for about 10% of the reads. Of all the protein-coding genes in Gencode, 65% are expressed in the retina. We exploited inter-individual variability in gene expression to infer a gene co-expression network and to identify genes specifically expressed in photoreceptor cells. We experimentally validated the photoreceptors localization of three genes in human retina that had not been previously reported. RNA-seq data and the gene co-expression network are available online (http://retina.tigem.it).

Lipedema is an often underdiagnosed chronic disorder that affects subcutaneous adipose tissue almost exclusively in women, which leads to disproportionate fat accumulation in the lower and upper body extremities. Common comorbidities include anxiety, depression, and pain. The correlation between mood disorder and subcutaneous fat deposition suggests the involvement of steroids metabolism and neurohormones signaling, however no clear association has been established so far. In this study, we report on a family with three patients affected by sex-limited autosomal dominant nonsyndromic lipedema. They had been screened by whole exome sequencing (WES) which led to the discovery of a missense variant p.(Leu213Gln) in AKR1C1, the gene encoding for an aldo-keto reductase catalyzing the reduction of progesterone to its inactive form, 20-α-hydroxyprogesterone. Comparative molecular dynamics simulations of the wild-type vs. variant enzyme, corroborated by a thorough structural and functional bioinformatic analysis, suggest a partial loss-of-function of the variant. This would result in a slower and less efficient reduction of progesterone to hydroxyprogesterone and an increased subcutaneous fat deposition in variant carriers. Overall, our results suggest that AKR1C1 is the first candidate gene associated with nonsyndromic lipedema.