Whole exome sequencing of low grade serous ovarian carcinoma identifies genomic events associated with clinical outcome.

OBJECTIVES: Low-grade serous ovarian carcinoma (LGSOC) is a distinct, rare, ovarian cancer type characterised by younger patient age and intrinsic chemoresistance. Understanding the molecular landscape is crucial for optimising targeted therapy. METHODS: Genomic data from whole exome sequencing of tumour tissue was analysed in a LGSOC cohort with detailed clinical annotation. RESULTS: 63 cases were analysed and three subgroups identified based on single nucleotide variants: canonical MAPK mutant (cMAPKm: 52%, KRAS/BRAF/NRAS), MAPK-associated gene mutation (MAPK-assoc: 27%) and MAPK wild-type (MAPKwt: 21%). NOTCH pathway disruption occurred across all subgroups. Tumour mutational burden (TMB), mutational signatures and recurrent copy number (CN) changes varied across the cohort with co-occurrence of chromosome 1p loss and 1q gain (CN Chr1pq) a recurrent feature. Low TMB and CN Chr1pq were associated with inferior disease-specific survival (HR 6.43; p < 0.001 and HR 3.29, p = 0.011 respectively). Stepwise genomic classification in relation to outcome resulted in four groups (TMB low; CN Chr1pq; MAPKwt/MAPKassoc; cMAPKm). 5 year disease-specific survival was 46%, 55%, 79% and 100% respectively for these groups. The two most favourable genomic subgroups were enriched for the SBS10b mutational signature, particularly the cMAPKm subgroup. CONCLUSIONS: LGSOC comprises multiple genomic subgroups with distinct clinical and molecular features. Chr1pq CN arm disruption and TMB represent promising methods to identify individuals with poorer prognosis. Further investigation of the molecular basis for these observations is required. MAPKwt cases represent around a fifth of patients. NOTCH inhibitors represent a candidate therapeutic strategy worthy of exploration across these cases.

Clinical and molecular characterization of ovarian carcinoma displaying isolated lymph node relapse.

BACKGROUND: Disease relapse is the primary cause of death from ovarian carcinoma. Isolated lymph node relapse is a rare pattern of ovarian carcinoma recurrence, with a reported median postrelapse survival of 2.5 to 4 years. To date, investigations have not compared isolated lymph node relapse ovarian carcinoma directly to a matched extranodal relapse cohort or performed molecular characterization of cases that subsequently experience isolated lymph node relapse. OBJECTIVE: Here we seek to compare the clinical outcome, tumor-infiltrating lymphocyte burden, and frequency of known prognostic genomic events in isolated lymph node relapse ovarian carcinoma vs extranodal relapse ovarian carcinoma. STUDY DESIGN: Forty-nine isolated lymph node relapse ovarian carcinoma patients were identified and matched to 49 extranodal relapse cases using the Edinburgh Ovarian Cancer Database, from which the clinical data for identified patients were retrieved. Matching criteria were disease stage, histologic subtype and grade, extent of residual disease following surgical debulking, and age at diagnosis. Clinicopathologic factors and survival data were compared between the isolated lymph node relapse and extranodal relapse cohorts. Genomic characterization of tumor material from diagnosis was performed using panel-based high-throughput sequencing and tumor-infiltrating T cell burden was assessed using immunohistochemistry for CD3+ and CD8+ cells. RESULTS: Isolated lymph node relapse cases demonstrated significantly prolonged postrelapse survival and overall survival vs extranodal relapse upon multivariable analysis (HRmulti = 0.52 [0.33-0.84] and 0.51 [0.31-0.84]). Diagnostic specimens from high-grade serous ovarian carcinomas that subsequently displayed isolated lymph node relapse harbored significantly greater CD3+ and CD8+ cell infiltration compared to extranodal relapse cases (P = .001 and P = .009, Bonferroni-adjusted P = .003 and P = .019). Isolated lymph node relapse high-grade serous ovarian carcinoma cases did not show marked enrichment or depletion of cases with BRCA1/2 mutation or CCNE1 copy number gain when compared to their extranodal relapse counterparts (24.4% vs 19.4% and 18.2% vs 22.6%, P = .865 and P = .900). CONCLUSION: Isolated lymph node relapse ovarian carcinoma represents a distinct clinical entity with favorable outcome compared to extranodal relapse. There was no clear enrichment or depletion of BRCA1/2 mutation or CCNE1 gain in the isolated lymph node relapse ovarian carcinoma cohort compared with extranodal relapse cases, suggesting that these known prognostic genomically defined subtypes of disease do not display markedly altered propensity for isolated lymph node relapse. Diagnostic tumor material from isolated lymph node relapse patients demonstrated greater CD3+ and CD8+ cell infiltration, indicating stronger tumor engagement by T cell populations, which may contribute to the more indolent disease course of isolated lymph node relapse.

