The first genetic landscape of inherited retinal dystrophies in Portuguese patients identifies recurrent homozygous mutations as a frequent cause of pathogenesis.

Inherited retinal diseases (IRDs) are a group of ocular conditions characterized by an elevated genetic and clinical heterogeneity. They are transmitted almost invariantly as monogenic traits. However, with more than 280 disease genes identified so far, association of clinical phenotypes with genotypes can be very challenging, and molecular diagnosis is essential for genetic counseling and correct management of the disease. In addition, the prevalence and the assortment of IRD mutations are often population-specific. In this work, we examined 230 families from Portugal, with individuals suffering from a variety of IRD diagnostic classes (270 subjects in total). Overall, we identified 157 unique mutations (34 previously unreported) in 57 distinct genes, with a diagnostic rate of 76%. The IRD mutational landscape was, to some extent, different from those reported in other European populations, including Spanish cohorts. For instance, the EYS gene appeared to be the most frequently mutated, with a prevalence of 10% among all IRD cases. This was, in part, due to the presence of a recurrent and seemingly founder mutation involving the deletion of exons 13 and 14 of this gene. Moreover, our analysis highlighted that as many as 51% of our cases had mutations in a homozygous state. To our knowledge, this is the first study assessing a cross-sectional genotype-phenotype landscape of IRDs in Portugal. Our data reveal a rather unique distribution of mutations, possibly shaped by a small number of rare ancestral events that have now become prevalent alleles in patients.

Identification of New Vulnerabilities in Conjunctival Melanoma Using Image-Based High Content Drug Screening.

Recent evidence suggests that numerous similarities exist between the genomic landscapes of both conjunctival and cutaneous melanoma. Since alterations of several components of the MAP kinases, PI3K/mTOR, and cell cycle pathways have been reported in conjunctival melanoma, we decided to assess the sensitivity of conjunctival melanoma to targeted inhibition mostly of kinase inhibitors. A high content drug screening assay based on automated fluorescence microscopy was performed in three conjunctival melanoma cell lines with different genomic backgrounds with 489 kinase inhibitors and 53 other inhibitors. IC50 and apoptosis induction were respectively assessed for 53 and 48 compounds. The genomic background influenced the response to MAK and PI3K/mTOR inhibition, more specifically cell lines with BRAF V600E mutations were more sensitive to BRAF/MEK inhibition, while CRMM2 bearing the NRASQ61L mutation was more sensitive to PI3k/mTOR inhibition. All cell lines demonstrated sensitivity to cell cycle inhibition, being more pronounced in CRMM2, especially with polo-like inhibitors. Our data also revealed new vulnerabilities to Hsp90 and Src inhibition. This study demonstrates that the genomic background partially influences the response to targeted therapy and uncovers a large panel of potential vulnerabilities in conjunctival melanoma that may expand available options for the management of this tumor.

Analysis of missense variants in the human genome reveals widespread gene-specific clustering and improves prediction of pathogenicity.

We used a machine learning approach to analyze the within-gene distribution of missense variants observed in hereditary conditions and cancer. When applied to 840 genes from the ClinVar database, this approach detected a significant non-random distribution of pathogenic and benign variants in 387 (46%) and 172 (20%) genes, respectively, revealing that variant clustering is widespread across the human exome. This clustering likely occurs as a consequence of mechanisms shaping pathogenicity at the protein level, as illustrated by the overlap of some clusters with known functional domains. We then took advantage of these findings to develop a pathogenicity predictor, MutScore, that integrates qualitative features of DNA substitutions with the new additional information derived from this positional clustering. Using a random forest approach, MutScore was able to identify pathogenic missense mutations with very high accuracy, outperforming existing predictive tools, especially for variants associated with autosomal-dominant disease and cancer. Thus, the within-gene clustering of pathogenic and benign DNA changes is an important and previously underappreciated feature of the human exome, which can be harnessed to improve the prediction of pathogenicity and disambiguation of DNA variants of uncertain significance.

A monoallelic SEC23A variant E599K associated with cranio-lenticulo-sutural dysplasia.

