Whole-exome sequencing reveals ZNF408 as a new gene associated with autosomal recessive retinitis pigmentosa with vitreal alterations

Retinitis pigmentosa (RP) is a group of progressive inherited retinal dystrophies that cause visual impairment as a result of photoreceptor cell death. RP is heterogeneous, both clinically and genetically making difficult to establish precise genotype–phenotype correlations. In a Spanish family with autosomal recessive RP (arRP), homozygosity mapping and whole-exome sequencing led to the identification of a homozygous mutation (c.358_359delGT; p.Ala122Leufs*2) in the ZNF408 gene. A screening performed in 217 additional unrelated families revealed another homozygous mutation (c.1621C>T; p.Arg541Cys) in an isolated RP case. ZNF408 encodes a transcription factor that harbors 10 predicted C2H2-type fingers thought to be implicated in DNA binding. To elucidate the ZNF408 role in the retina and the pathogenesis of these mutations we have performed different functional studies. By immunohistochemical analysis in healthy human retina, we identified that ZNF408 is expressed in both cone and rod photoreceptors, in a specific type of amacrine and ganglion cells, and in retinal blood vessels. ZNF408 revealed a cytoplasmic localization and a nuclear distribution in areas corresponding with the euchromatin fraction. Immunolocalization studies showed a partial mislocalization of the p.Arg541Cys mutant protein retaining part of the WT protein in the cytoplasm. Our study demonstrates that ZNF408, previously associated with Familial Exudative Vitreoretinopathy (FEVR), is a new gene causing arRP with vitreous condensations supporting the evidence that this protein plays additional functions into the human retina.

Whole Exome Sequencing widens the genetic landscape of Hereditary Optic Neuropathies

Hereditary optic neuropathies (HON) are a genetic cause of visual impairment characterized by degeneration of retinal ganglion cells. The majority of HON are caused by pathogenic variants in mtDNA genes and in gene OPA1. However, several other genes can cause optic atrophy and can only be identified by high throughput genetic analysis. Whole Exome Sequencing (WES) is becoming the primary choice in rare disease molecular diagnosis, being both cost effective and informative. We performed WES on a cohort of 106 cases, of which 74 isolated ON patients (ON) and 32 syndromic ON patients (sON). The total diagnostic yield amounts to 27%, slightly higher for syndromic ON (31%) than for isolated ON (26%). The majority of genes found are related to mitochondrial function and already reported for harbouring HON pathogenic variants: ACO2, AFG3L2, C19orf12, DNAJC30, FDXR, MECR, MTFMT, NDUFAF2, NDUFB11, NDUFV2, OPA1, PDSS1, SDHA, SSBP1, and WFS1. Among these OPA1, ACO2, and WFS1 were confirmed as the most relevant genetic causes of ON. Moreover, several genes were identified, especially in sON patients, with direct impairment of non-mitochondrial molecular pathways: from autophagy and ubiquitin system (LYST, SNF8, WDR45, UCHL1), to neural cells development and function (KIF1A, GFAP, EPHB2, CACNA1A, CACNA1F), but also vitamin metabolism (SLC52A2, BTD), cilia structure (USH2A), and nuclear pore shuttling (NUTF2). Functional validation on yeast model was performed for pathogenic variants detected in MECR, MTFMT, SDHA, and UCHL1 genes. For SDHA and UCHL1 also muscle biopsy and fibroblast cell lines from patients were analysed, pointing to possible pathogenic mechanisms that will be investigated in further studies. In conclusion, WES proved to be an efficient tool when applied to our ON cohort, for both common disease-genes identification and novel genes discovery. It is therefore recommended to consider WES in ON molecular diagnostic pipeline, as for other rare genetic diseases.

Whole genome sequencing in patients with retinitis pigmentosa reveals pathogenic DNA structural changes and NEK2 as a new disease gene.

We performed whole genome sequencing in 16 unrelated patients with autosomal recessive retinitis pigmentosa (ARRP), a disease characterized by progressive retinal degeneration and caused by mutations in over 50 genes, in search of pathogenic DNA variants. Eight patients were from North America, whereas eight were Japanese, a population for which ARRP seems to have different genetic drivers. Using a specific workflow, we assessed both the coding and noncoding regions of the human genome, including the evaluation of highly polymorphic SNPs, structural and copy number variations, as well as 69 control genomes sequenced by the same procedures. We detected homozygous or compound heterozygous mutations in 7 genes associated with ARRP (USH2A, RDH12, CNGB1, EYS, PDE6B, DFNB31, and CERKL) in eight patients, three Japanese and five Americans. Fourteen of the 16 mutant alleles identified were previously unknown. Among these, there was a 2.3-kb deletion in USH2A and an inverted duplication of ∼446 kb in EYS, which would have likely escaped conventional screening techniques or exome sequencing. Moreover, in another Japanese patient, we identified a homozygous frameshift (p.L206fs), absent in more than 2,500 chromosomes from ethnically matched controls, in the ciliary gene NEK2, encoding a serine/threonine-protein kinase. Inactivation of this gene in zebrafish induced retinal photoreceptor defects that were rescued by human NEK2 mRNA. In addition to identifying a previously undescribed ARRP gene, our study highlights the importance of rare structural DNA variations in Mendelian diseases and advocates the need for screening approaches that transcend the analysis of the coding sequences of the human genome.