Enhanced response rate to pegylated liposomal doxorubicin in high grade serous ovarian carcinomas harbouring BRCA1 and BRCA2 aberrations.

BACKGROUND: Approximately 10-15% of ovarian carcinomas (OC) are attributed to inherited susceptibility, the majority of which are due to mutations in BRCA1 or BRCA2 (BRCA1/2). These patients display superior clinical outcome, including enhanced sensitivity to platinum-based chemotherapy. Here, we seek to investigate whether BRCA1/2 status influences the response rate to single-agent pegylated liposomal doxorubicin (PLD) in high grade serous (HGS) OC. METHODS: One hundred and forty-eight patients treated with single-agent PLD were identified retrospectively from the Edinburgh Ovarian Cancer Database. DNA was extracted from formalin-fixed paraffin-embedded (FFPE) archival tumour material and sequenced using the Ion Ampliseq BRCA1 and BRCA2 panel. A minimum variant allele frequency threshold was applied to correct for sequencing artefacts associated with formalin fixation. RESULTS: A superior response rate to PLD was observed in patients with HGS OC who harboured variants likely to affect BRCA1 or BRCA2 function compared to the BRCA1/2 wild-type population (36%, 9 of 25 patients versus 12.1%, 7 of 58 patients; p = 0.016). An enhanced response rate was also seen in patients harbouring only the BRCA1 SNP rs1799950, predicted to be detrimental to BRCA1 function (50%, 3 of 6 patients versus 12.1%, 7 of 58 patients; p = 0.044). CONCLUSIONS: These data demonstrate that HGS OC patients with BRCA1/2 variants predicted damaging to protein function experience superior sensitivity to PLD, consistent with impaired DNA repair. Further characterisation of rs1799950 is now warranted in relation to chemosensitivity and susceptibility to developing ovarian carcinoma.

Genetic heterogeneity in Cornelia de Lange syndrome (CdLS) and CdLS-like phenotypes with observed and predicted levels of mosaicism.

BACKGROUND: Cornelia de Lange syndrome (CdLS) is a multisystem disorder with distinctive facial appearance, intellectual disability and growth failure as prominent features. Most individuals with typical CdLS have de novo heterozygous loss-of-function mutations in NIPBL with mosaic individuals representing a significant proportion. Mutations in other cohesin components, SMC1A, SMC3, HDAC8 and RAD21 cause less typical CdLS. METHODS: We screened 163 affected individuals for coding region mutations in the known genes, 90 for genomic rearrangements, 19 for deep intronic variants in NIPBL and 5 had whole-exome sequencing. RESULTS: Pathogenic mutations [including mosaic changes] were identified in: NIPBL 46 [3] (28.2%); SMC1A 5 [1] (3.1%); SMC3 5 [1] (3.1%); HDAC8 6 [0] (3.6%) and RAD21 1 [0] (0.6%). One individual had a de novo 1.3 Mb deletion of 1p36.3. Another had a 520 kb duplication of 12q13.13 encompassing ESPL1, encoding separase, an enzyme that cleaves the cohesin ring. Three de novo mutations were identified in ANKRD11 demonstrating a phenotypic overlap with KBG syndrome. To estimate the number of undetected mosaic cases we used recursive partitioning to identify discriminating features in the NIPBL-positive subgroup. Filtering of the mutation-negative group on these features classified at least 18% as 'NIPBL-like'. A computer composition of the average face of this NIPBL-like subgroup was also more typical in appearance than that of all others in the mutation-negative group supporting the existence of undetected mosaic cases. CONCLUSIONS: Future diagnostic testing in 'mutation-negative' CdLS thus merits deeper sequencing of multiple DNA samples derived from different tissues.

Variant detection sensitivity and biases in whole genome and exome sequencing.

