Implementation of non-invasive prenatal testing by semiconductor sequencing in a genetic laboratory

Objectives: To implement non-invasive prenatal testing (NIPT) for fetal aneuploidies with semiconductor sequencing in an academic cytogenomic laboratory and to evaluate the first 15-month experience on clinical samples. Methods: We validated a NIPT protocol for cell-free fetal DNA sequencing from maternal plasma for the detection of trisomy 13, 18 and 21 on a semiconductor sequencing instrument. Fetal DNA fraction calculation for all samples and several quality parameters were implemented in the workflow. One thousand eighty-one clinical NIPT samples were analysed, following the described protocol. Results: Non-invasive prenatal testing was successfully implemented and validated on 201 normal and 74 aneuploid samples. From 1081 clinical samples, 17 samples showed an abnormal result: 14 trisomy 21 samples, one trisomy 18 and one trisomy 16 were detected. Also a maternal copy number variation on chromosome 13 was observed, which could potentially lead to a false positive trisomy 13 result. One sex discordant result was reported, possibly attributable to a vanishing twin. Moreover, our combined fetal fraction calculation enabled a more reliable risk estimate for trisomy 13, 18 and 21. Conclusions: Non-invasive prenatal testing for trisomy 21, 18 and 13 has a very high specificity and sensitivity. Because of several biological phenomena, diagnostic invasive confirmation of abnormal results remains required

Microhomology-mediated mechanisms underlie non-recurrent disease-causing microdeletions of the FOXL2 gene or its regulatory domain

Genomic disorders are often caused by recurrent copy number variations (CNVs), with nonallelic homologous recombination (NAHR) as the underlying mechanism. Recently, several microhomology-mediated repair mechanisms-such as microhomology-mediated end-joining (MMEJ), fork stalling and template switching (FoSTeS), microhomology-mediated break-induced replication (MMBIR), serial replication slippage (SRS), and break-induced SRS (BISRS)-were described in the etiology of non-recurrent CNVs in human disease. In addition, their formation may be stimulated by genomic architectural features. It is, however, largely unexplored to what extent these mechanisms contribute to rare, locus-specific pathogenic CNVs. Here, fine-mapping of 42 microdeletions of the FOXL2 locus, encompassing FOXL2 (32) or its regulatory domain (10), serves as a model for rare, locus-specific CNVs implicated in genetic disease. These deletions lead to blepharophimosis syndrome (BPES), a developmental condition affecting the eyelids and the ovary. For breakpoint mapping we used targeted array-based comparative genomic hybridization (aCGH), quantitative PCR (qPCR), long-range PCR, and Sanger sequencing of the junction products. Microhomology, ranging from 1 bp to 66 bp, was found in 91.7% of 24 characterized breakpoint junctions, being significantly enriched in comparison with a random control sample. Our results show that microhomology-mediated repair mechanisms underlie at least 50% of these microdeletions. Moreover, genomic architectural features, like sequence motifs, non-B DNA conformations, and repetitive elements, were found in all breakpoint regions. In conclusion, the majority of these microdeletions result from microhomology-mediated mechanisms like MMEJ, FoSTeS, MMBIR, SRS, or BISRS. Moreover, we hypothesize that the genomic architecture might drive their formation by increasing the susceptibility for DNA breakage or promote replication fork stalling. Finally, our locus-centered study, elucidating the etiology of a large set of rare microdeletions involved in a monogenic disorder, can serve as a model for other clustered, non-recurrent microdeletions in genetic disease

arrEYE : a customized platform for high-resolution copy number analysis of coding and noncoding regions of known and candidate retinal dystrophy genes and retinal noncoding RNAs