New insights into the performance of human whole-exome capture platforms.

Whole exome sequencing (WES) is increasingly used in research and diagnostics. WES users expect coverage of the entire coding region of known genes as well as sufficient read depth for the covered regions. It is, however, unknown which recent WES platform is most suitable to meet these expectations. We present insights into the performance of the most recent standard exome enrichment platforms from Agilent, NimbleGen and Illumina applied to six different DNA samples by two sequencing vendors per platform. Our results suggest that both Agilent and NimbleGen overall perform better than Illumina and that the high enrichment performance of Agilent is stable among samples and between vendors, whereas NimbleGen is only able to achieve vendor- and sample-specific best exome coverage. Moreover, the recent Agilent platform overall captures more coding exons with sufficient read depth than NimbleGen and Illumina. Due to considerable gaps in effective exome coverage, however, the three platforms cannot capture all known coding exons alone or in combination, requiring improvement. Our data emphasize the importance of evaluation of updated platform versions and suggest that enrichment-free whole genome sequencing can overcome the limitations of WES in sufficiently covering coding exons, especially GC-rich regions, and in characterizing structural variants.

Characterizing Genetic Drivers of Lymphoma through High-Throughput Sequencing

The advent of next-generation sequencing, now nearing a decade in age, has enabled, among other capabilities, measurement of genome-wide sequence features at unprecedented scale and resolution.

In this dissertation, I describe work to understand the genetic underpinnings of non-Hodgkin’s lymphoma through exploration of the epigenetics of its cell of origin, initial characterization and interpretation of driver mutations, and finally, a larger-scale, population-level study that incorporates mutation interpretation with clinical outcome.

In the first research chapter, I describe genomic characteristics of lymphomas through the lens of their cells of origin. Just as many other cancers, such as breast cancer or lung cancer, are categorized based on their cell of origin, lymphoma subtypes can be examined through the context of their normal B Cells of origin, Naïve, Germinal Center, and post-Germinal Center. By applying integrative analysis of the epigenetics of normal B Cells of origin through chromatin-immunoprecipitation sequencing, we find that differences in normal B Cell subtypes are reflected in the mutational landscapes of the cancers that arise from them, namely Mantle Cell, Burkitt, and Diffuse Large B-Cell Lymphoma.

In the next research chapter, I describe our first endeavor into understanding the genetic heterogeneity of Diffuse Large B Cell Lymphoma, the most common form of non-Hodgkin’s lymphoma, which affects 100,000 patients in the world. Through whole-genome sequencing of 1 case as well as whole-exome sequencing of 94 cases, we characterize the most recurrent genetic features of DLBCL and lay the groundwork for a larger study.

In the last research chapter, I describe work to characterize and interpret the whole exomes of 1001 cases of DLBCL in the largest single-cancer study to date. This highly-powered study enabled sub-gene, gene-level, and gene-network level understanding of driver mutations within DLBCL. Moreover, matched genomic and clinical data enabled the connection of these driver mutations to clinical features such as treatment response or overall survival. As sequencing costs continue to drop, whole-exome sequencing will become a routine clinical assay, and another diagnostic dimension in addition to existing methods such as histology. However, to unlock the full utility of sequencing data, we must be able to interpret it. This study undertakes a first step in developing the understanding necessary to uncover the genomic signals of DLBCL hidden within its exomes. However, beyond the scope of this one disease, the experimental and analytical methods can be readily applied to other cancer sequencing studies.

Thus, this dissertation leverages next-generation sequencing analysis to understand the genetic underpinnings of lymphoma, both by examining its normal cells of origin as well as through a large-scale study to sensitively identify recurrently mutated genes and their relationship to clinical outcome.

A machine learning based method to detect genomic imbalances exploiting X chromosome exome reads