BACKGROUND: Less than two percent of the human genome is protein coding, yet that small fraction harbours the majority of known disease causing mutations. Despite rapidly falling whole genome sequencing (WGS) costs, much research and increasingly the clinical use of sequence data is likely to remain focused on the protein coding exome. We set out to quantify and understand how WGS compares with the targeted capture and sequencing of the exome (exome-seq), for the specific purpose of identifying single nucleotide polymorphisms (SNPs) in exome targeted regions. RESULTS: We have compared polymorphism detection sensitivity and systematic biases using a set of tissue samples that have been subject to both deep exome and whole genome sequencing. The scoring of detection sensitivity was based on sequence down sampling and reference to a set of gold-standard SNP calls for each sample. Despite evidence of incremental improvements in exome capture technology over time, whole genome sequencing has greater uniformity of sequence read coverage and reduced biases in the detection of non-reference alleles than exome-seq. Exome-seq achieves 95% SNP detection sensitivity at a mean on-target depth of 40 reads, whereas WGS only requires a mean of 14 reads. Known disease causing mutations are not biased towards easy or hard to sequence areas of the genome for either exome-seq or WGS. CONCLUSIONS: From an economic perspective, WGS is at parity with exome-seq for variant detection in the targeted coding regions. WGS offers benefits in uniformity of read coverage and more balanced allele ratio calls, both of which can in most cases be offset by deeper exome-seq, with the caveat that some exome-seq targets will never achieve sufficient mapped read depth for variant detection due to technical difficulties or probe failures. As WGS is intrinsically richer data that can provide insight into polymorphisms outside coding regions and reveal genomic rearrangements, it is likely to progressively replace exome-seq for many applications.

Quantifying single nucleotide variant detection sensitivity in exome sequencing.

BACKGROUND: The targeted capture and sequencing of genomic regions has rapidly demonstrated its utility in genetic studies. Inherent in this technology is considerable heterogeneity of target coverage and this is expected to systematically impact our sensitivity to detect genuine polymorphisms. To fully interpret the polymorphisms identified in a genetic study it is often essential to both detect polymorphisms and to understand where and with what probability real polymorphisms may have been missed. RESULTS: Using down-sampling of 30 deeply sequenced exomes and a set of gold-standard single nucleotide variant (SNV) genotype calls for each sample, we developed an empirical model relating the read depth at a polymorphic site to the probability of calling the correct genotype at that site. We find that measured sensitivity in SNV detection is substantially worse than that predicted from the naive expectation of sampling from a binomial. This calibrated model allows us to produce single nucleotide resolution SNV sensitivity estimates which can be merged to give summary sensitivity measures for any arbitrary partition of the target sequences (nucleotide, exon, gene, pathway, exome). These metrics are directly comparable between platforms and can be combined between samples to give "power estimates" for an entire study. We estimate a local read depth of 13X is required to detect the alleles and genotype of a heterozygous SNV 95% of the time, but only 3X for a homozygous SNV. At a mean on-target read depth of 20X, commonly used for rare disease exome sequencing studies, we predict 5-15% of heterozygous and 1-4% of homozygous SNVs in the targeted regions will be missed. CONCLUSIONS: Non-reference alleles in the heterozygote state have a high chance of being missed when commonly applied read coverage thresholds are used despite the widely held assumption that there is good polymorphism detection at these coverage levels. Such alleles are likely to be of functional importance in population based studies of rare diseases, somatic mutations in cancer and explaining the "missing heritability" of quantitative traits.

Distinct histopathological features are associated with molecular subtypes and outcome in low grade serous ovarian carcinoma.

Low grade serous ovarian carcinoma (LGSOC) demonstrates unique clinical and molecular features compared to other ovarian cancer types. The relationship between common histological features of LGSOC and molecular events, such as hormone receptor expression patterns and MAPK gene mutation status, remains poorly understood. Recent data suggest some of these molecular features may be biomarkers of response to recently introduced biologically-targeted therapies, namely endocrine therapy and MEK inhibitors. We utilize a cohort of 63 pathologically-confirmed LGSOC cases with whole exome sequencing and hormone receptor expression data to investigate these relationships. LGSOC cases demonstrated uniformly high oestrogen receptor (ER) expression, but variable progesterone receptor (PR) expression intensity. 60% and 37% of cases demonstrated micropapillary and macropapillary patterns of stromal invasion, respectively. 63% of cases demonstrated desmoplasia, which was significantly associated with advanced disease stage and visible residual disease after cytoreductive surgery. MAPK-mutant cases (KRAS, BRAF, NRAS) more frequently demonstrated macropapillary stromal invasion, while Chr1p loss was associated with desmoplasia and low PR expression. Presence of micropapillary stromal invasion and low PR expression were associated with significantly poorer survival after accounting for stage and residual disease status. Together, these data identify novel relationships between histopathological features and molecularly-defined subgroups in LGSOC.