Purpose: Our goal was to design a customized microarray, arrEYE, for high-resolution copy number variant (CNV) analysis of known and candidate genes for inherited retinal dystrophy (iRD) and retina expressed noncoding RNAs (ncRNAs). Methods: arrEYE contains probes for the full genomic region of 106 known iRD genes, including those implicated in retinitis pigmentosa (RP) (the most frequent iRD), cone rod dystrophies, macular dystrophies, and an additional 60 candidate iRD genes and 196 ncRNAs. Eight CNVs in iRD genes identified by other techniques were used as positive controls. The test cohort consisted of 57 patients with autosomal dominant, X-linked, or simplex RP. Results: In an RP patient, a novel heterozygous deletion of exons 7 and 8 of the HGSNAT gene was identified: c.634-408_820+338delins AGAATATG, p.(G1u2 I 2Glyfs*2). A known variant was found on the second allele: c.1843G>A, p.(A1a615Thr). Furthermore, we expanded the allelic spectrum of USH2A and RCBTB1 with novel CNVs. Conclusion: The arrEYE platform revealed subtle single-exon to larger CNVs in iRD genes that could be characterized at the nucleotide level, facilitated by the high resolution of the platform. We report the first CNV in HGSNAT that, combined with another mutation, leads to RP, further supporting its recently identified role in nonsyndromic iRD

Leigh syndrome followed by parkinsonism in an adult with homozygous c.626C > T mutation in MTFMT

Objective : To report the clinical, radiologic, biochemical, and molecular characteristics in a 46-year-old participant with adult-onset Leigh syndrome (LS), followed by parkinsonism. Methods : Case description with diagnostic workup included blood and CSF analysis, skeletal muscle investigations, blue native polyacrylamide gel electrophoresis, whole exome sequencing targeting nuclear genes involved in mitochondrial transcription and translation, cerebral MRI, 123I-FP-CIT brain single-photon emission computed tomography (SPECT), and C-11 raclopride positron emission tomography (PET). Results : The participant was found to have a defect in the oxidative phosphorylation caused by a c.626C>T mutation in the gene coding for mitochondrial methionyl-tRNA formyltransferase (MTFMT), which is a pathogenic mutation affecting intramitochondrial protein translation. The proband had a normal concentration of lactate in blood and no abnormal microscopic findings in skeletal muscle. Cerebral MRI showed bilateral lesions in the striatum, mesencephalon, pons, and medial thalamus. Lactate concentration in CSF was increased. FP-CIT SPECT and C-11 raclopride PET demonstrated a defect in the dopaminergic system. Conclusions : We report on a case with adult-onset LS related to a MTFMT mutation. Two years after the onset of symptoms of LS, the proband developed a parkinson-like disease. The c.626C>T mutation is the most common pathogenic mutation found in 22 patients reported earlier in the literature with a defect in MTFMT. The age of the previously reported cases varied between 14 months and 24 years. Our report expands the phenotypical spectrum of MTFMT-related neurologic disease and provides clinical evidence for involvement of MTFMT in extrapyramidal syndromes

Whole genome amplification with SurePlex results in better copy number alteration detection using sequencing data compared to the MALBAC method

Current whole genome amplification (WGA) methods lead to amplification bias resulting in over-and under-represented regions in the genome. Nevertheless, certain WGA methods, such as SurePlex and subsequent arrayCGH analysis, make it possible to detect copy number alterations (CNAs) at a 10 Mb resolution. A more uniform WGA combined with massive parallel sequencing (MPS), however, could allow detection at higher resolution and lower cost. Recently, MALBAC, a new WGA method, claims unparalleled performance. Here, we compared the well-established SurePlex and MALBAC WGA for their ability to detect CNAs in MPS generated data and, in addition, compared PCR-free MPS library preparation with the standard enrichment PCR library preparation. Results showed that SurePlex amplification led to more uniformity across the genome, allowing for a better CNA detection with less false positives compared to MALBAC amplified samples. An even more uniform coverage was observed in samples following a PCR-free library preparation. In general, the combination of SurePlex and MPS led to the same chromosomal profile compared to a reference arrayCGH from unamplified genomic DNA, underlining the large potential of MPS techniques in CNA detection from a limited number of DNA material

Focal DNA copy number changes in neuroblastoma target MYCN regulated genes

