Morad Ansari

Variant detection sensitivity and biases in whole genome and exome sequencing

BackgroundLess 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.ResultsWe 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.ConclusionsFrom 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.

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.

Loss-of-function HDAC8 mutations cause a phenotypic spectrum of Cornelia de Lange syndrome-like features, ocular hypertelorism, large fontanelle and X-linked inheritance

Cornelia de Lange syndrome (CdLS) is a multisystem genetic disorder with distinct facies, growth failure, intellectual disability, distal limb anomalies, gastrointestinal and neurological disease. Mutations in NIPBL, encoding a cohesin regulatory protein, account for >80% of cases with typical facies. Mutations in the core cohesin complex proteins, encoded by the SMC1A, SMC3 and RAD21 genes, together account for ∼5% of subjects, often with atypical CdLS features. Recently, we identified mutations in the X-linked gene HDAC8 as the cause of a small number of CdLS cases. Here, we report a cohort of 38 individuals with an emerging spectrum of features caused by HDAC8 mutations. For several individuals, the diagnosis of CdLS was not considered prior to genomic testing. Most mutations identified are missense and de novo. Many cases are heterozygous females, each with marked skewing of X-inactivation in peripheral blood DNA. We also identified eight hemizygous males who are more severely affected. The craniofacial appearance caused by HDAC8 mutations overlaps that of typical CdLS but often displays delayed anterior fontanelle closure, ocular hypertelorism, hooding of the eyelids, a broader nose and dental anomalies, which may be useful discriminating features. HDAC8 encodes the lysine deacetylase for the cohesin subunit SMC3 and analysis of the functional consequences of the missense mutations indicates that all cause a loss of enzymatic function. These data demonstrate that loss-of-function mutations in HDAC8 cause a range of overlapping human developmental phenotypes, including a phenotypically distinct subgroup of CdLS.

Clinical and mutation analysis of 51 probands with anophthalmia and/or severe microphthalmia from a single center.

Clinical evaluation and mutation analysis was performed in 51 consecutive probands with severe eye malformations – anophthalmia and/or severe microphthalmia – seen in a single specialist ophthalmology center. The mutation analysis consisted of bidirectional sequencing of the coding regions of SOX2, OTX2, PAX6 (paired domain), STRA6, BMP4, SMOC1, FOXE3, and RAX, and genome‐wide array‐based copy number assessment. Fifteen (29.4%) of the 51 probands had likely causative mutations affecting SOX2 (9/51), OTX2 (5/51), and STRA6 (1/51). Of the cases with bilateral anophthalmia, 9/12 (75%) were found to be mutation positive. Three of these mutations were large genomic deletions encompassing SOX2 (one case) or OTX2 (two cases). Familial inheritance of three intragenic, plausibly pathogenic, and heterozygous mutations was observed. An unaffected carrier parent of an affected child with an identified OTX2 mutation confirmed the previously reported nonpenetrance for this disorder. Two families with SOX2 mutations demonstrated a parent and child both with significant but highly variable eye malformations. Heterozygous loss‐of‐function mutations in SOX2 and OTX2 are the most common genetic pathology associated with severe eye malformations and bi‐allelic loss‐of‐function in STRA6 is confirmed as an emerging cause of nonsyndromal eye malformations.

Discovery of four recessive developmental disorders using probabilistic genotype and phenotype matching among 4,125 families.

Discovery of most autosomal recessive disease-associated genes has involved analysis of large, often consanguineous multiplex families or small cohorts of unrelated individuals with a well-defined clinical condition. Discovery of new dominant causes of rare, genetically heterogeneous developmental disorders has been revolutionized by exome analysis of large cohorts of phenotypically diverse parent-offspring trios. Here we analyzed 4,125 families with diverse, rare and genetically heterogeneous developmental disorders and identified four new autosomal recessive disorders. These four disorders were identified by integrating Mendelian filtering (selecting probands with rare, biallelic and putatively damaging variants in the same gene) with statistical assessments of (i) the likelihood of sampling the observed genotypes from the general population and (ii) the phenotypic similarity of patients with recessive variants in the same candidate gene. This new paradigm promises to catalyze the discovery of novel recessive disorders, especially those with less consistent or nonspecific clinical presentations and those caused predominantly by compound heterozygous genotypes.