In Vivo Modeling of Patient Genetic Heterogeneity Identifies New Ways to Target Cholangiocarcinoma.

Intrahepatic cholangiocarcinoma (ICC) is an aggressive malignancy of the bile ducts within the liver characterized by high levels of genetic heterogeneity. In the context of such genetic variability, determining which oncogenic mutations drive ICC growth has been difficult, and developing modes of patient stratification and targeted therapies remains challenging. Here we model the interactions between rare mutations with more common driver genes and combine in silico analysis of patient data with highly multiplexed in vivo CRISPR-spCas9 screens to perform a functional in vivo study into the role genetic heterogeneity plays in driving ICC. Novel tumor suppressors were uncovered, which, when lost, cooperate with the RAS oncoprotein to drive ICC growth. Focusing on a set of driver mutations that interact with KRAS to initiate aggressive, sarcomatoid-type ICC revealed that tumor growth relies on Wnt and PI3K signaling. Pharmacologic coinhibition of Wnt and PI3K in vivo impeded ICC growth regardless of mutational profile. Therefore, Wnt and PI3K activity should be considered as a signature by which patients can be stratified for treatment independent of tumor genotype, and inhibitors of these pathways should be levied to treat ICC. SIGNIFICANCE: This work shows that, despite significant genetic heterogeneity, intrahepatic cholangiocarcinoma relies on a limited number of signaling pathways to grow, suggesting common therapeutic vulnerabilities across patients.

Multiomic Characterization of High-Grade Serous Ovarian Carcinoma Enables High-Resolution Patient Stratification.

PURPOSE: High-grade serous ovarian carcinoma (HGSOC) is the most common ovarian cancer type; most patients experience disease recurrence that accumulates chemoresistance, leading to treatment failure. Genomic and transcriptomic features have been associated with differential outcome and treatment response. However, the relationship between events at the gene sequence, copy number, and gene-expression levels remains poorly defined. EXPERIMENTAL DESIGN: We perform multiomic characterization of a large HGSOC cohort (n = 362) with detailed clinical annotation to interrogate the relationship between patient subgroups defined by specific molecular events. RESULTS: BRCA2-mutant (BRCA2m) and EMSY-overexpressing cases demonstrated prolonged survival [multivariable hazard ratios (HR) 0.40 and 0.51] and significantly higher first- and second-line chemotherapy response rate. CCNE1-gained (CCNE1g) cases demonstrated underrepresentation of FIGO stage IV cases, with shorter survival but no significant difference in treatment response. We demonstrate marked overlap between the TCGA- and Tothill-derived subtypes. IMR/C2 cases displayed higher BRCA1/2m frequency (25.5%, 32.5%) and significantly greater immune cell infiltration, whereas PRO/C5 cases had the highest CCNE1g rate (23.9%, 22.2%) and were uniformly low in immune cell infiltration. The survival benefit for cases with aberrations in homologous recombination repair (HRR) genes was apparent across all transcriptomic subtypes (HR range, 0.48-0.68). There was significant co-occurrence of RB loss and HRR gene aberrations; RB loss was further associated with favorable survival within HRR-aberrant cases (multivariable HR, 0.50). CONCLUSIONS: These data paint a high-resolution picture of the molecular landscape in HGSOC, better defining patients who may benefit most from specific molecular therapeutics and highlighting those for whom novel treatment strategies are needed to improve outcomes.

Integrated molecular characterisation of endometrioid ovarian carcinoma identifies opportunities for stratification.

