NanoImprint: A DNA methylation tool for clinical interpretation and diagnosis of common imprinting disorders using nanopore long-read sequencing.

INTRODUCTION: Long-read whole genome sequencing like Oxford Nanopore Technology, is increasingly being introduced in clinical settings. With its ability to simultaneously call sequence variation and DNA modifications including 5-methylcytosine, nanopore is a promising technology to improve diagnostics of imprinting disorders. METHODS: Currently, no tools to analyze DNA methylation patterns at known clinically relevant imprinted regions are available. Here we present NanoImprint, which generates an easily interpretable report, based on long-read nanopore sequencing, to use for identifying clinical relevant abnormalities in methylation levels at 14 imprinted regions and diagnosis of common imprinting disorders. RESULTS AND CONCLUSION: NanoImprint outputs a summarizing table and visualization plots displays methylation frequency (%) and chromosomal positions for all regions, with phased data color-coded for the two alleles. We demonstrate the utility of NanoImprint using three imprinting disorder samples from patients with Beckwith-Wiedemann syndrome (BWS), Angelman syndrome (AS) and Prader-Willi syndrome (PWS). NanoImprint script is available from https://github.com/carolinehey/NanoImprint.

Comprehensive prenatal diagnostics: Exome versus genome sequencing.

OBJECTIVE: This study aimed to assess the diagnostic yield of prenatal genetic testing using trio whole exome sequencing (WES) and trio whole genome sequencing (WGS) in pregnancies with fetal anomalies by comparing the results with conventional chromosomal microarray (CMA) analysis. METHODS: A total of 40 pregnancies with fetal anomalies or increased nuchal translucency (NT ≥ 5 mm) were included between the 12th and 21st week of gestation. Trio WES/WGS and CMA were performed in all cases. RESULTS: The trio WES/WGS analysis increased the diagnostic yield by 25% in cases with negative CMA results. Furthermore, all six chromosomal aberrations identified by CMA were independently detected by WES/WGS analysis. In total, 16 out of 40 cases obtained a genetic sequence variant, copy number variant, or aneuploidy explaining the phenotype, resulting in an overall WES/WGS diagnostic yield of 40%. WES analysis provided a more reliable identification of mosaic sequence variants than WGS because of its higher sequencing depth. CONCLUSIONS: Prenatal WES/WGS proved to be powerful diagnostic tools for fetal anomalies, surpassing the diagnostic yield of CMA. They have the potential to serve as standalone methods for prenatal diagnosis. The study highlighted the limitations of WGS in accurately detecting mosaic variants, which is particularly relevant when analyzing chorionic villus samples.

Ryanodine receptor 1 related myasthenia like myopathy responsive to pyridostigmine.

Disease causing variants in the Ryanodine receptor 1 (RYR1) gene are a common cause for congenital myopathy and for malignant hyperthermia susceptibility. We report a 17 year old boy with congenital muscle weakness progressing to a myasthenia like myopathy with muscle weakness, fatigability, ptosis, and ophthalmoplegia. Muscle biopsy showed predominance and atrophy of type 1 fibers. Whole-exome trio sequencing revealed three variants in the RYR1-gene in the patient: c.6721C > T,p.(Arg2241*) and c.2122G > A,p.(Asp708Asn) in cis position, and the c.325C > T,p.(Arg109Trp) variant in trans. Treatment with pyridostigmine improved symptoms. This case supports that a myasthenia like phenotype is part of the phenotypic spectrum of RYR1 related disorders, and that treatment with pyridostigmine can be beneficial for patients with this phenotype.

Whole exome sequencing of 28 families of Danish descent reveals novel candidate genes and pathways in developmental dysplasia of the hip.

Developmental dysplasia of the hip (DDH) is a common condition involving instability of the hip with multifactorial etiology. Early diagnosis and treatment are critical as undetected DDH is an important cause of long-term hip complications. Better diagnostics may be achieved through genetic methods, especially for patients with positive family history. Several candidate genes have been reported but the exact molecular etiology of the disease is yet unknown. In the present study, we performed whole exome sequencing of DDH patients from 28 families with at least two affected first-degree relatives. Four genes previously not associated with DDH (METTL21B, DIS3L2, PPP6R2, and TM4SF19) were identified with the same variants shared among affected family members, in more than two families. Among known association genes, we found damaging variants in DACH1, MYH10, NOTCH2, TBX4, EVC2, OTOG, and SHC3. Mutational burden analysis across the families identified 322 candidate genes, and enriched pathways include the extracellular matrix, cytoskeleton, ion-binding, and detection of mechanical stimulus. Taken altogether, our data suggest a polygenic mode of inheritance for DDH, and we propose that an impaired transduction of the mechanical stimulus is involved in the etiopathological mechanism. Our findings refine our current understanding of candidate causal genes in DDH, and provide a foundation for downstream functional studies.