Neuroblastoma is an embryonic tumor arising from immature sympathetic nervous system cells. Recurrent genomic alterations include MYCN and ALK amplification as well as recurrent patterns of gains and losses of whole or large partial chromosome segments. A recent whole genome sequencing effort yielded no frequently recurring mutations in genes other than those affecting ALK. However, the study further stresses the importance of DNA copy number alterations in this disease, in particular for genes implicated in neuritogenesis. Here we provide additional evidence for the importance of focal DNA copy number gains and losses, which are predominantly observed in MYCN amplified tumors. A focal 5 kb gain encompassing the MYCN regulated miR-17,92 cluster as sole gene was detected in a neuroblastoma cell line and further analyses of the array CGH data set demonstrated enrichment for other MYCN target genes in focal gains and amplifications. Next we applied an integrated genomics analysis to prioritize MYCN down regulated genes mediated by MYCN driven miRNAs within regions of focal heterozygous or homozygous deletion. We identified RGS5, a negative regulator of G-protein signaling implicated in vascular normalization, invasion and metastasis, targeted by a focal homozygous deletion, as a new MYCN target gene, down regulated through MYCN activated miRNAs. In addition, we expand the miR-17,92 regulatory network controlling TGFß signaling in neuroblastoma with the ring finger protein 11 encoding gene RNF11, which was previously shown to be targeted by the miR-17,92 member miR-19b. Taken together, our data indicate that focal DNA copy number imbalances in neuroblastoma (1) target genes that are implicated in MYCN signaling, possibly selected to reinforce MYCN oncogene addiction and (2) serve as a resource for identifying new molecular targets for treatment

Finishing the euchromatic sequence of the human genome

The sequence of the human genome encodes the genetic instructions for human physiology, as well as rich information about human evolution. In 2001, the International Human Genome Sequencing Consortium reported a draft sequence of the euchromatic portion of the human genome. Since then, the international collaboration has worked to convert this draft into a genome sequence with high accuracy and nearly complete coverage. Here, we report the result of this finishing process. The current genome sequence (Build 35) contains 2.85 billion nucleotides interrupted by only 341 gaps. It covers ∼99% of the euchromatic genome and is accurate to an error rate of ∼1 event per 100,000 bases. Many of the remaining euchromatic gaps are associated with segmental duplications and will require focused work with new methods. The near-complete sequence, the first for a vertebrate, greatly improves the precision of biological analyses of the human genome including studies of gene number, birth and death. Notably, the human enome seems to encode only 20,000-25,000 protein-coding genes. The genome sequence reported here should serve as a firm foundation for biomedical research in the decades ahead

Development of a sequence analysis pipeline for variant detection and interpretation

Developing a sequencing analysis pipeline for variant discovery and interpretation and the necessary supporting infrastructure is a great challenge. The objectives presented in Chapter 2 were essential for handling the growing volume and velocity at which sequencing data is being generated these last few years. There is a great demand for better platforms to drive our analyses, and as introduced in section 2.1 and discussed at the end of the thesis, a solid but elastic infrastructure is essential to support such platform. The early efforts focused on moving away from custom-built physical application servers to a more flexible private cloud platform that can offer compute power as a service to researchers, eliminating the need to handle lower level system management duties. At the same time we implemented a distributed network filesystem to handle the volume and performance expected for genome analysis, and keep the data replicated to protect against loss. Building on that infrastructure, the following work focused on streamlining an analysis pipeline for structural variation. The introduction of new technologies such as mate-pair/paired- end sequencing and coverage based analysis, spurred many different tools but no clear way to tackle the complete process, starting from reads until we can produce a genetic report. The web based platform, called ViVar, handles intelligent processing of the data, provides careful filtering and a diverse set of visualizations, to distill sequence reads to a set of clinically rel- evant structural variants. Additionally, a SNP analysis focused platform was build. It runs a pipeline conforming to the GATK best practice guide and helps users process the results aided by extensive filtering options to interpret and select the most relevant variants. As should be clear after reading the discussion, several focus areas remain to lift the per- formance and flexibility of the analysis tools and the researchers to the next level. The most important aspects are: building scalable, cost-effective dependable storage systems; optimizing cycle times with better infrastructure and efficient tools improving developer and user produc- tivity; and respecting information security when working with sensitive data. These will be essential drivers of research innovation and novel clinical applications based on sequencing data