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

Genotype-phenotype correlations of TGFBI p.Leu509Pro, p.Leu509Arg, p.Val613Gly, and the allelic association of p.Met502Val-p.Arg555Gln mutations.

International audiencePURPOSE: Investigate the genotype-phenotype correlations for five TGFBI (transforming growth factor, beta-induced) mutations including one novel pathogenic variant and one complex allele affecting the fourth FAS1 domain of keratoepithelin, and their potential effects on the protein's structure. METHODS: Three unrelated families were clinically diagnosed with lattice corneal dystrophy (CD) and one with an unclassified CD of Bowman's layer. Mutations in the TGFBI gene were detected by direct sequencing, and the functional impact of each variant was predicted using in silico algorithms. Corneal phenotypes, including histological examinations, were compared with the literature data. Furthermore, molecular modeling studies of these mutations were performed. RESULTS: Two distinct missense mutations affecting the same residue at position 509 of keratoepithelin: p.Leu509Pro (c.1526T>C) and p.Leu509Arg (c.1526T>G) were found to be associated with a lattice-type CD. The novel p.Val613Gly (c.1828T>G) TGFBI mutation was found in a sporadic case of an Algerian individual affected by lattice CD. Finally, the Bowman's layer CD was linked to the association in cis of the p.Met502Val and p.Arg555Gln variants, leading to the reclassification of this CD as atypical Thiel-Behnke CD. Structural modeling of these TGFBI mutations argues in favor of these mutations being responsible for instability and/or incorrect folding of keratoepithelin, predictions that are compatible with the clinical diagnoses. CONCLUSIONS: Description of a novel TGFBI mutation and a complex TGFBI allele further extends the mutational spectrum of TGFBI. Moreover, we show convincing evidence that TGFBI mutations affecting Leu509 are linked to the lattice phenotype in two unrelated French families, contrasting with findings previously reported. The p.Leu509Pro was reported to be associated with both amyloid and non-amyloid aggregates, whereas p.Leu509Arg has been described as being responsible for Epithelial Basement Membrane Dystrophy (EBMD)

The Immune Signaling Adaptor LAT Contributes to the Neuroanatomical Phenotype of 16p11.2 BP2-BP3 CNVs

Copy-number changes in 16p11.2 contribute significantly to neuropsychiatric traits. Besides the 600 kb BP4-BP5 CNV found in 0.5%–1% of individuals with autism spectrum disorders and schizophrenia and whose rearrangement causes reciprocal defects in head size and body weight, a second distal 220 kb BP2-BP3 CNV is likewise a potent driver of neuropsychiatric, anatomical, and metabolic pathologies. These two CNVs are engaged in complex reciprocal chromatin looping, intimating a functional relationship between genes in these regions that might be relevant to pathomechanism. We assessed the drivers of the distal 16p11.2 duplication by overexpressing each of the nine encompassed genes in zebrafish. Only overexpression of LAT induced a reduction of brain proliferating cells and concomitant microcephaly. Consistently, suppression of the zebrafish ortholog induced an increase of proliferation and macrocephaly. These phenotypes were not unique to zebrafish; Lat knockout mice show brain volumetric changes. Consistent with the hypothesis that LAT dosage is relevant to the CNV pathology, we observed similar effects upon overexpression of CD247 and ZAP70, encoding members of the LAT signalosome. We also evaluated whether LAT was interacting with KCTD13, MVP, and MAPK3, major driver and modifiers of the proximal 16p11.2 600 kb BP4-BP5 syndromes, respectively. Co-injected embryos exhibited an increased microcephaly, suggesting the presence of genetic interaction. Correspondingly, carriers of 1.7 Mb BP1-BP5 rearrangements that encompass both the BP2-BP3 and BP4-BP5 loci showed more severe phenotypes. Taken together, our results suggest that LAT, besides its well-recognized function in T cell development, is a major contributor of the 16p11.2 220 kb BP2-BP3 CNV-associated neurodevelopmental phenotypes

Chromosomal contacts connect loci associated with autism, BMI and head circumference phenotypes

Copy number variants (CNVs) are major contributors to genomic imbalance disorders. Phenotyping of 137 unrelated deletion and reciprocal duplication carriers of the distal 16p11.2 220 kb BP2-BP3 interval showed that these rearrangements are associated with autism spectrum disorders and mirror phenotypes of obesity/underweight and macrocephaly/microcephaly. Such phenotypes were previously associated with rearrangements of the non-overlapping proximal 16p11.2 600 kb BP4-BP5 interval. These two CNV-prone regions at 16p11.2 are reciprocally engaged in complex chromatin looping, as successfully confirmed by 4C-seq, fluorescence in situ hybridization and Hi-C, as well as coordinated expression and regulation of encompassed genes. We observed that genes differentially expressed in 16p11.2 BP4-BP5 CNV carriers are concomitantly modified in their chromatin interactions, suggesting that disruption of chromatin interplays could participate in the observed phenotypes. We also identified cis- and trans-acting chromatin contacts to other genomic regions previously associated with analogous phenotypes. For example, we uncovered that individuals with reciprocal rearrangements of the trans-contacted 2p15 locus similarly display mirror phenotypes on head circumference and weight. Our results indicate that chromosomal contacts’ maps could uncover functionally and clinically related genes