Whole genome sequencing identifies rare genetic variants in familial pancreatic cancer patients.

Pancreatic ductal adenocarcinoma (PDAC) represents one of the most lethal malignancies with very high mortality and short survival time. About 5-10% of the PDAC patients have a familial predisposition to the disease designated as familial pancreatic cancer (FPC), suggesting genetic modulation of FPC pathogenesis. It is estimated that currently identified sequence variants account for less than 20% of the genetic basis of FPC leaving the majority of the genetic architecture unclarified. We performed whole genome sequencing (WGS) analysis on benign formalin-fixed paraffin-embedded (FFPE) tissues from 35 FPC patients focusing on genes enriched by rare and functional sequence variants. We identified 40 genes hosting at least 2 protein truncating variants (PTVs). Significant overlaps of the 40 genes were found (p < 1 × 10-22 ) with cancer genes, cancer driver genes and genes found in previous studies on cancer, including ATM, POLE, BRCA2, TYR03, PABPC1 and SSC5D. The PTV genes are significantly overrepresented in biological pathways in cancer development and progression including extracellular matrix organization, signaling by RHO GTPases and RHO GTPase cycle. Association analysis using external controls detected 6 genes with p < 0.05. The WGS analysis revealed high heterogeneity in the detected rare variants among FPC patients and provides novel genes harboring potential mutational hotspots for future validation and replication.

Prenatal cases with rare RIT1 variants causing severe fetal hydrops and death.

We describe two clinical prenatal cases with rare de novo RIT1 variants, which showed more severe clinical manifestations than other Noonan Syndrome genotypes, resulting in fetal death. Extra attention is recommended when these variants are detected.

Total number of reads affects the accuracy of fetal fraction estimates in NIPT.

BACKGROUND: Sufficient fetal fraction (FF) is crucial for quality control of NIPT (Non-Invasive Prenatal Test) results. Different factors influencing bioinformatic estimation of FF should be considered when implementing NIPT. To what extent the total number of sequencing reads influences FF estimate has been unexplored. In this study, to test the robustness of SeqFF FF estimation and provide additional recommendations for NIPT analysis quality control, we compared the SeqFF FF estimates with two other methods and investigated how the number of sequencing reads and FF level affects the accuracy and precision of FF estimates. METHODS: WGS data of 516 NIPT samples from a prenatal screening program was obtained. Sample data were randomly downsampled by the read count, and FF was calculated by SeqFF software. Then, the outcome was compared with FF estimates from SNP- and chrY-based methods. FF estimated with different read counts and FF levels were compared with FF at 30 M reads as a reference. RESULTS: SeqFF FF highly correlates with SNP- and chrY-based FF estimates. Raising read count from 2 M to 10 M drastically increased the accuracy of FF estimates. After adding more reads, we saw a further improvement in FF accuracy, reaching a plateau at 20 M reads. Precision of SeqFF FF estimate is independent of FF level in the sample. CONCLUSION: SeqFF is a robust method for FF estimation for both genders and for any FF level in range 2-13%. Accuracy of FF estimates highly depends on the read count. We recommend using no less than 10 M reads to achieve accurate FF estimates for NIPT analysis in clinical settings.

Is MED13L-related intellectual disability a recognizable syndrome?

INTRODUCTION: MED13L-related intellectual disability is characterized by moderate intellectual disability (ID), speech impairment, and dysmorphic facial features. We present 8 patients with MED13L-related intellectual disability and review the literature for phenotypical and genetic aspects of previously described patients. MATERIALS AND METHODS: In the search for genetic aberrations in individuals with ID, two of the patients were identified by chromosomal microarray analysis, and five by exome sequencing. One of the individuals, suspected of MED13L-related intellectual disability, based on clinical features, was identified by Sanger sequencing. RESULTS: All 8 individuals had de novo MED13L aberrations, including two intragenic microdeletions, two frameshift, three nonsense variants, and one missense variant. Phenotypically, they all had intellectual disability, speech and motor delay, and features of the mouth (open mouth appearance, macroglossia, and/or macrostomia). Two individuals were diagnosed with autism, and one had autistic features. One had complex congenital heart defect, and one had persistent foramen ovale. The literature was reviewed with respect to clinical and dysmorphic features, and genetic aberrations. CONCLUSIONS: Even if most clinical features of MED13L-related intellectual disability are rather non-specific, the syndrome may be suspected in some individuals based on the association of developmental delay, speech impairment, bulbous nasal tip, and macroglossia, macrostomia, or open mouth appearance.