Whole Exome Sequencing (WES) is rapidly becoming the first-tier test in clinics, both thanks to its declining costs and the development of new platforms that help clinicians in the analysis and interpretation of SNV and InDels. However, we still know very little on how CNV detection could increase WES diagnostic yield. A plethora of exome CNV callers have been published over the years, all showing good performances towards specific CNV classes and sizes, suggesting that the combination of multiple tools is needed to obtain an overall good detection performance. Here we present TrainX, a ML-based method for calling heterozygous CNVs in WES data using EXCAVATOR2 Normalized Read Counts. We select males and females’ non pseudo-autosomal chromosome X alignments to construct our dataset and train our model, make predictions on autosomes target regions and use HMM to call CNVs. We compared TrainX against a set of CNV tools differing for the detection method (GATK4 gCNV, ExomeDepth, DECoN, CNVkit and EXCAVATOR2) and found that our algorithm outperformed them in terms of stability, as we identified both deletions and duplications with good scores (0.87 and 0.82 F1-scores respectively) and for sizes reaching the minimum resolution of 2 target regions. We also evaluated the method robustness using a set of WES and SNP array data (n=251), part of the Italian cohort of Epi25 collaborative, and were able to retrieve all clinical CNVs previously identified by the SNP array. TrainX showed good accuracy in detecting heterozygous CNVs of different sizes, making it a promising tool to use in a diagnostic setting.

Mutations in NOTCH2 cause Hajdu-Cheney syndrome, a disorder of severe and progressive bone loss

We used an exome-sequencing strategy and identified an allelic series of NOTCH2 mutations in Hajdu-Cheney syndrome, an autosomal dominant multisystem disorder characterized by severe and progressive bone loss. The Hajdu-Cheney syndrome mutations are predicted to lead to the premature truncation of NOTCH2 with either disruption or loss of the C-terminal proline-glutamate-serine-threonine-rich proteolytic recognition sequence, the absence of which has previously been shown to increase Notch signaling.

A Novel Autosomal Recessive GJA1 Missense Mutation Linked to Craniometaphyseal Dysplasia

