Scientific Publications

Abstract +

Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent{\textendash}child trio data. We present 64 assembled haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average contig N50: 26 Mbp) integrate all forms of genetic variation even across complex loci. We identify 107,590 structural variants (SVs), of which 68\% are not discovered by short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterize 130 of the most active mobile element source elements and find that 63\% of all SVs arise by homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1,526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population.

Abstract +

We report monozygotic twin girls with syndromic intellectual disability who underwent exome sequencing but with negative pathogenic variants. To search for variants that are unrecognized by exome sequencing, high-fidelity long-read genome sequencing (HiFi LR-GS) was applied. A 12-kb copy-neutral inversion was precisely identified by HiFi LR-GS after trio-based variant filtering. This inversion directly disrupted two genes, CPNE9 and BRPF1, the latter of which attracted our attention because pathogenic BRPF1 variants have been identified in autosomal dominant intellectual developmental disorder with dysmorphic facies and ptosis (IDDDFP), which later turned out to be clinically found in the twins. Trio-based HiFi LR-GS together with haplotype phasing revealed that the 12-kb inversion occurred de novo on the maternally transmitted chromosome. This study clearly indicates that submicroscopic copy-neutral inversions are important but often uncharacterized culprits in monogenic disorders and that long-read sequencing is highly advantageous for detecting such inversions involved in genetic diseases.

Abstract +

Long-read sequencing (LRS) has the potential to comprehensively identify all medically relevant genome variation, including variation commonly missed by short-read sequencing (SRS) approaches. To determine this potential, we performed LRS around 15×–40× genome coverage using the Pacific Biosciences Sequel I System for five trios. The respective probands were diagnosed with intellectual disability (ID) whose etiology remained unresolved after SRS exomes and genomes. Systematic assessment of LRS coverage showed that ~35 Mb of the human reference genome was only accessible by LRS and not SRS. Genome-wide structural variant (SV) calling yielded on average 28,292 SV calls per individual, totaling 12.9 Mb of sequence. Trio-based analyses which allowed to study segregation, showed concordance for up to 95% of these SV calls across the genome, and 80% of the LRS SV calls were not identified by SRS. De novo mutation analysis did not identify any de novo SVs, confirming that these are rare events. Because of high sequence coverage, we were also able to call single nucleotide substitutions. On average, we identified 3 million substitutions per genome, with a Mendelian inheritance concordance of up to 97%. Of these, ~100,000 were located in the ~35 Mb of the genome that was only captured by LRS. Moreover, these variants affected the coding sequence of 64 genes, including 32 known Mendelian disease genes. Our data show the potential added value of LRS compared to SRS for identifying medically relevant genome variation.

Abstract +

Purpose Exome and genome sequencing have proven to be effective tools for the diagnosis of neurodevelopmental disorders (NDDs), but large fractions of NDDs cannot be attributed to currently detectable genetic variation. This is likely, at least in part, a result of the fact that many genetic variants are difficult or impossible to detect through typical short-read sequencing approaches.Methods Here, we describe a genomic analysis using Pacific Biosciences circular consensus sequencing (CCS) reads, which are both long (>10 kb) and accurate (>99% bp accuracy). We used CCS on six proband-parent trios with NDDs that were unexplained despite extensive testing, including genome sequencing with short reads.Results We identified variants and created de novo assemblies in each trio, with global metrics indicating these data sets are more accurate and comprehensive than those provided by short-read data. In one proband, we identified a likely pathogenic (LP), de novo L1-mediated insertion in CDKL5 that results in duplication of exon 3, leading to a frameshift. In a second proband, we identified multiple large de novo structural variants, including insertion-translocations affecting DGKB and MLLT3, which we show disrupt MLLT3 transcript levels. We consider this extensive structural variation likely pathogenic.Conclusion The breadth and quality of variant detection, coupled to finding variants of clinical and research interest in two of six probands with unexplained NDDs strongly support the value of long-read genome sequencing for understanding rare disease.Competing Interest StatementThe authors have declared no competing interest.

Abstract +

Potato virus Y (PVY) is a major potato (Solanum tuberosum L.) pathogen that causes severe annual crop losses worth billions of dollars worldwide. PVY is transmitted by aphids, and successful control of virus transmission requires the extensive use of environmentally damaging insecticides to reduce vector populations. Rysto , from the wild relative S. stoloniferum, confers extreme resistance (ER) to PVY and related viruses and is a valuable trait that is widely employed in potato resistance breeding programmes. Rysto was previously mapped to a region of potato chromosome XII, but the specific gene has not been identified to date. In this study, we isolated Rysto using resistance gene enrichment sequencing (RenSeq) and PacBio SMRT (Pacific Biosciences single-molecule real-time sequencing). Rysto was found to encode a nucleotide-binding leucine-rich repeat (NLR) protein with an N-terminal TIR domain and was sufficient for PVY perception and ER in transgenic potato plants. Rysto -dependent extreme resistance was temperature-independent and requires EDS1 and NRG1 proteins. Rysto may prove valuable for creating PVY-resistant cultivars of potato and other Solanaceae crops. © 2019 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.

Abstract +