Identification of metastasis driver genes by massive parallel sequencing of successive steps of breast cancer progression.

Cancer results from alterations at essential genomic sites and is characterized by uncontrolled cell proliferation, invasion and metastasis. Identification of driver genes of metastatic progression is essential, as metastases, not primary tumors, are fatal. To gain insight into the mutational concordance between different steps of malignant progression we performed exome sequencing and validation with targeted deep sequencing of successive steps of malignant progression from pre-invasive stages to asynchronous distant metastases in six breast cancer patients. Using the ratio of non-synonymous to synonymous mutations, a surprisingly large number of cancer driver genes, ranging between 3 and 145, were estimated to confer a selective advantage in the studied primary tumors. We report a substantial amount of metastasis specific mutations and a number of novel putative metastasis driver genes. Most notable are the DCC, ABCA13, TIAM2, CREBBP, BCL6B and ZNF185 genes, mainly mutated exclusively in metastases and highly likely driver genes of metastatic progression. We find different genes and pathways to be affected at different steps of malignant progression. The Adherens junction pathway is affected in four of the six studied patients and this pathway most likely plays a vital role in the metastatic process.

Genomic Analyses of Breast Cancer Progression Reveal Distinct Routes of Metastasis Emergence.

A main controversy in cancer research is whether metastatic abilities are present in the most advanced clone of the primary tumor or result from independently acquired aberrations in early disseminated cancer cells as suggested by the linear and the parallel progression models, respectively. The genetic concordance between different steps of malignant progression is mostly unexplored as very few studies have included cancer samples separated by both space and time. We applied whole exome sequencing and targeted deep sequencing to 26 successive samples from six patients with metastatic estrogen receptor (ER)-positive breast cancer. Our data provide support for both linear and parallel progression towards metastasis. We report for the first time evidence of metastasis-to-metastasis seeding in breast cancer. Our results point to three distinct routes of metastasis emergence. This may have profound clinical implications and provides substantial novel molecular insights into the timing and mutational evolution of breast cancer metastasis.

Bone structure in two adult subjects with impaired minor spliceosome function resulting from RNU4ATAC mutations causing microcephalic osteodysplastic primordial dwarfism type 1 (MOPD1).

Microcephalic osteodysplastic primordial dwarfism type 1 (MOPD1), or Taybi-Linder syndrome is characterized by distinctive skeletal dysplasia, severe intrauterine and postnatal growth retardation, microcephaly, dysmorphic features, and neurological malformations. It is an autosomal recessive disorder caused by homozygous or compound heterozygous mutations in the RNU4ATAC gene resulting in impaired function of the minor spliceosome. Here, we present the first report on bone morphology, bone density and bone microstructure in two adult MOPD1 patients and applied radiographs, dual energy X-ray absorptiometry, high-resolution peripheral quantitative computed tomography and biochemical evaluation. The MOPD1 patients presented with short stature, low BMI but normal macroscopic bone configuration. Bone mineral density was low. Compared to Danish reference data, total bone area, cortical bone area, cortical thickness, total bone density, cortical bone density, trabecular bone density and trabecular bone volume per tissue volume (BV/TV) were all low. These findings may correlate to the short stature and low body weight of the MOPD1 patients. Our findings suggest that minor spliceosome malfunction may be associated with altered bone modelling.

DamX Controls Reversible Cell Morphology Switching in Uropathogenic Escherichia coli.