Endometrioid ovarian carcinoma (EnOC) is an under-investigated ovarian cancer type. Recent studies have described disease subtypes defined by genomics and hormone receptor expression patterns; here, we determine the relationship between these subtyping layers to define the molecular landscape of EnOC with high granularity and identify therapeutic vulnerabilities in high-risk cases. Whole exome sequencing data were integrated with progesterone and oestrogen receptor (PR and ER) expression-defined subtypes in 90 EnOC cases following robust pathological assessment, revealing dominant clinical and molecular features in the resulting integrated subtypes. We demonstrate significant correlation between subtyping approaches: PR-high (PR + /ER + , PR + /ER-) cases were predominantly CTNNB1-mutant (73.2% vs 18.4%, P < 0.001), while PR-low (PR-/ER + , PR-/ER-) cases displayed higher TP53 mutation frequency (38.8% vs 7.3%, P = 0.001), greater genomic complexity (P = 0.007) and more frequent copy number alterations (P = 0.001). PR-high EnOC patients experience favourable disease-specific survival independent of clinicopathological and genomic features (HR = 0.16, 95% CI 0.04-0.71). TP53 mutation further delineates the outcome of patients with PR-low tumours (HR = 2.56, 95% CI 1.14-5.75). A simple, routinely applicable, classification algorithm utilising immunohistochemistry for PR and p53 recapitulated these subtypes and their survival profiles. The genomic profile of high-risk EnOC subtypes suggests that inhibitors of the MAPK and PI3K-AKT pathways, alongside PARP inhibitors, represent promising candidate agents for improving patient survival. Patients with PR-low TP53-mutant EnOC have the greatest unmet clinical need, while PR-high tumours-which are typically CTNNB1-mutant and TP53 wild-type-experience excellent survival and may represent candidates for trials investigating de-escalation of adjuvant chemotherapy to agents such as endocrine therapy.

Molecular stratification of endometrioid ovarian carcinoma predicts clinical outcome.

Endometrioid ovarian carcinoma (EnOC) demonstrates substantial clinical and molecular heterogeneity. Here, we report whole exome sequencing of 112 EnOC cases following rigorous pathological assessment. We detect a high frequency of mutation in CTNNB1 (43%), PIK3CA (43%), ARID1A (36%), PTEN (29%), KRAS (26%), TP53 (26%) and SOX8 (19%), a recurrently-mutated gene previously unreported in EnOC. POLE and mismatch repair protein-encoding genes were mutated at lower frequency (6%, 18%) with significant co-occurrence. A molecular taxonomy is constructed, identifying clinically distinct EnOC subtypes: cases with TP53 mutation demonstrate greater genomic complexity, are commonly FIGO stage III/IV at diagnosis (48%), are frequently incompletely debulked (44%) and demonstrate inferior survival; conversely, cases with CTNNB1 mutation, which is mutually exclusive with TP53 mutation, demonstrate low genomic complexity and excellent clinical outcome, and are predominantly stage I/II at diagnosis (89%) and completely resected (87%). Moreover, we identify the WNT, MAPK/RAS and PI3K pathways as good candidate targets for molecular therapeutics in EnOC.

An actionable KCNH2 Long QT Syndrome variant detected by sequence and haplotype analysis in a population research cohort.

The Viking Health Study Shetland is a population-based research cohort of 2,122 volunteer participants with ancestry from the Shetland Isles in northern Scotland. The high kinship and detailed phenotype data support a range of approaches for associating rare genetic variants, enriched in this isolate population, with quantitative traits and diseases. As an exemplar, the c.1750G > A; p.Gly584Ser variant within the coding sequence of the KCNH2 gene implicated in Long QT Syndrome (LQTS), which occurred once in 500 whole genome sequences from this population, was investigated. Targeted sequencing of the KCNH2 gene in family members of the initial participant confirmed the presence of the sequence variant and identified two further members of the same family pedigree who shared the variant. Investigation of these three related participants for whom single nucleotide polymorphism (SNP) array genotypes were available allowed a unique shared haplotype of 1.22 Mb to be defined around this locus. Searching across the full cohort for this haplotype uncovered two additional apparently unrelated individuals with no known genealogical connection to the original kindred. All five participants with the defined haplotype were shown to share the rare variant by targeted Sanger sequencing. If this result were verified in a healthcare setting, it would be considered clinically actionable, and has been actioned in relatives ascertained independently through clinical presentation. The General Practitioners of four study participants with the rare variant were alerted to the research findings by letters outlining the phenotype (prolonged electrocardiographic QTc interval). A lack of detectable haplotype sharing between c.1750G > A; p.Gly584Ser chromosomes from previously reported individuals from Finland and those in this study from Shetland suggests that this mutation has arisen more than once in human history. This study showcases the potential value of isolate population-based research resources for genomic medicine. It also illustrates some challenges around communication of actionable findings in research participants in this context.

CAG-encoded polyglutamine length polymorphism in the human genome.