ESBL and AmpC ß-lactamases are an increasing concern for public health. Studies suggest that ESBL/pAmpC-producing Escherichia coli and their plasmids carrying antibiotic resistance genes can spread from broilers to humans working or living on broiler farms. These studies used traditional typing methods, which may not have provided sufficient resolution to reliably assess the relatedness of these isolates.Eleven suspected transmission events among broilers and humans living/working on eight broiler farms were investigated using whole-genome short-read (Illumina) and long-read sequencing (PacBio). Core genome MLST (cgMLST) was performed to investigate the occurrence of strain transmission. Horizontal plasmid and gene transfer were analysed using BLAST.Of eight suspected strain transmission events, six were confirmed. The isolate pairs had identical ESBL/AmpC genes and fewer than eight allelic differences according to the cgMLST, and five had an almost identical plasmid composition. On one of the farms, cgMLST revealed that the isolate pairs belonging to ST10 from a broiler and a household member of the farmer had 475 different alleles, but that the plasmids were identical, indicating horizontal transfer of mobile elements rather than strain transfer. Of three suspected horizontal plasmid transmission events, one was confirmed. In addition, gene transfer between plasmids was found.The present study confirms transmission of strains as well as horizontal plasmid and gene transfer between broilers and farmers and household members on the same farm. WGS is an important tool to confirm suspected zoonotic strain and resistance gene transmission. © The Author(s) 2019. Published by Oxford University Press on behalf of the British Society for Antimicrobial Chemotherapy. All rights reserved. For permissions, please email: journals.permissions@oup.com.

Abstract +

Long-read sequencing has substantial advantages for structural variant discovery and phasing of vari- ants compared to short-read technologies, but the required and optimal read length has not been as- sessed. In this work, we used long reads simulated from human genomes and evaluated structural vari- ant discovery and variant phasing using current best practicebioinformaticsmethods.Wedeterminedthatoptimal discovery of structural variants from human genomes can be obtained with reads of minimally 20 kb. Haplotyping variants across genes only reaches its optimum from reads of 100 kb. These findings are important for the design of future long-read sequenc- ing projects.

Abstract +

Cheese is a fermented dairy product that is popular for its unique flavor and nutritional value. Recent studies have shown that microorganisms in cheese play an important role in the fermentation process and determine the quality of the cheese. We collected 12 cheese samples from different regions and studied the composition of their bacterial communities using PacBio small-molecule real-time sequencing (Pacific Biosciences, Menlo Park, CA). Our data revealed 144 bacterial genera (including Lactobacillus, Streptococcus, Lactococcus, and Staphylococcus) and 217 bacterial species (including Lactococcus lactis, Streptococcus thermophilus, Staphylococcus equorum, and Streptococcus uberis). We investigated the flavor quality of the cheese samples using an electronic nose system and we found differences in flavor-quality indices among samples from different regions. We found a clustering tendency based on flavor quality using principal component analysis. We found correlations between lactic acid bacteria and the flavor quality of the cheese samples. Biodegradation and metabolism of xenobiotics, and lipid-metabolism-related pathways, were predicted to contribute to differences in cheese flavor using Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt). This preliminary study explored the bacterial communities in cheeses collected from different regions and their potential genome functions from the perspective of flavor quality.Copyright © 2020 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.

Abstract +

Carbapenem-resistant K. pneumoniae 2297, isolated from a patient treated with tigecycline for pneumonia, developed tigecycline resistance, in contrast to carbapenem-resistant isolate 1215, which was collected four months prior to the 2297 isolate. Mechanisms underlying tigecycline resistance were elucidated for the clinical isolates.The tigecycline minimum inhibitory concentration (MIC) was determined using the broth microdilution method, with or without phenylalanine-arginine ß-naphthylamide (PABN), and whole-genome sequencing was carried out by single-molecule real-time sequencing. The expression levels of the genes acrA, oqxA, ramA, rarA, and rpoB were determined by reverse-transcription quantitative PCR.Both isolates presented identical antibiograms, except for tigecycline, which showed an MIC of 0.5 mg/L in 1215 and 2 mg/L in 2297. The addition of PABN to tigecycline-resistant 2297 caused a four-fold decrease in the tigecycline MIC to 0.5 mg/L, although acrA expression (encoding the AcrAB efflux pump) was upregulated by 2.5 fold and ramA expression (encoding the pump activator RamA) was upregulated by 1.4 fold. We identified a 6,096-bp fragment insertion flanking direct TATAT repeats that disrupted the romA gene located upstream of ramA in the chromosome of K. pneumoniae 2297; the insertion led the ramA gene promoter replacement resulting in stronger activation of the gene.The K. pneumoniae isolate developed tigecycline resistance during tigecycline treatment. It was related to the overexpression of the AcrAB resistance-nodulation-cell division efflux system due to promoter replacement. © The Korean Society for Laboratory Medicine.

Abstract +

Olive (Olea europaea) is a rich source of valuable bioactive polyphenols, which has attracted widespread interest. In this study, we combined targeted metabolome, Pacbio ISOseq transcriptome, and Illumina RNA-seq transcriptome to investigate the association between polyphenols and gene expression in the developing olive fruits and leaves. A total of 12 main polyphenols were measured, and 122 transcripts of 17 gene families, 101 transcripts of 9 gene families, and 106 transcripts of 6 gene families that encode for enzymes involved in flavonoid, oleuropein, and hydroxytyrosol biosynthesis were separately identified. Additionally, 232 alternative splicing events of 18 genes related to polyphenol synthesis were analyzed. This is the first time that the third generations of full-length transcriptome technology were used to study the gene expression pattern of olive fruits and leaves. The results of transcriptome combined with targeted metabolome can help us better understand the polyphenol biosynthesis pathways in the olive.Copyright © 2019 Elsevier Ltd. All rights reserved.