Cranio-lenticulo-sutural dysplasia (CLSD; MIM 607812) is a rare or underdiagnosed condition, as only two families have been reported. The original family (Boyadjiev et al., Human Genetics, 2003, 113, 1-9 and Boyadjiev et al., Nature Genetics, 2006, 38, 1192-1197) showed recessive inheritance of the condition with a biallelic SEC23A missense variant in affected individuals. In contrast, another child with sporadic CLSD had a monoallelic SEC23A variant inherited from the reportedly unaffected father (Boyadjiev et al., Clinical Genetics, 2011, 80, 169-176), raising questions on possible digenism. Here, we report a 2-month-old boy seen because of large fontanels with wide cranial sutures, a large forehead, hypertelorism, a thin nose, a high arched palate, and micrognathia. His mother was clinically unremarkable, while his father had a history of large fontanels in infancy who had closed only around age 10 years; he also had a large forehead, hypertelorism, a thin, beaked nose and was operated for bilateral glaucoma with exfoliation of the lens capsule. Trio genome sequencing and familial segregation revealed a monoallelic c.1795G > A transition in SEC23A that was de novo in the father and transmitted to the proband. The variant predicts a nonconservative substitution (p.E599K) in an ultra-conserved residue that is seen in 3D models of yeast SEC23 to be involved in direct binding between SEC23 and SAR1 subunits of the coat protein complex II coat. This observation confirms the link between SEC23A variants and CLSD but suggests that in addition to the recessive inheritance described in the original family, SEC23A variants may result in dominant inheritance of CLSD, possibly by a dominant-negative disruptive effect on the SEC23 multimer.

De novo variants in CACNA1E found in patients with intellectual disability, developmental regression and social cognition deficit but no seizures.

BACKGROUND: De novo variants in the voltage-gated calcium channel subunit α1 E gene (CACNA1E) have been described as causative of epileptic encephalopathy with contractures, macrocephaly and dyskinesias. METHODS: Following the observation of an index patient with developmental delay and autism spectrum disorder (ASD) without seizures who had a de novo deleterious CACNA1E variant, we screened GeneMatcher for other individuals with CACNA1E variants and neurodevelopmental phenotypes without epilepsy. The spectrum of pathogenic CACNA1E variants was compared to the mutational landscape of variants in the gnomAD control population database. RESULTS: We identified seven unrelated individuals with intellectual disability, developmental regression and ASD-like behavioral profile, and notably without epilepsy, who had de novo heterozygous putatively pathogenic variants in CACNA1E. Age of onset of clinical manifestation, presence or absence of regression and degree of severity were variable, and no clear-cut genotype-phenotype association could be recognized. The analysis of disease-associated variants and their comparison to benign variants from the control population allowed for the identification of regions in the CACNA1E protein that seem to be intolerant to substitutions and thus more likely to harbor pathogenic variants. As in a few reported cases with CACNA1E variants and epilepsy, one patient showed a positive clinical behavioral response to topiramate, a specific calcium channel modulator. LIMITATIONS: The significance of our study is limited by the absence of functional experiments of the effect of identified variants, the small sample size and the lack of systematic ASD assessment in all participants. Moreover, topiramate was given to one patient only and for a short period of time. CONCLUSIONS: Our results indicate that CACNA1E variants may result in neurodevelopmental disorders without epilepsy and expand the mutational and phenotypic spectrum of this gene. CACNA1E deserves to be included in gene panels for non-specific developmental disorders, including ASD, and not limited to patients with seizures, to improve diagnostic recognition and explore the possible efficacy of topiramate.

CNV Detection from Exome Sequencing Data in Routine Diagnostics of Rare Genetic Disorders: Opportunities and Limitations.

To assess the potential of detecting copy number variations (CNVs) directly from exome sequencing (ES) data in diagnostic settings, we developed a CNV-detection pipeline based on ExomeDepth software and applied it to ES data of 450 individuals. Initially, only CNVs affecting genes in the requested diagnostic gene panels were scored and tested against arrayCGH results. Pathogenic CNVs were detected in 18 individuals. Most detected CNVs were larger than 400 kb (11/18), but three individuals had small CNVs impacting one or a few exons only and were thus not detectable by arrayCGH. Conversely, two pathogenic CNVs were initially missed, as they impacted genes not included in the original gene panel analysed, and a third one was missed as it was in a poorly covered region. The overall combined diagnostic rate (SNVs + CNVs) in our cohort was 36%, with wide differences between clinical domains. We conclude that (1) the ES-based CNV pipeline detects efficiently large and small pathogenic CNVs, (2) the detection of CNV relies on uniformity of sequencing and good coverage, and (3) in patients who remain unsolved by the gene panel analysis, CNV analysis should be extended to all captured genes, as diagnostically relevant CNVs may occur everywhere in the genome.

A live single-cell reporter assay links intratumor heterogeneity to metastatic proclivity in Ewing sarcoma.