Craniometaphyseal dysplasia (CMD) is a rare sclerosing skeletal disorder with progressive hyperostosis of craniofacial bones. CMD can be inherited in an autosomal dominant (AD) trait or occur after de novo mutations in the pyrophosphate transporter ANKH. Although the autosomal recessive (AR)form of CMD had been mapped to 6q21-22 the mutation has been elusive. In this study, we performed whole-exome sequencing for one subject with AR CMD and identified a novel missense mutation (c.716G>A, p.Arg239Gln) in the C-terminus of the gap junction protein alpha-1 (GJA1) coding for connexin 43 (Cx43). We confirmed this mutation in 6 individuals from 3 additional families. The homozygous mutation cosegregated only with affected family members. Connexin 43 is a major component of gap junctions in osteoblasts, osteocytes, osteoclasts and chondrocytes. Gap junctions are responsible for the diffusion of low molecular weight molecules between cells. Mutations in Cx43 cause several dominant and recessive disorders involving developmental abnormalities of bone such as dominant and recessive oculodentodigital dysplasia (ODDD; MIM #164200, 257850) and isolated syndactyly type III (MIM #186100), the characteristic digital anomaly in ODDD. However, characteristic ocular and dental features of ODDD as well as syndactyly are absent in patients with the recessive Arg239Gln Cx43 mutation. Bone remodeling mechanisms disrupted by this novel Cx43 mutation remain to be elucidated.

Host factors and genetic susceptibility to infections due to intracellular bacteria and fastidious organisms.

While genetic polymorphisms play a paramount role in tuberculosis (TB), less is known about their contribution to the severity of diseases caused by other intracellular bacteria and fastidious microorganisms. We searched electronic databases for observational studies reporting on host factors and genetic predisposition to infections caused by intracellular fastidious bacteria published up to 30 May 2014. The contribution of genetic polymorphisms was documented for TB. This includes genetic defects in the mononuclear phagocyte/T helper cell type 1 (Th1) pathway contributing to disseminated TB disease in children and genome-wide linkage analysis (GWAS) in reactivated pulmonary TB in adults. Similarly, experimental studies supported the role of host genetic factors in the clinical presentation of illnesses resulting from other fastidious intracellular bacteria. These include IL-6 -174G/C or low mannose-binding (MBL) polymorphisms, which are incriminated in chronic pulmonary conditions triggered by C. pneumoniae, type 2-like cytokine secretion polymorphisms, which are correlated with various clinical patterns of M. pneumoniae infections, and genetic variation in the NOD2 gene, which is an indicator of tubal pathology resulting from Chamydia trachomatis infections. Monocyte/macrophage migration and T lymphocyte recruitment defects are corroborated to ineffective granuloma formation observed among patients with chronic Q fever. Similar genetic polymorphisms have also been suggested for infections caused by T. whipplei although not confirmed yet. In conclusion, this review supports the paramount role of genetic factors in clinical presentations and severity of infections caused by intracellular fastidious bacteria. Genetic predisposition should be further explored through such as exome sequencing.

Chondrodysplasia and Abnormal Joint Development Associated with Mutations in IMPAD1, Encoding the Golgi-Resident Nucleotide Phosphatase, gPAPP.

We used whole-exome sequencing to study three individuals with a distinct condition characterized by short stature, chondrodysplasia with brachydactyly, congenital joint dislocations, cleft palate, and facial dysmorphism. Affected individuals carried homozygous missense mutations in IMPAD1, the gene coding for gPAPP, a Golgi-resident nucleotide phosphatase that hydrolyzes phosphoadenosine phosphate (PAP), the byproduct of sulfotransferase reactions, to AMP. The mutations affected residues in or adjacent to the phosphatase active site and are predicted to impair enzyme activity. A fourth unrelated patient was subsequently found to be homozygous for a premature termination codon in IMPAD1. Impad1 inactivation in mice has previously been shown to produce chondrodysplasia with abnormal joint formation and impaired proteoglycan sulfation. The human chondrodysplasia associated with gPAPP deficiency joins a growing number of skeletoarticular conditions associated with defective synthesis of sulfated proteoglycans, highlighting the importance of proteoglycans in the development of skeletal elements and joints.

Mutations in the colony stimulating factor 1 receptor (CSF1R) gene cause hereditary diffuse leukoencephalopathy with spheroids.

Hereditary diffuse leukoencephalopathy with spheroids (HDLS) is an autosomal-dominant central nervous system white-matter disease with variable clinical presentations, including personality and behavioral changes, dementia, depression, parkinsonism, seizures and other phenotypes. We combined genome-wide linkage analysis with exome sequencing and identified 14 different mutations affecting the tyrosine kinase domain of the colony stimulating factor 1 receptor (encoded by CSF1R) in 14 families with HDLS. In one kindred, we confirmed the de novo occurrence of the mutation. Follow-up sequencing identified an additional CSF1R mutation in an individual diagnosed with corticobasal syndrome. In vitro, CSF-1 stimulation resulted in rapid autophosphorylation of selected tyrosine residues in the kinase domain of wild-type but not mutant CSF1R, suggesting that HDLS may result from partial loss of CSF1R function. As CSF1R is a crucial mediator of microglial proliferation and differentiation in the brain, our findings suggest an important role for microglial dysfunction in HDLS pathogenesis.

Mutations at a single codon in Mad homology 2 domain of SMAD4 cause Myhre syndrome.

Myhre syndrome (MIM 139210) is a developmental disorder characterized by short stature, short hands and feet, facial dysmorphism, muscular hypertrophy, deafness and cognitive delay. Using exome sequencing of individuals with Myhre syndrome, we identified SMAD4 as a candidate gene that contributes to this syndrome on the basis of its pivotal role in the bone morphogenetic pathway (BMP) and transforming growth factor (TGF)-β signaling. We identified three distinct heterozygous missense SMAD4 mutations affecting the codon for Ile500 in 11 individuals with Myhre syndrome. All three mutations are located in the region of SMAD4 encoding the Mad homology 2 (MH2) domain near the site of monoubiquitination at Lys519, and we found a defect in SMAD4 ubiquitination in fibroblasts from affected individuals. We also observed decreased expression of downstream TGF-β target genes, supporting the idea of impaired TGF-β-mediated transcriptional control in individuals with Myhre syndrome.

Prostaglandin transporter mutations cause pachydermoperiostosis with myelofibrosis.

Pachydermoperiostosis, or primary hypertrophic osteoarthropathy (PHO), is an inherited multisystem disorder, whose features closely mimic the reactive osteoarthropathy that commonly accompanies neoplastic and inflammatory pathologies. We previously described deficiency of the prostaglandin-degrading enzyme 15-hydroxyprostaglandin dehydrogenase (HPGD) as a cause of this condition, implicating elevated circulating prostaglandin E(2) (PGE(2) ) as causative of PHO, and perhaps also as the principal mediator of secondary HO. However, PHO is genetically heterogeneous. Here, we use whole-exome sequencing to identify recessive mutations of the prostaglandin transporter SLCO2A1, in individuals lacking HPGD mutations. We performed exome sequencing of four probands with severe PHO, followed by conventional mutation analysis of SLCO2A1 in nine others. Biallelic SLCO2A1 mutations were identified in 12 of the 13 families. Affected individuals had elevated urinary PGE(2) , but unlike HPGD-deficient patients, also excreted considerable quantities of the PGE(2) metabolite, PGE-M. Clinical differences between the two groups were also identified, notably that SLCO2A1-deficient individuals have a high frequency of severe anemia due to myelofibrosis. These findings reinforce the key role of systemic or local prostaglandin excess as the stimulus to HO. They also suggest that the induction or maintenance of hematopoietic stem cells by prostaglandin may depend upon transporter activity. Hum Mutat 33:1175-1181, 2012. © 2012 Wiley Periodicals, Inc.