UNLABELLED: The ability to change cell morphology is an advantageous characteristic adopted by multiple pathogenic bacteria in order to evade host immune detection and assault during infection. Uropathogenic Escherichia coli (UPEC) exhibits such cellular dynamics and has been shown to transition through a series of distinct morphological phenotypes during a urinary tract infection. Here, we report the first systematic spatio-temporal gene expression analysis of the UPEC transition through these phenotypes by using a flow chamber-based in vitro infection model that simulates conditions in the bladder. This analysis revealed a novel association between the cell division gene damX and reversible UPEC filamentation. We demonstrate a lack of reversible bacterial filamentation in a damX deletion mutant in vitro and absence of a filamentous response by this mutant in a murine model of cystitis. While deletion of damX abrogated UPEC filamentation and secondary surface colonization in tissue culture and in mouse infections, transient overexpression of damX resulted in reversible UPEC filamentation. In this study, we identify a hitherto-unknown damX-mediated mechanism underlying UPEC morphotypical switching. Murine infection studies showed that DamX is essential for establishment of a robust urinary tract infection, thus emphasizing its role as a mediator of virulence. Our study demonstrates the value of an in vitro methodology, in which uroepithelium infection is closely simulated, when undertaking targeted investigations that are challenging to perform in animal infection models. IMPORTANCE: Urinary tract infections (UTIs) are most often caused by uropathogenic Escherichia coli (UPEC) and account for a considerable health care burden. UPEC exhibits a dynamic lifestyle in the course of infection, in which the bacterium transiently adopts alternative morphologies ranging from rod shaped to coccoid and filamentous, rendering it better at immune evasion and host epithelium adhesion. This penchant for morphotype switching might in large measure account for UPEC's success as a pathogen. In aiming to uncover genes underlying the phenomenon of UPEC morphotype switching, this study identifies damX, a cell division gene, as a mediator of reversible filamentation during UTI. DamX-mediated filamentation represents an additional pathway for bacterial cell shape control, an alternative to SulA-mediated FtsZ sequestration during E. coli uropathogenesis, and hence represents a potential target for combating UTI.

Molecular Concordance Between Primary Breast Cancer and Matched Metastases.

Clinical management of breast cancer is increasingly personalized and based on molecular profiling. Often, primary tumors are used as proxies for systemic disease at the time of recurrence. However, recent studies have revealed substantial discordances between primary tumors and metastases, both with respect to traditional clinical treatment targets and on the genomic and transcriptomic level. With the increasing use of molecularly targeted therapy, discordance of actionable molecular targets between primary tumors and recurrences can result in nonoptimal treatment or unnecessary side effects. The purpose of this review is to illuminate the extent of cancer genome evolution through disease progression and the degree of molecular concordance between primary breast cancers and matched metastases. We present an overview of the most prominent studies investigating the expression of endocrine receptors, transcriptomics, and genome aberrations in primary tumors and metastases. In conclusion, biopsy of metastatic lesions at recurrence of breast cancer is encouraged to provide optimal treatment of the disease. Furthermore, molecular profiling of metastatic tissue provides invaluable mechanistic insight into the biology underlying metastatic progression and has the potential to identify novel, potentially druggable, drivers of progression.

Evaluation of Nine Somatic Variant Callers for Detection of Somatic Mutations in Exome and Targeted Deep Sequencing Data.

Next generation sequencing is extensively applied to catalogue somatic mutations in cancer, in research settings and increasingly in clinical settings for molecular diagnostics, guiding therapy decisions. Somatic variant callers perform paired comparisons of sequencing data from cancer tissue and matched normal tissue in order to detect somatic mutations. The advent of many new somatic variant callers creates a need for comparison and validation of the tools, as no de facto standard for detection of somatic mutations exists and only limited comparisons have been reported. We have performed a comprehensive evaluation using exome sequencing and targeted deep sequencing data of paired tumor-normal samples from five breast cancer patients to evaluate the performance of nine publicly available somatic variant callers: EBCall, Mutect, Seurat, Shimmer, Indelocator, Somatic Sniper, Strelka, VarScan 2 and Virmid for the detection of single nucleotide mutations and small deletions and insertions. We report a large variation in the number of calls from the nine somatic variant callers on the same sequencing data and highly variable agreement. Sequencing depth had markedly diverse impact on individual callers, as for some callers, increased sequencing depth highly improved sensitivity. For SNV calling, we report EBCall, Mutect, Virmid and Strelka to be the most reliable somatic variant callers for both exome sequencing and targeted deep sequencing. For indel calling, EBCall is superior due to high sensitivity and robustness to changes in sequencing depths.

Global gene expression profiling of telangiectasial tissue from patients with hereditary hemorrhagic telangiectasia.