Targeting of the most aggressive tumor cell subpopulations is key for effective management of most solid malignancies. However, the metastable nature of tumor heterogeneity, which allows cells to transition between strong and weak tumorigenic phenotypes, and the lack of reliable markers of tumor-promoting properties hamper identification of the most relevant cells. To overcome these obstacles, we designed a functional microRNA (miR)-based live-cell reporter assay to identify highly tumorigenic cells in xenotransplants of primary Ewing sarcoma (EwS) 3D cultures. Leveraging the inverse relationship between cell pluripotency and miR-145 expression, we successfully separated highly tumorigenic, metastasis-prone (miR-145low) cells from poorly tumorigenic, nonmetastatic (miR-145high) counterparts. Gene expression and functional studies of the two cell populations identified the EPHB2 receptor as a prognostic biomarker in patients with EwS and a major promoter of metastasis. Our study provides a simple and powerful means to identify and isolate tumor cells that display aggressive behavior.

A hypomorphic variant in EYS detected by genome-wide association study contributes toward retinitis pigmentosa.

The genetic basis of Japanese autosomal recessive retinitis pigmentosa (ARRP) remains largely unknown. Herein, we applied a 2-step genome-wide association study (GWAS) in 640 Japanese patients. Meta-GWAS identified three independent peaks at P < 5.0 × 10-8, all within the major ARRP gene EYS. Two of the three were each in linkage disequilibrium with a different low frequency variant (allele frequency < 0.05); a known founder Mendelian mutation (c.4957dupA, p.S1653Kfs*2) and a non-synonymous variant (c.2528 G > A, p.G843E) of unknown significance. mRNA harboring c.2528 G > A failed to restore rhodopsin mislocalization induced by morpholino-mediated knockdown of eys in zebrafish, consistent with the variant being pathogenic. c.2528 G > A solved an additional 7.0% of Japanese ARRP cases. The third peak was in linkage disequilibrium with a common non-synonymous variant (c.7666 A > T, p.S2556C), possibly representing an unreported disease-susceptibility signal. GWAS successfully unraveled genetic causes of a rare monogenic disorder and identified a high frequency variant potentially linked to development of local genome therapeutics.

AutoMap is a high performance homozygosity mapping tool using next-generation sequencing data.

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.

Genomic and transcriptomic landscape of conjunctival melanoma.

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.

A frequent variant in the Japanese population determines quasi-Mendelian inheritance of rare retinal ciliopathy.

Hereditary retinal degenerations (HRDs) are Mendelian diseases characterized by progressive blindness and caused by ultra-rare mutations. In a genomic screen of 331 unrelated Japanese patients, we identify a disruptive Alu insertion and a nonsense variant (p.Arg1933*) in the ciliary gene RP1, neither of which are rare alleles in Japan. p.Arg1933* is almost polymorphic (frequency = 0.6%, amongst 12,000 individuals), does not cause disease in homozygosis or heterozygosis, and yet is significantly enriched in HRD patients (frequency = 2.1%, i.e., a 3.5-fold enrichment; p-value = 9.2 × 10-5). Familial co-segregation and association analyses show that p.Arg1933* can act as a Mendelian mutation in trans with the Alu insertion, but might also associate with disease in combination with two alleles in the EYS gene in a non-Mendelian pattern of heredity. Our results suggest that rare conditions such as HRDs can be paradoxically determined by relatively common variants, following a quasi-Mendelian model linking monogenic and complex inheritance.

DOMINO: Using Machine Learning to Predict Genes Associated with Dominant Disorders.

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.

Comprehensive Genetic Landscape of Uveal Melanoma by Whole-Genome Sequencing.

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.

Mutations in CEP78 Cause Cone-Rod Dystrophy and Hearing Loss Associated with Primary-Cilia Defects.

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.

A missense mutation in the rabbit melanocortin 4 receptor (MC4R) gene is associated with finishing weight in a meat rabbit line.

In this study we resequenced 1729 bp of the rabbit melanocortin 4 receptor (MC4 R) gene in 31 rabbits from different breeds/lines and identified ten polymorphisms: one was an indel and 9 were single nucleotide polymorphisms (SNPs). The indel and 5 SNPs were in the 5'-flanking region, 3 were synonymous SNPs and one was a missense mutation (c.101G>A; p.G34D), located in a conserved position of the extracellular tail of the MC4 R protein. The missense mutation was analyzed in a panel of 74 rabbits of different breeds and in 516 performance tested rabbits of a commercial paternal line under selection for growth efficiency. Association analysis indicated that rabbits with the less frequent genotype in this population (DD) had a lighter weight at 70 postnatal days than animals with genotype GD (P < 0.10) and animals with genotype GG (P < 0.05). This is the third study on candidate genes, after those on GH1 and IGF2 that reported a marker associated with finishing weight. Therefore, it seems that a candidate gene approach in rabbit based on previous information accumulated in other livestock species could be useful to identify genes explaining a fraction of variability of performance traits with potential application on rabbit breeding and selection.