In this PAGBioDay 2021 presentation, Michelle Vierra of PacBio provides the official welcome to PAGBio Day including a presentation on the latest in HiFi sequencing. Michelle is be followed by…
During the past decade, the search for pathogenic mutations in rare human genetic diseases has involved huge efforts to sequence coding regions, or the entire genome, using massively parallel short-read sequencers. However, the approximate current diagnostic rate is <50% using these approaches, and there remain many rare genetic diseases with unknown cause. There may be many reasons for this, but one plausible explanation is that the responsible mutations are in regions of the genome that are difficult to sequence using conventional technologies (e.g., tandem-repeat expansion or complex chromosomal structural aberrations). Despite the drawbacks of high cost and a shortage of standard analytical methods, several studies have analyzed pathogenic changes in the genome using long-read sequencers. The results of these studies provide hope that further application of long-read sequencers to identify the causative mutations in unsolved genetic diseases may expand our understanding of the human genome and diseases. Such approaches may also be applied to molecular diagnosis and therapeutic strategies for patients with genetic diseases in the future.
The sequence and assembly of human genomes using long-read sequencing technologies has revolutionized our understanding of structural variation and genome organization. We compared the accuracy, continuity, and gene annotation of genome assemblies generated from either high-fidelity (HiFi) or continuous long-read (CLR) datasets from the same complete hydatidiform mole human genome. We find that the HiFi sequence data assemble an additional 10% of duplicated regions and more accurately represent the structure of tandem repeats, as validated with orthogonal analyses. As a result, an additional 5 Mbp of pericentromeric sequences are recovered in the HiFi assembly, resulting in a 2.5-fold increase in the NG50 within 1 Mbp of the centromere (HiFi 480.6 kbp, CLR 191.5 kbp). Additionally, the HiFi genome assembly was generated in significantly less time with fewer computational resources than the CLR assembly. Although the HiFi assembly has significantly improved continuity and accuracy in many complex regions of the genome, it still falls short of the assembly of centromeric DNA and the largest regions of segmental duplication using existing assemblers. Despite these shortcomings, our results suggest that HiFi may be the most effective standalone technology for de novo assembly of human genomes. © 2019 John Wiley & Sons Ltd/University College London.
Zygotic genome activation (ZGA) following fertilization is accomplished through a process termed the maternal-to-zygotic transition, during which the maternal RNAs and proteins are degraded and zygotic genome is transcriptionally activated.1 In mice, minor ZGA occurs from S phase of the zygote to G1 phase of the two-cell (2C) embryo, while major ZGA takes place during the middle-to-late 2C stage with a burst of transcription of totipotent cleavage stage-specific genes and retrotransposons.2Dux has been recently identified and considered as a master inducer that regulates the ZGA process.3–5Dux can directly bind and robustly activate 2C stage-specific ZGA transcripts and convert mouse embryonic stem cells (mESCs) to a 2C-like state with unique features that resembles the 2C embryos.4Intriguingly, ~20% embryos with zygotic depletion of Dux unexpectedly reached morula or blastocyst stage even though defective ZGA program was detected.
N6-methyladenosine (m6A) is a widespread RNA modification that influences nearly every aspect of the messenger RNA lifecycle. Our understanding of m6A has been facilitated by the development of global m6A mapping methods, which use antibodies to immunoprecipitate methylated RNA. However, these methods have several limitations, including high input RNA requirements and cross-reactivity to other RNA modifications. Here, we present DART-seq (deamination adjacent to RNA modification targets), an antibody-free method for detecting m6A sites. In DART-seq, the cytidine deaminase APOBEC1 is fused to the m6A-binding YTH domain. APOBEC1-YTH expression in cells induces C-to-U deamination at sites adjacent to m6A residues, which are detected using standard RNA-seq. DART-seq identifies thousands of m6A sites in cells from as little as 10?ng of total RNA and can detect m6A accumulation in cells over time. Additionally, we use long-read DART-seq to gain insights into m6A distribution along the length of individual transcripts.
The widespread occurrence of repetitive stretches of DNA in genomes of organisms across the tree of life imposes fundamental challenges for sequencing, genome assembly, and automated annotation of genes and proteins. This multi-level problem can lead to errors in genome and protein databases that are often not recognized or acknowledged. As a consequence, end users working with sequences with repetitive regions are faced with ‘ready-to-use’ deposited data whose trustworthiness is difficult to determine, let alone to quantify. Here, we provide a review of the problems associated with tandem repeat sequences that originate from different stages during the sequencing-assembly-annotation-deposition workflow, and that may proliferate in public database repositories affecting all downstream analyses. As a case study, we provide examples of the Atlantic cod genome, whose sequencing and assembly were hindered by a particularly high prevalence of tandem repeats. We complement this case study with examples from other species, where mis-annotations and sequencing errors have propagated into protein databases. With this review, we aim to raise the awareness level within the community of database users, and alert scientists working in the underlying workflow of database creation that the data they omit or improperly assemble may well contain important biological information valuable to others. © The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.
Domestication of clonally propagated crops such as pineapple from South America was hypothesized to be a ‘one-step operation’. We sequenced the genome of Ananas comosus var. bracteatus CB5 and assembled 513?Mb into 25 chromosomes with 29,412 genes. Comparison of the genomes of CB5, F153 and MD2 elucidated the genomic basis of fiber production, color formation, sugar accumulation and fruit maturation. We also resequenced 89 Ananas genomes. Cultivars ‘Smooth Cayenne’ and ‘Queen’ exhibited ancient and recent admixture, while ‘Singapore Spanish’ supported a one-step operation of domestication. We identified 25 selective sweeps, including a strong sweep containing a pair of tandemly duplicated bromelain inhibitors. Four candidate genes for self-incompatibility were linked in F153, but were not functional in self-compatible CB5. Our findings support the coexistence of sexual recombination and a one-step operation in the domestication of clonally propagated crops. This work guides the exploration of sexual and asexual domestication trajectories in other clonally propagated crops.
A Gram-stain-negative bacterial strain, designated CA10T, was isolated from bovine raw milk sampled in Anseong, Republic of Korea. Cells were yellow-pigmented, aerobic, non-motile bacilli and grew optimally at 30?°C and pH 7.0 on tryptic soy agar without supplementation of NaCl. Phylogenetic analysis based on the 16S rRNA gene sequences revealed that strain CA10T belonged to the genus Chryseobacterium, family Flavobacteriaceae, and was most closely related to Chryseobacterium indoltheticum ATCC 27950T (98.75?% similarity). The average nucleotide identity and digital DNA-DNA hybridization values of strain CA10T were 94.4 and 56.9?%, respectively, relative to Chryseobacterium scophthalmum DSM 16779T, being lower than the cut-off values of 95-96?and 70?%, respectively. The predominant respiratory quinone was menaquinone-6; major polar lipid, phosphatidylethanolamine; major fatty acids, iso-C15?:?0, summed feature 9 (iso-C17?:?1?9c and/or C16?:?0 10-methyl), summed feature 3 (iso-C15?:?0 2-OH and/or C16?:?1?7c) and iso-C17?:?0 3-OH. The results of physiological, chemotaxonomic and biochemical analyses suggested that strain CA10T is a novel species of genus Chryseobacterium, for which the name Chryseobacterium mulctrae sp. nov. is proposed. The type strain is CA10T (=KACC 21234T=JCM 33443T).
A Gram-stain-negative, rod-shaped and red-pigmented strain, HME7025T, was isolated from freshwater sampled in the Republic of Korea. Phylogenetic analysis based on its 16S rRNA gene sequence revealed that strain HME7025T formed a lineage within the family Cytophagaceae of the phylum Bacteroidetes. Strain HME7025T was closely related to the genera Pseudarcicella, Arcicella and Flectobacillus. The 16S rRNA gene sequence similarity values of strain HME7025T were under 94.5?% to its closest phylogenetic neighbours. The major fatty acids of strain HME7025T were iso-C15?:?0 (41.9?%), summed feature 3 (comprising C16?:?1?7c and/or C16?:?1?6c; 12.2?%) and anteiso-C15?:?0 (10.8?%). The major respiratory quinone was menaquinone-7. The major polar lipids were phosphatidylethanolamine, two unidentified aminophospholipids and one unidentified polar lipid. The DNA G+C content of strain HME7025T was 37.9?mol%. On the basis of the evidence presented in this study, strain HME7025T represents a novel species of a novel genus within the family Cytophagaceae, for which the name Allopseudarcicella aquatilis gen. nov., sp. nov. is proposed. The type strain is HME7025T (=KCTC 23617T=CECT 7957T).
Pharmacogenetic testing increasingly is available from clinical and research laboratories. However, only a limited number of quality control and other reference materials currently are available for the complex rearrangements and rare variants that occur in the CYP2D6 gene. To address this need, the Division of Laboratory Systems, CDC-based Genetic Testing Reference Material Coordination Program, in collaboration with members of the pharmacogenetic testing and research communities and the Coriell Cell Repositories (Camden, NJ), has characterized 179 DNA samples derived from Coriell cell lines. Testing included the recharacterization of 137 genomic DNAs that were genotyped in previous Genetic Testing Reference Material Coordination Program studies and 42 additional samples that had not been characterized previously. DNA samples were distributed to volunteer testing laboratories for genotyping using a variety of commercially available and laboratory-developed tests. These publicly available samples will support the quality-assurance and quality-control programs of clinical laboratories performing CYP2D6 testing.Published by Elsevier Inc.
Antarcticibacterium flavum JB01H24T was isolated from a marine sediment of the Ross Sea, Antarctica. Whole-genome sequencing of the strain Antarcticibacterium flavum JB01H24T was achieved using PacBio RS II platform. The resulting complete genome comprised of one closed, complete chromosome of 4,319,074 base pairs with a 40.87% G?+?C content, where genomic analyses demonstrated that it is constituted mostly by putative ORFs with unknown functions, representing a novel genetic feature. It is the first complete genome sequence of the Antarcticibacterium strain.
Bacillus cereus is an opportunistic human pathogen causing food-borne gastrointestinal infections and non-gastrointestinal infections worldwide. The strain B. cereus FORC_013 was isolated from fried eel. Its genome was completely sequenced by PacBio technology, analyzed and compared with other complete genome sequences of Bacillus to elucidate the distinct pathogenic features of the strain isolated in South Korea. Genomic analysis revealed pathogenesis and host immune evasion-associated genes encoding tissue-destructive exoenzymes, and pore-forming toxins. In particular, tissue-destructive (hemolysin BL, nonhaemolytic enterotoxins) and cytolytic proteins (cytolysin) were observed in the genome, which damage the plasma membrane of the epithelial cells of the small intestine causing diarrhea in humans. Capsule biosynthesis gene found in both chromosome and plasmid, which might be responsible for protecting the pathogen from the host cell immune defense system after host cell invasion. Additionally, multidrug resistance operon and efflux pumps were identified in the genome, which play a prominent role in multi-antibiotic resistance. Comparative phylogenetic tree analysis of the strain FORC_013 and other B. cereus strains revealed that the closest strains are ATCC 14579 and B4264. This genome data can be used to identify virulence factors that can be applied for the development of novel biomarkers for the rapid detection of this pathogen in foods.Copyright © 2018. Published by Elsevier Inc.
To track stepwise changes in genetic diversity and antimicrobial resistance in rapidly evolving OXA-232-producing Klebsiella pneumoniae ST14, an emerging carbapenem-resistant high-risk clone, in clinical settings.Twenty-six K. pneumoniae ST14 isolates were collected by the Korean Nationwide Surveillance of Antimicrobial Resistance system over the course of 1 year. Isolates were subjected to whole-genome sequencing and MIC determinations using 33 antibiotics from 14 classes.Single-nucleotide polymorphism (SNP) typing identified 72 unique SNP sites spanning the chromosomes of the isolates, dividing them into three clusters (I, II and III). The initial isolate possessed two plasmids with 18 antibiotic-resistance genes, including blaOXA-232, and exhibited resistance to 11 antibiotic classes. Four other plasmids containing 12 different resistance genes, including blaCTX-M-15 and strA/B, were introduced over time, providing additional resistance to aztreonam and streptomycin. Moreover, chromosomal integration of insertion sequence Ecp1-blaCTX-M-15 mediated the inactivation of mgrB responsible for colistin resistance in four isolates from cluster III. To the best of our knowledge, this is the first description of K. pneumoniae ST14 resistant to both carbapenem and colistin in South Korea. Furthermore, although some acquired genes were lost over time, the retention of 12 resistance genes and inactivation of mgrB provided resistance to 13 classes of antibiotics.We describe stepwise changes in OXA-232-producing K. pneumoniae ST14 in vivo over time in terms of antimicrobial resistance. Our findings contribute to our understanding of the evolution of emerging high-risk K. pneumoniae clones and provide reference data for future outbreaks.Copyright © 2019 European Society of Clinical Microbiology and Infectious Diseases. Published by Elsevier Ltd. All rights reserved.
Circoviruses are found in many species, including mammals, birds, lower vertebrates and invertebrates. To date, there are no reports of circovirus-induced diseases in chickens. In this study, we identified a new strain of chicken circovirus (CCV) by PacBio third-generation sequencing samples from chickens with acute gastroenteritis in a Shandong commercial broiler farm in China. The complete genome of CCV was verified by inverse PCR. Genomic analysis revealed that CCV codes two inverse open reading frames (ORFs), and a potential stem-loop structure was present at the 5′ end with a structure typical of a circular virus. Phylogenetic tree analysis showed that CCV formed an independent branch between mammalian and avian circovirus, and homology analysis indicated that the homology of CCV with 21 other known circoviruses was less than 40%. Thus, this CCV strain represents a new species in the genus Circovirus. The infection rate of CCV in 12 chickens with diarrhoea was 100%, but no CCV was found in healthy chickens, thereby indicating that the novel CCV strain is highly associated with acute infectious gastroenteritis in chickens. The emergence of a novel CCV in commercial broiler chickens is highly concerning for the broiler industry. © 2019 Blackwell Verlag GmbH.
The genomes of classical inbred mouse strains include genes derived from all three major subspecies of the house mouse, Mus musculus. We recently posited that genetic diversity in the immunoglobulin heavy chain (IGH) gene loci of C57BL/6 and BALB/c mice reflect differences in subspecies origin. To investigate this hypothesis, we conducted high-throughput sequencing of IGH gene rearrangements to document IGH variable (IGHV), joining (IGHJ), and diversity (IGHD) genes in four inbred wild-derived mouse strains (CAST/EiJ, LEWES/EiJ, MSM/MsJ, and PWD/PhJ), and a single disease model strain (NOD/ShiLtJ), collectively representing genetic backgrounds of several major mouse subspecies. A total of 341 germline IGHV sequences were inferred in the wild-derived strains, including 247 not curated in the International Immunogenetics Information System. In contrast, 83/84 inferred NOD IGHV genes had previously been observed in C57BL/6 mice. Variability among the strains examined was observed for only a single IGHJ gene, involving a description of a novel allele. In contrast, unexpected variation was found in the IGHD gene loci, with four previously unreported IGHD gene sequences being documented. Very few IGHV sequences of C57BL/6 and BALB/c mice were shared with strains representing major subspecies, suggesting that their IGH loci may be complex mosaics of genes of disparate origins. This suggests a similar level of diversity is likely present in the IGH loci of other classical inbred strains. This must now be documented if we are to properly understand inter-strain variation in models of antibody-mediated disease. This article is protected by copyright. All rights reserved.This article is protected by copyright. All rights reserved.
If you have a question, need to check the status of an order, or are interested in purchasing an instrument, we're here to help.