Genetic influences on plasma CFH and CFHR1 concentrations and their role in susceptibility to age-related macular degeneration

It is a longstanding puzzle why non-coding variants in the complement factor H (CFH) gene are more strongly associated with age-related macular degeneration (AMD) than functional coding variants that directly influence the alternative complement pathway. The situation is complicated by tight genetic associations across the region, including the adjacent CFH-related genes CFHR3 and CFHR1, which may themselves influence the alternative complement pathway and are contained within a common deletion (CNP147) which is associated with protection against AMD. It is unclear whether this association is mediated through a protective effect of low plasma CFHR1 concentrations, high plasma CFH or both. We examined the triangular relationships of CFH/CFHR3/CFHR1 genotype, plasma CFH or CFHR1 concentrations and AMD susceptibility in combined case–control (1256 cases, 1020 controls) and cross-sectional population (n = 1004) studies and carried out genome-wide association studies of plasma CFH and CFHR1 concentrations. A non-coding CFH SNP (rs6677604) and the CNP147 deletion were strongly correlated both with each other and with plasma CFH and CFHR1 concentrations. The plasma CFH-raising rs6677604 allele and raised plasma CFH concentration were each associated with AMD protection. In contrast, the protective association of the CNP147 deletion with AMD was not mediated by low plasma CFHR1, since AMD-free controls showed increased plasma CFHR1 compared with cases, but it may be mediated by the association of CNP147 with raised plasma CFH concentration. The results are most consistent with a regulatory locus within a 32 kb region of the CFH gene, with a major effect on plasma CFH concentration and AMD susceptibility.

Concordant association of insulin degrading enzyme gene (IDE) variants with IDE mRNA, Abeta, and Alzheimer's disease.

Background The insulin-degrading enzyme gene (IDE) is a strong functional and positional candidate for late onset Alzheimers disease (LOAD). Methodology/Principal Findings We examined conserved regions of IDE and its 10 kb flanks in 269 AD cases and 252 controls thereby identifying 17 putative functional polymorphisms. These variants formed eleven haplotypes that were tagged with ten variants. Four of these showed significant association with IDE transcript levels in samples from 194 LOAD cerebella. The strongest, rs6583817, which has not previously been reported, showed unequivocal association (p = 1.5×10−8, fold-increase = 2.12,); the eleven haplotypes were also significantly associated with transcript levels (global p = 0.003). Using an in vitro dual luciferase reporter assay, we found that rs6583817 increases reporter gene expression in Be(2)-C (p = 0.006) and HepG2 (p = 0.02) cell lines. Furthermore, using data from a recent genome-wide association study of two Croatian isolated populations (n = 1,879), we identified a proxy for rs6583817 that associated significantly with decreased plasma Aβ40 levels (ß = −0.124, p = 0.011) and total measured plasma Aβ levels (b = −0.130, p = 0.009). Finally, rs6583817 was associated with decreased risk of LOAD in 3,891 AD cases and 3,605 controls. (OR = 0.87, p = 0.03), and the eleven IDE haplotypes (global p = 0.02) also showed significant association. Conclusions Thus, a previously unreported variant unequivocally associated with increased IDE expression was also associated with reduced plasma Aß40 and decreased LOAD susceptibility. Genetic association between LOAD and IDE has been difficult to replicate. Our findings suggest that targeted testing of expression SNPs (eSNPs) strongly associated with altered transcript levels in autopsy brain samples may be a powerful way to identify genetic associations with LOAD that would otherwise be difficult to detect.

Identification of Novel Craniofacial Regulatory Domains Located far Upstream of SOX9 and Disrupted in Pierre Robin Sequence