UNLABELLED: Hereditary hemorrhagic telangiectasia (HHT), the most common inherited vascular disorder, is predominantly caused by mutations in ENG and ACVRL1, which are part of the transforming growth factor beta (TGF-ß) signaling pathway. HHT is characterized by the presence of mucocutaneous telangiectases and arteriovenous malformations in visceral organs, primarily the lungs, brain and liver. The most common symptom in HHT is epistaxis originating from nasal telangiectasia, which can be difficult to prevent and can lead to severe anemia. The clinical manifestations of HHT are extremely variable, even within family members, and the exact mechanism of how endoglin and ALK1 haploinsufficiency leads to HHT manifestations remains to be identified. OBJECTIVES: The purpose of this study was to detect significantly differentially regulated genes in HHT, and try to elucidate the pathways and regulatory mechanisms occurring in the affected tissue of HHT patients, in order to further characterize this disorder and hypothesize on how telangiectases develop. By microarray technology (Agilent G3 Human GE 8x60), we performed global gene expression profiling of mRNA transcripts from HHT nasal telangiectasial (n = 40) and non-telangiectasial (n = 40) tissue using a paired design. Comparing HHT telangiectasial and non-telangiectasial tissue, significantly differentially expressed genes were detected using a paired t-test. Gene set analysis was performed using GSA-SNP. In the group of ENG mutation carriers, we detected 67 differentially expressed mRNAs, of which 62 were down-regulated in the telangiectasial tissue. Gene set analysis identified the gene ontology (GO) terms vasculogenesis, TGF-ß signaling, and Wnt signaling as differentially expressed in HHT1. Altered Wnt signaling might be related to HHT pathogenesis and a greater understanding of this may lead to the discovery of therapeutic targets in HHT.

Clonal expansion and linear genome evolution through breast cancer progression from pre-invasive stages to asynchronous metastasis.

Evolution of the breast cancer genome from pre-invasive stages to asynchronous metastasis is complex and mostly unexplored, but highly demanded as it may provide novel markers for and mechanistic insights in cancer progression. The increasing use of personalized therapy of breast cancer necessitates knowledge of the degree of genomic concordance between different steps of malignant progression as primary tumors often are used as surrogates of systemic disease. Based on exome sequencing we performed copy number profiling and point mutation detection on successive steps of breast cancer progression from one breast cancer patient, including two different regions of Ductal Carcinoma In Situ (DCIS), primary tumor and an asynchronous metastasis. We identify a remarkable landscape of somatic mutations, retained throughout breast cancer progression and with new mutational events emerging at each step. Our data, contrary to the proposed model of early dissemination of metastatic cells and parallel progression of primary tumors and metastases, provide evidence of linear progression of breast cancer with relatively late dissemination from the primary tumor. The genomic discordance between the different stages of tumor evolution in this patient emphasizes the importance of molecular profiling of metastatic tissue directing molecularly targeted therapy at recurrence.

Immunodeficiency associated with a nonsense mutation of IKBKB.

We report an infant of consanguineous parents of Turkish decent with a novel immunodeficiency associated with homozygosity for a nonsense mutation of the gene encoding Inhibitor of nuclear factor kappa-B (NF-κB) kinase subunit beta (IKKß). At five months, she presented with respiratory insufficiency and Pneumocystis jirovecii pneumonia which was successfully treated. At nine months, iatrogenic systemic infection with Mycobacterium bovis was found and eventually led to her death at age 14 months. Laboratory findings were reminiscent of hyper-IgM syndrome, but genetic testing gave no explanation before whole exome sequencing revealed a novel mutation abrogating signaling through the canonical NF-κB pathway.

Long non-coding RNA expression profiles in hereditary haemorrhagic telangiectasia.

Hereditary Haemorrhagic Telangiectasia (HHT) is an autosomal dominantly inherited vascular disease characterized by the presence of mucocutaneous telangiectasia and arteriovenous malformations in visceral organs. HHT is predominantly caused by mutations in ENG and ACVRL1, which both belong to the TGF-ß signalling pathway. The exact mechanism of how haploinsufficiency of ENG and ACVRL1 leads to HHT manifestations remains to be identified. As long non-coding RNAs (lncRNAs) are increasingly recognized as key regulators of gene expression and constitute a sizable fraction of the human transcriptome, we wanted to assess whether lncRNAs play a role in the molecular pathogenesis of HHT manifestations. By microarray technology, we profiled lncRNA transcripts from HHT nasal telangiectasial and non-telangiectasial tissue using a paired design. The microarray probes were annotated using the GENCODE v.16 dataset, identifying 4,810 probes mapping to 2,811 lncRNAs. Comparing HHT telangiectasial tissue with HHT non-telangiectasial tissue, we identified 42 lncRNAs that are differentially expressed (q<0.001). Using GREAT, a tool that assumes cis-regulation, we showed that differently expressed lncRNAs are enriched for genomic loci involved in key pathways concerning HHT. Our study identified lncRNAs that are aberrantly expressed in HHT telangiectasia and indicates that lncRNAs may contribute to regulate protein-coding loci in HHT. These results suggest that the lncRNA component of the transcriptome deserves more attention in HHT. A deeper understanding of lncRNAs and their role in telangiectasia formation possesses potential for discovering therapeutic targets in HHT.