BACKGROUND: Expansion of polyglutamine-encoding CAG trinucleotide repeats has been identified as the pathogenic mutation in nine different genes associated with neurodegenerative disorders. The majority of individuals clinically diagnosed with spinocerebellar ataxia do not have mutations within known disease genes, and it is likely that additional ataxias or Huntington disease-like disorders will be found to be caused by this common mutational mechanism. We set out to determine the length distributions of CAG-polyglutamine tracts for the entire human genome in a set of healthy individuals in order to characterize the nature of polyglutamine repeat length variation across the human genome, to establish the background against which pathogenic repeat expansions can be detected, and to prioritize candidate genes for repeat expansion disorders. RESULTS: We found that repeats, including those in known disease genes, have unique distributions of glutamine tract lengths, as measured by fragment analysis of PCR-amplified repeat regions. This emphasizes the need to characterize each distribution and avoid making generalizations between loci. The best predictors of known disease genes were occurrence of a long CAG-tract uninterrupted by CAA codons in their reference genome sequence, and high glutamine tract length variance in the normal population. We used these parameters to identify eight priority candidate genes for polyglutamine expansion disorders. Twelve CAG-polyglutamine repeats were invariant and these can likely be excluded as candidates. We outline some confusion in the literature about this type of data, difficulties in comparing such data between publications, and its application to studies of disease prevalence in different populations. Analysis of Gene Ontology-based functions of CAG-polyglutamine-containing genes provided a visual framework for interpretation of these genes' functions. All nine known disease genes were involved in DNA-dependent regulation of transcription or in neurogenesis, as were all of the well-characterized priority candidate genes. CONCLUSION: This publication makes freely available the normal distributions of CAG-polyglutamine repeats in the human genome. Using these background distributions, against which pathogenic expansions can be identified, we have begun screening for mutations in individuals clinically diagnosed with novel forms of spinocerebellar ataxia or Huntington disease-like disorders who do not have identified mutations within the known disease-associated genes.

Novel PEX11B Mutations Extend the Peroxisome Biogenesis Disorder 14B Phenotypic Spectrum and Underscore Congenital Cataract as an Early Feature.

Purpose: Peroxisomes perform complex metabolic and catabolic functions essential for normal growth and development. Mutations in 14 genes cause a spectrum of peroxisomal disease in humans. Most recently, PEX11B was associated with an atypical peroxisome biogenesis disorder (PBD) in a single individual. In this study, we identify further PEX11B cases and delineate associated phenotypes. Methods: Probands from three families underwent next generation sequencing (NGS) for diagnosis of a multisystem developmental disorder. Autozygosity mapping was conducted in one affected sibling pair. ExomeDepth was used to identify copy number variants from NGS data and confirmed by dosage analysis. Biochemical profiling was used to investigate the metabolic signature of the condition. Results: All patients presented with bilateral cataract at birth but the systemic phenotype was variable, including short stature, skeletal abnormalities, and dysmorphism-features not described in the original case. Next generation sequencing identified biallelic loss-of-function mutations in PEX11B as the underlying cause of disease in each case (PEX11B c.235C>T p.(Arg79Ter) homozygous; PEX11B c.136C>T p.(Arg46Ter) homozygous; PEX11B c.595C>T p.(Arg199Ter) heterozygous, PEX11B ex1-3 del heterozygous). Biochemical studies identified very low plasmalogens in one patient, whilst a mildly deranged very long chain fatty acid profile was found in another. Conclusions: Our findings expand the phenotypic spectrum of the condition and underscore congenital cataract as the consistent primary presenting feature. We also find that biochemical measurements of peroxisome function may be disturbed in some cases. Furthermore, diagnosis by NGS is proficient and may circumvent the requirement for an invasive skin biopsy for disease identification from fibroblast cells.

Loss of ALDH18A1 function is associated with a cellular lipid droplet phenotype suggesting a link between autosomal recessive cutis laxa type 3A and Warburg Micro syndrome.

Autosomal recessive cutis laxa type 3A is caused by mutations in ALDH18A1, a gene encoding the mitochondrial enzyme Δ(1)-pyrroline-5-carboxylate synthase (P5CS). It is a rare disorder with only six pathogenic mutations and 10 affected individuals from five families previously described in the literature. Here we report the identification of novel compound heterozygous missense mutations in two affected siblings from a Lebanese family by whole-exome sequencing. The mutations alter a conserved C-terminal domain of the encoded protein and reduce protein stability as determined through Western blot analysis of patient fibroblasts. Patient fibroblasts exhibit a lipid droplet phenotype similar to that recently reported in Warburg Micro syndrome, a disorder with similar features but hitherto unrelated cellular etiology.