Mutations in the coding sequence of SOX9 cause campomelic dysplasia (CD), a disorder of skeletal development associated with 46,XY disorders of sex development (DSDs). Translocations, deletions, and duplications within a ∼2 Mb region upstream of SOX9 can recapitulate the CD–DSD phenotype fully or partially, suggesting the existence of an unusually large cis‐regulatory control region. Pierre Robin sequence (PRS) is a craniofacial disorder that is frequently an endophenotype of CD and a locus for isolated PRS at ∼1.2–1.5 Mb upstream of SOX9 has been previously reported. The craniofacial regulatory potential within this locus, and within the greater genomic domain surrounding SOX9, remains poorly defined. We report two novel deletions upstream of SOX9 in families with PRS, allowing refinement of the regions harboring candidate craniofacial regulatory elements. In parallel, ChIP‐Seq for p300 binding sites in mouse craniofacial tissue led to the identification of several novel craniofacial enhancers at the SOX9 locus, which were validated in transgenic reporter mice and zebrafish. Notably, some of the functionally validated elements fall within the PRS deletions. These studies suggest that multiple noncoding elements contribute to the craniofacial regulation of SOX9 expression, and that their disruption results in PRS.

Disruption of SATB2 or its long-range cis-regulation by SOX9 causes a syndromic form of Pierre Robin Sequence

Heterozygous loss-of-function (LOF) mutations in the gene encoding the DNA-binding protein, SATB2, result in micrognathia and cleft palate in both humans and mice. In three unrelated individuals, we show that translocation breakpoints (BPs) up to 896 kb 3′ of SATB2 polyadenylation site cause a phenotype which is indistinguishable from that caused by SATB2 LOF mutations. This syndrome comprises long nose, small mouth, micrognathia, cleft palate, arachnodactyly and intellectual disability. These BPs map to a gene desert between PLCL1 and SATB2. We identified three putative cis-regulatory elements (CRE1–3) using a comparative genomic approach each of which would be placed in trans relative to SATB2 by all three BPs. CRE1–3 each bind p300 and mono-methylated H3K4 consistent with enhancer function. In silico analysis suggested that CRE1–3 contain one or more conserved SOX9-binding sites, and this binding was confirmed using chromatin immunoprecipitation on cells derived from mouse embryonic pharyngeal arch. Interphase bacterial artificial chromosome fluorescence in situ hybridization measurements in embryonic craniofacial tissues showed that the orthologous region in mice exhibits Satb2 expression-dependent chromatin decondensation consistent with Satb2 being a target gene of CRE1–3. To assess their in vivo function, we made multiple stable reporter transgenic lines for each enhancer in zebrafish. CRE2 was shown to drive SATB2-like expression in the embryonic craniofacial region. This expression could be eliminated by mutating the SOX9-binding site of CRE2. These observations suggest that SATB2 and SOX9 may be acting together via complex cis-regulation to coordinate the growth of the developing jaw.

Miller (Genée–Wiedemann) syndrome represents a clinically and biochemically distinct subgroup of postaxial acrofacial dysostosis associated with partial deficiency of DHODH

Biallelic mutations in the gene encoding DHOdehase [dihydroorotate dehydrogenase (DHODH)], an enzyme required for de novo pyrimidine biosynthesis, have been identified as the cause of Miller (Genée-Weidemann or postaxial acrofacial dysostosis) syndrome (MIM 263750). We report compound heterozygous DHODH mutations in four additional families with typical Miller syndrome. Complementation in auxotrophic yeast demonstrated reduced pyrimidine synthesis and in vitro enzymatic analysis confirmed reduced DHOdehase activity in 11 disease-associated missense mutations, with 7 alleles showing discrepant activity between the assays. These discrepancies are partly explained by the domain structure of DHODH and suggest both assays are useful for interpretation of individual alleles. However, in all affected individuals, the genotype predicts that there should be significant residual DHOdehase activity. Urine samples obtained from two mutation-positive cases showed elevated levels of orotic acid (OA) but not dihydroorotate (DHO), an unexpected finding since these represent the product and the substrate of DHODH enzymatic activity, respectively. Screening of four unrelated cases with overlapping but atypical clinical features showed no mutations in either DHODH or the other de novo pyrimidine biosynthesis genes (CAD, UMPS), with these cases also showing normal levels of urinary OA and DHO. In situ analysis of mouse embryos showed Dhodh, Cad and Umps to be strongly expressed in the pharyngeal arch and limb bud, supporting a site- and stage-specific requirement for de novo pyrimidine synthesis. The developmental sensitivity to reduced pyrimidine synthesis capacity may reflect the requirement for an exceptional mitogenic response to growth factor signalling in the affected tissues.