Uveal melanoma (UM) is a rare intraocular tumor that, similar to cutaneous melanoma, originates from melanocytes. To gain insights into its genetics, we performed whole-genome sequencing at very deep coverage of tumor-control pairs in 33 samples (24 primary and 9 metastases). Genome-wide, the number of coding mutations was rather low (only 17 variants per tumor on average; range 7-28), thus radically different from cutaneous melanoma, where hundreds of exonic DNA insults are usually detected. Furthermore, no UV light-induced mutational signature was identified. Recurrent coding mutations were found in the known UM drivers GNAQ, GNA11, BAP1, EIF1AX, and SF3B1. Other genes, i.e., TP53BP1, CSMD1, TTC28, DLK2, and KTN1, were also found to harbor somatic mutations in more than one individual, possibly indicating a previously undescribed association with UM pathogenesis. De novo assembly of unmatched reads from non-coding DNA revealed peculiar copy-number variations defining specific UM subtypes, which in turn could be associated with metastatic transformation. Mutational-driven comparison with other tumor types showed that UM is very similar to pediatric tumors, characterized by very few somatic insults and, possibly, important epigenetic changes. Through the analysis of whole-genome sequencing data, our findings shed new light on the molecular genetics of uveal melanoma, delineating it as an atypical tumor of the adult for which somatic events other than mutations in exonic DNA shape its genetic landscape and define its metastatic potential.

In contrast to recessive conditions with biallelic inheritance, identification of dominant (monoallelic) mutations for Mendelian disorders is more difficult, because of the abundance of benign heterozygous variants that act as massive background noise (typically, in a 400:1 excess ratio). To reduce this overflow of false positives in next-generation sequencing (NGS) screens, we developed DOMINO, a tool assessing the likelihood for a gene to harbor dominant changes. Unlike commonly-used predictors of pathogenicity, DOMINO takes into consideration features that are the properties of genes, rather than of variants. It uses a machine-learning approach to extract discriminant information from a broad array of features (N ¼ 432), including: genomic data, intra-, and interspecies conservation, gene expression, protein-protein interactions, protein structure, etc. DOMINO's iterative architecture includes a training process on 985 genes with well-established inheritance patterns for Mendelian conditions, and repeated cross-validation that optimizes its discriminant power. When validated on 99 newly-discovered genes with pathogenic mutations, the algorithm displays an excellent final performance, with an area under the curve (AUC) of 0.92. Furthermore, unsupervised analysis by DOMINO of real sets of NGS data from individuals with intellectual disability or epilepsy correctly recognizes known genes and predicts 9 new candidates, with very high confidence. In summary, DOMINO is a robust and reliable tool that can infer dominance of candidate genes with high sensitivity and specificity, making it a useful complement to any NGS pipeline dealing with the analysis of the morbid human genome.

Cone-rod degeneration (CRD) belongs to the disease spectrum of retinal degenerations, a group of hereditary disorders characterized by an extreme clinical and genetic heterogeneity. It mainly differentiates from other retinal dystrophies, and in particular from the more frequent disease retinitis pigmentosa, because cone photoreceptors degenerate at a higher rate than rod photoreceptors, causing severe deficiency of central vision. After exome analysis of a cohort of individuals with CRD, we identified biallelic mutations in the orphan gene CEP78 in three subjects from two families: one from Greece and another from Sweden. The Greek subject, from the island of Crete, was homozygous for the c.499+1G>T (IVS3+1G>T) mutation in intron 3. The Swedish subjects, two siblings, were compound heterozygotes for the nearby mutation c.499+5G>A (IVS3+5G>A) and for the frameshift-causing variant c.633delC (p.Trp212Glyfs(∗)18). In addition to CRD, these three individuals had hearing loss or hearing deficit. Immunostaining highlighted the presence of CEP78 in the inner segments of retinal photoreceptors, predominantly of cones, and at the base of the primary cilium of fibroblasts. Interaction studies also showed that CEP78 binds to FAM161A, another ciliary protein associated with retinal degeneration. Finally, analysis of skin fibroblasts derived from affected individuals revealed abnormal ciliary morphology, as compared to that of control cells. Altogether, our data strongly suggest that mutations in CEP78 cause a previously undescribed clinical entity of a ciliary nature characterized by blindness and deafness but clearly distinct from Usher syndrome, a condition for which visual impairment is due to retinitis pigmentosa.

Homozygosity mapping is a powerful method for identifying mutations in patients with recessive conditions, especially in consanguineous families or isolated populations. Historically, it has been used in conjunction with genotypes from highly polymorphic markers, such as DNA microsatellites or common SNPs. Traditional software performs rather poorly with data from Whole Exome Sequencing (WES) and Whole Genome Sequencing (WGS), which are now extensively used in medical genetics. We develop AutoMap, a tool that is both web-based or downloadable, to allow performing homozygosity mapping directly on VCF (Variant Call Format) calls from WES or WGS projects. Following a training step on WES data from 26 consanguineous families and a validation procedure on a matched cohort, our method shows higher overall performances when compared with eight existing tools. Most importantly, when tested on real cases with negative molecular diagnosis from an internal set, AutoMap detects three gene-disease and multiple variant-disease associations that were previously unrecognized, projecting clear benefits for both molecular diagnosis and research activities in medical genetics.

Conjunctival melanoma (CJM) is a rare but potentially lethal and highly-recurrent cancer of the eye. Similar to cutaneous melanoma (CM), it originates from melanocytes. Unlike CM, however, CJM is relatively poorly characterized from a genomic point of view. To fill this knowledge gap and gain insight into the genomic nature of CJM, we performed whole-exome (WES) or whole-genome sequencing (WGS) of tumor-normal tissue pairs in 14 affected individuals, as well as RNA sequencing in a subset of 11 tumor tissues. Our results show that, similarly to CM, CJM is also characterized by a very high mutation load, composed of approximately 500 somatic mutations in exonic regions. This, as well as the presence of a UV light-induced mutational signature, are clear signs of the role of sunlight in CJM tumorigenesis. In addition, the genomic classification of CM proposed by TCGA seems to be well-applicable to CJM, with the presence of four typical subclasses defined on the basis of the most frequently mutated genes: BRAF, NF1, RAS, and triple wild-type. In line with these results, transcriptomic analyses revealed similarities with CM as well, namely the presence of a transcriptomic subtype enriched for immune genes and a subtype enriched for genes associated with keratins and epithelial functions. Finally, in seven tumors we detected somatic mutations in ACSS3, a possible new candidate oncogene. Transfected conjunctival melanoma cells overexpressing mutant ACSS3 showed higher proliferative activity, supporting the direct involvement of this gene in the tumorigenesis of CJM. Altogether, our results provide the first unbiased and complete genomic and transcriptomic classification of CJM.