1
|
Agustinho DP, Fu Y, Menon VK, Metcalf GA, Treangen TJ, Sedlazeck FJ. Unveiling microbial diversity: harnessing long-read sequencing technology. Nat Methods 2024; 21:954-966. [PMID: 38689099 DOI: 10.1038/s41592-024-02262-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2022] [Accepted: 03/29/2024] [Indexed: 05/02/2024]
Abstract
Long-read sequencing has recently transformed metagenomics, enhancing strain-level pathogen characterization, enabling accurate and complete metagenome-assembled genomes, and improving microbiome taxonomic classification and profiling. These advancements are not only due to improvements in sequencing accuracy, but also happening across rapidly changing analysis methods. In this Review, we explore long-read sequencing's profound impact on metagenomics, focusing on computational pipelines for genome assembly, taxonomic characterization and variant detection, to summarize recent advancements in the field and provide an overview of available analytical methods to fully leverage long reads. We provide insights into the advantages and disadvantages of long reads over short reads and their evolution from the early days of long-read sequencing to their recent impact on metagenomics and clinical diagnostics. We further point out remaining challenges for the field such as the integration of methylation signals in sub-strain analysis and the lack of benchmarks.
Collapse
Affiliation(s)
- Daniel P Agustinho
- Human Genome Sequencing center, Baylor College of Medicine, Houston, TX, USA
| | - Yilei Fu
- Department of Computer Science, Rice University, Houston, TX, USA
| | - Vipin K Menon
- Human Genome Sequencing center, Baylor College of Medicine, Houston, TX, USA
- Senior research project manager, Human Genetics, Genentech, South San Francisco, CA, USA
| | - Ginger A Metcalf
- Human Genome Sequencing center, Baylor College of Medicine, Houston, TX, USA
| | - Todd J Treangen
- Department of Computer Science, Rice University, Houston, TX, USA
- Department of Bioengineering, Rice University, Houston, TX, USA
| | - Fritz J Sedlazeck
- Human Genome Sequencing center, Baylor College of Medicine, Houston, TX, USA.
- Department of Computer Science, Rice University, Houston, TX, USA.
| |
Collapse
|
2
|
Gofton AW, Popa-Baez A, Takano A, Soennichsen K, Michie M, Short M, Supriyono S, Pascoe J, Cusbert S, Mulley R. Characterisation and comparative genomics of three new Varanus-associated Borrelia spp. from Indonesia and Australia. Parasit Vectors 2023; 16:317. [PMID: 37670353 PMCID: PMC10481545 DOI: 10.1186/s13071-023-05937-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Accepted: 08/16/2023] [Indexed: 09/07/2023] Open
Abstract
BACKGROUND Borrelia are important disease-causing tick- and louse-borne spirochaetes than can infect a wide variety of vertebrates, including humans and reptiles. Reptile-associated (REP) Borrelia, once considered a peculiarity, are now recognised as a distinct and important evolutionary lineage, and are increasingly being discovered worldwide in association with novel hosts. Numerous novel Borrelia spp. associated with monitor lizards (Varanus spp.) have been recently identified throughout the Indo-Pacific region; however, there is a lack of genomic data on these Borrelia. METHODS We used metagenomic techniques to sequence almost complete genomes of novel Borrelia spp. from Varanus varius and Varanus giganteus from Australia, and used long- and short-read technologies to sequence the complete genomes of two strains of a novel Borrelia sp. previously isolated from ticks infesting Varanus salvator from Indonesia. We investigated intra- and interspecies genomic diversity, including plasmid diversity and relatedness, among Varanus-associated Borrelia and other available REP Borrelia and, based on 712 whole genome orthologues, produced the most complete phylogenetic analysis, to the best of our knowledge, of REP Borrelia to date. RESULTS The genomic architecture of Varanus-associated Borrelia spp. is similar to that of Borrelia spp. that cause relapsing fever (RF), and includes a highly conserved megaplasmid and numerous smaller linear and circular plasmids that lack structural consistency between species. Analysis of PF32 and PF57/62 plasmid partitioning genes indicated that REP Borrelia plasmids fall into at least six distinct plasmid families, some of which are related to previously defined Borrelia plasmid families, whereas the others appear to be unique. REP Borrelia contain immunogenic variable major proteins that are homologous to those found in Borrelia spp. that cause RF, although they are limited in copy number and variability and have low sequence identities to RF variable major proteins. Phylogenetic analyses based on single marker genes and 712 single copy orthologs also definitively demonstrated the monophyly of REP Borrelia as a unique lineage. CONCLUSIONS In this work we present four new genomes from three novel Borrelia, and thus double the number of REP Borrelia genomes publicly available. The genomic characterisation of these Borrelia clearly demonstrates their distinctiveness as species, and we propose the names Borrelia salvatorii, 'Candidatus Borrelia undatumii', and 'Candidatus Borrelia rubricentralis' for them.
Collapse
Affiliation(s)
- Alexander William Gofton
- Health and Biosecurity, Commonwealth Scientific and Industrial Research Organisation (CSIRO), Canberra, Australia
| | - Angel Popa-Baez
- Health and Biosecurity, Commonwealth Scientific and Industrial Research Organisation (CSIRO), Canberra, Australia
| | - Ai Takano
- Department of Veterinary Medicine, Joint Faculty of Veterinary Medicine, Yamaguchi University, Yamaguchi, Japan
| | - Kari Soennichsen
- Institute for Applied Ecology, University of Canberra, Canberra, Australia
| | - Michelle Michie
- Health and Biosecurity, Commonwealth Scientific and Industrial Research Organisation (CSIRO), Canberra, Australia
| | - Makenna Short
- Health and Biosecurity, Commonwealth Scientific and Industrial Research Organisation (CSIRO), Canberra, Australia
| | - Supriyono Supriyono
- Department of Animal Diseases and Veterinary Health, Bogor Agricultural University, Bogor, Indonesia
| | - Jack Pascoe
- School of Agricultural and Ecosystem Sciences, University of Melbourne, Melbourne, Australia
| | - Sue Cusbert
- School of Science and Health, Western Sydney University, Penrith, Australia
| | - Robert Mulley
- School of Science and Health, Western Sydney University, Penrith, Australia
| |
Collapse
|
3
|
Lee YC, Ke HM, Liu YC, Lee HH, Wang MC, Tseng YC, Kikuchi T, Tsai IJ. Single-worm long-read sequencing reveals genome diversity in free-living nematodes. Nucleic Acids Res 2023; 51:8035-8047. [PMID: 37526286 PMCID: PMC10450198 DOI: 10.1093/nar/gkad647] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2023] [Revised: 07/10/2023] [Accepted: 07/21/2023] [Indexed: 08/02/2023] Open
Abstract
Obtaining sufficient genetic material from a limited biological source is currently the primary operational bottleneck in studies investigating biodiversity and genome evolution. In this study, we employed multiple displacement amplification (MDA) and Smartseq2 to amplify nanograms of genomic DNA and mRNA, respectively, from individual Caenorhabditis elegans. Although reduced genome coverage was observed in repetitive regions, we produced assemblies covering 98% of the reference genome using long-read sequences generated with Oxford Nanopore Technologies (ONT). Annotation with the sequenced transcriptome coupled with the available assembly revealed that gene predictions were more accurate, complete and contained far fewer false positives than de novo transcriptome assembly approaches. We sampled and sequenced the genomes and transcriptomes of 13 nematodes from early-branching species in Chromadoria, Dorylaimia and Enoplia. The basal Chromadoria and Enoplia species had larger genome sizes, ranging from 136.6 to 738.8 Mb, compared with those in the other clades. Nine mitogenomes were fully assembled, and displayed a complete lack of synteny to other species. Phylogenomic analyses based on the new annotations revealed strong support for Enoplia as sister to the rest of Nematoda. Our result demonstrates the robustness of MDA in combination with ONT, paving the way for the study of genome diversity in the phylum Nematoda and beyond.
Collapse
Affiliation(s)
- Yi-Chien Lee
- Biodiversity Research Center, Academia Sinica, Taipei 115, Taiwan
- Biodiversity Program, Taiwan International Graduate Program, Academia Sinica and National Taiwan Normal University, Taipei, Taiwan
- Department of Life Science, National Taiwan Normal University, 116 Wenshan, Taipei, Taiwan
| | - Huei-Mien Ke
- Department of Microbiology, Soochow University, Taipei, Taiwan
| | - Yu-Ching Liu
- Biodiversity Research Center, Academia Sinica, Taipei 115, Taiwan
| | - Hsin-Han Lee
- Biodiversity Research Center, Academia Sinica, Taipei 115, Taiwan
| | - Min-Chen Wang
- Marine Research Station (MRS), Institute of Cellular and Organismic Biology, Academia Sinica, 262 I-Lan County, Taiwan
| | - Yung-Che Tseng
- Marine Research Station (MRS), Institute of Cellular and Organismic Biology, Academia Sinica, 262 I-Lan County, Taiwan
| | - Taisei Kikuchi
- Department of Integrated Biosciences, Graduate School of Frontier Sciences, The University of Tokyo, Chiba 277-8562, Japan
| | - Isheng Jason Tsai
- Biodiversity Research Center, Academia Sinica, Taipei 115, Taiwan
- Biodiversity Program, Taiwan International Graduate Program, Academia Sinica and National Taiwan Normal University, Taipei, Taiwan
| |
Collapse
|
4
|
Tomaž Š, Petek M, Lukan T, Pogačar K, Stare K, Teixeira Prates E, Jacobson DA, Zrimec J, Bajc G, Butala M, Pompe Novak M, Dudley Q, Patron N, Taler-Verčič A, Usenik A, Turk D, Prat S, Coll A, Gruden K. A mini-TGA protein modulates gene expression through heterogeneous association with transcription factors. PLANT PHYSIOLOGY 2023; 191:1934-1952. [PMID: 36517238 PMCID: PMC10022624 DOI: 10.1093/plphys/kiac579] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/09/2022] [Accepted: 11/24/2022] [Indexed: 06/17/2023]
Abstract
TGA (TGACG-binding) transcription factors, which bind their target DNA through a conserved basic region leucine zipper (bZIP) domain, are vital regulators of gene expression in salicylic acid (SA)-mediated plant immunity. Here, we investigated the role of StTGA2.1, a potato (Solanum tuberosum) TGA lacking the full bZIP, which we named a mini-TGA. Such truncated proteins have been widely assigned as loss-of-function mutants. We, however, confirmed that StTGA2.1 overexpression compensates for SA-deficiency, indicating a distinct mechanism of action compared with model plant species. To understand the underlying mechanisms, we showed that StTGA2.1 can physically interact with StTGA2.2 and StTGA2.3, while its interaction with DNA was not detected. We investigated the changes in transcriptional regulation due to StTGA2.1 overexpression, identifying direct and indirect target genes. Using in planta transactivation assays, we confirmed that StTGA2.1 interacts with StTGA2.3 to activate StPRX07, a member of class III peroxidases (StPRX), which are known to play role in immune response. Finally, via structural modeling and molecular dynamics simulations, we hypothesized that the compact molecular architecture of StTGA2.1 distorts DNA conformation upon heterodimer binding to enable transcriptional activation. This study demonstrates how protein truncation can lead to distinct functions and that such events should be studied carefully in other protein families.
Collapse
Affiliation(s)
| | - Marko Petek
- Department of Biotechnology and Systems Biology, National Institute of Biology, 1000 Ljubljana, Slovenia
| | - Tjaša Lukan
- Department of Biotechnology and Systems Biology, National Institute of Biology, 1000 Ljubljana, Slovenia
| | - Karmen Pogačar
- Department of Biotechnology and Systems Biology, National Institute of Biology, 1000 Ljubljana, Slovenia
| | - Katja Stare
- Department of Biotechnology and Systems Biology, National Institute of Biology, 1000 Ljubljana, Slovenia
| | - Erica Teixeira Prates
- Biosciences Division, Oak Ridge National Laboratory,, Oak Ridge, Tennessee 37831, USA
| | - Daniel A Jacobson
- Biosciences Division, Oak Ridge National Laboratory,, Oak Ridge, Tennessee 37831, USA
| | - Jan Zrimec
- Department of Biotechnology and Systems Biology, National Institute of Biology, 1000 Ljubljana, Slovenia
| | - Gregor Bajc
- Department of Biology, Biotechnical Faculty, University of Ljubljana, 1000 Ljubljana, Slovenia
| | - Matej Butala
- Department of Biology, Biotechnical Faculty, University of Ljubljana, 1000 Ljubljana, Slovenia
| | - Maruša Pompe Novak
- Department of Biotechnology and Systems Biology, National Institute of Biology, 1000 Ljubljana, Slovenia
- School for Viticulture and Enology, University of Nova Gorica, 5271 Vipava, Slovenia
| | - Quentin Dudley
- Earlham Institute, Norwich Research Park, Norwich NR4 7UZ, UK
| | - Nicola Patron
- Earlham Institute, Norwich Research Park, Norwich NR4 7UZ, UK
| | - Ajda Taler-Verčič
- Department of Biochemistry and Molecular and Structural Biology, Jožef Stefan Institute, 1000 Ljubljana, Slovenia
- Faculty of Medicine, Institute of Biochemistry and Molecular Genetics, University of Ljubljana, 1000 Ljubljana, Slovenia
| | - Aleksandra Usenik
- Department of Biochemistry and Molecular and Structural Biology, Jožef Stefan Institute, 1000 Ljubljana, Slovenia
- Centre of Excellence for Integrated Approaches in Chemistry and Biology of Proteins, 1000 Ljubljana, Slovenia
| | - Dušan Turk
- Department of Biochemistry and Molecular and Structural Biology, Jožef Stefan Institute, 1000 Ljubljana, Slovenia
- Centre of Excellence for Integrated Approaches in Chemistry and Biology of Proteins, 1000 Ljubljana, Slovenia
| | - Salomé Prat
- Department of Plant Development and Signal Transduction, Centre for Research in Agricultural Genomics, 08193 Cerdanyola, Barcelona, Spain
| | | | | |
Collapse
|
5
|
Chen H, Durinck S, Patel H, Foreman O, Mesh K, Eastham J, Caothien R, Newman RJ, Roose-Girma M, Darmanis S, Warming S, Lattanzi A, Liang Y, Haley B. Population-wide gene disruption in the murine lung epithelium via AAV-mediated delivery of CRISPR-Cas9 components. MOLECULAR THERAPY - METHODS & CLINICAL DEVELOPMENT 2022; 27:431-449. [DOI: 10.1016/j.omtm.2022.10.016] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/14/2022] [Accepted: 10/27/2022] [Indexed: 11/13/2022]
|
6
|
Nye DB, Tanner NA. Chimeric DNA byproducts in strand displacement amplification using the T7 replisome. PLoS One 2022; 17:e0273979. [PMID: 36121810 PMCID: PMC9484634 DOI: 10.1371/journal.pone.0273979] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2022] [Accepted: 08/19/2022] [Indexed: 11/30/2022] Open
Abstract
Recent advances in next generation sequencing technologies enable reading DNA molecules hundreds of kilobases in length and motivate development of DNA amplification methods capable of producing long amplicons. In vivo, DNA replication is performed not by a single polymerase enzyme, but multiprotein complexes called replisomes. Here, we investigate strand-displacement amplification reactions using the T7 replisome, a macromolecular complex of a helicase, a single-stranded DNA binding protein, and a DNA polymerase. The T7 replisome may initiate processive DNA synthesis from DNA nicks, and the reaction of a 48 kilobase linear double stranded DNA substrate with the T7 replisome and nicking endonucleases is shown to produce discrete DNA amplicons. To gain a mechanistic understanding of this reaction, we utilized Oxford Nanopore long-read sequencing technology. Sequence analysis of the amplicons revealed chimeric DNA reads and uncovered a connection between template switching and polymerase exonuclease activity. Nanopore sequencing provides insight to guide the further development of isothermal amplification methods for long DNA, and our results highlight the need for high-specificity, high-turnover nicking endonucleases to initiate DNA amplification without thermal denaturation.
Collapse
Affiliation(s)
- Dillon B. Nye
- Nucleic Acid Replication Division, New England Biolabs Inc., Ipswich, Massachusetts, United States of America
| | - Nathan A. Tanner
- Nucleic Acid Replication Division, New England Biolabs Inc., Ipswich, Massachusetts, United States of America
- * E-mail:
| |
Collapse
|
7
|
Zhang Z, An HH, Vege S, Hu T, Zhang S, Mosbruger T, Jayaraman P, Monos D, Westhoff CM, Chou ST. Accurate long-read sequencing allows assembly of the duplicated RHD and RHCE genes harboring variants relevant to blood transfusion. Am J Hum Genet 2022; 109:180-191. [PMID: 34968422 DOI: 10.1016/j.ajhg.2021.12.003] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2021] [Accepted: 12/07/2021] [Indexed: 12/18/2022] Open
Abstract
Next-generation sequencing (NGS) technologies have transformed medical genetics. However, short-read lengths pose a limitation on identification of structural variants, sequencing repetitive regions, phasing of distant nucleotide changes, and distinguishing highly homologous genomic regions. Long-read sequencing technologies may offer improvements in the characterization of genes that are currently difficult to assess. We used a combination of targeted DNA capture, long-read sequencing, and a customized bioinformatics pipeline to fully assemble the RH region, which harbors variation relevant to red cell donor-recipient mismatch, particularly among patients with sickle cell disease. RHD and RHCE are a pair of duplicated genes located within an ∼175 kb region on human chromosome 1 that have high sequence similarity and frequent structural variations. To achieve the assembly, we utilized palindrome repeats in PacBio SMRT reads to obtain consensus sequences of 2.1 to 2.9 kb average length with over 99% accuracy. We used these long consensus sequences to identify 771 assembly markers and to phase the RHD-RHCE region with high confidence. The dataset enabled direct linkage between coding and intronic variants, phasing of distant SNPs to determine RHD-RHCE haplotypes, and identification of known and novel structural variations along with the breakpoints. A limiting factor in phasing is the frequency of heterozygous assembly markers and therefore was most successful in samples from African Black individuals with increased heterogeneity at the RH locus. Overall, this approach allows RH genotyping and de novo assembly in an unbiased and comprehensive manner that is necessary to expand application of NGS technology to high-resolution RH typing.
Collapse
Affiliation(s)
- Zhe Zhang
- Department of Biomedical and Health Informatics, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA
| | - Hyun Hyung An
- Division of Hematology, Department of Pediatrics, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA
| | - Sunitha Vege
- Immunohematology and Genomics, New York Blood Center, New York, NY 11101, USA
| | - Taishan Hu
- Immunogenetics Laboratory, Department of Pathology and Laboratory Medicine, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA
| | - Shiping Zhang
- Department of Biomedical and Health Informatics, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA
| | - Timothy Mosbruger
- Immunogenetics Laboratory, Department of Pathology and Laboratory Medicine, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA
| | - Pushkala Jayaraman
- Immunogenetics Laboratory, Department of Pathology and Laboratory Medicine, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA
| | - Dimitri Monos
- Immunogenetics Laboratory, Department of Pathology and Laboratory Medicine, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA; Department of Pathology and Laboratory Medicine, Perelman Schools of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Connie M Westhoff
- Immunohematology and Genomics, New York Blood Center, New York, NY 11101, USA
| | - Stella T Chou
- Division of Hematology, Department of Pediatrics, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA; Division of Transfusion Medicine, Department of Pathology and Laboratory Medicine, The Children's Hospital of Philadelphia, Philadelphia, PA 19104, USA.
| |
Collapse
|
8
|
Kiguchi Y, Nishijima S, Kumar N, Hattori M, Suda W. Long-read metagenomics of multiple displacement amplified DNA of low-biomass human gut phageomes by SACRA pre-processing chimeric reads. DNA Res 2021; 28:6377780. [PMID: 34586399 DOI: 10.1093/dnares/dsab019] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2020] [Indexed: 01/21/2023] Open
Abstract
The human gut bacteriophage community (phageome) plays an important role in the host's health and disease; however, the entire structure is poorly understood, partly owing to the generation of many incomplete genomes in conventional short-read metagenomics. Here, we show long-read metagenomics of amplified DNA of low-biomass phageomes with multiple displacement amplification (MDA), involving the development of a novel bioinformatics tool, split amplified chimeric read algorithm (SACRA), that efficiently pre-processed numerous chimeric reads generated through MDA. Using five samples, SACRA markedly reduced the average chimera ratio from 72% to 1.5% in PacBio reads with an average length of 1.8 kb. De novo assembly of chimera-less PacBio long reads reconstructed contigs of ≥5 kb with an average proportion of 27%, which was 1% in contigs from MiSeq short reads, thereby dramatically improving contig length and genome completeness. Comparison of PacBio and MiSeq contigs found MiSeq contig fragmentations frequently near local repeats and hypervariable regions in the phage genomes, and those caused by multiple homologous phage genomes coexisting in the community. We also developed a reference-independent method to assess the completeness of the linear phage genomes. Overall, we established a SACRA-coupled long-read metagenomics robust to highly diverse gut phageomes, identifying high-quality circular and linear phage genomes with adequate sequence quantity.
Collapse
Affiliation(s)
- Yuya Kiguchi
- Cooperative Major in Advanced Health Science, Graduate School of Advanced Science and Engineering, Waseda University, Tokyo 169-8555, Japan
- Computational Bio Big-Data Open Innovation Laboratory (CBBD-OIL), National Institute of Advanced Industrial Science and Technology, Tokyo 169-8555, Japan
- Laboratory for Microbiome Sciences, RIKEN Center for Integrative Medical Sciences, Yokohama 230-0045, Japan
| | - Suguru Nishijima
- Computational Bio Big-Data Open Innovation Laboratory (CBBD-OIL), National Institute of Advanced Industrial Science and Technology, Tokyo 169-8555, Japan
- Integrated Institute for Regulatory Science, Waseda University, Tokyo 169-8555, Japan
- Structural and Computational Biology Unit, European Molecular Biology Laboratory, Heidelberg, Germany
| | - Naveen Kumar
- Laboratory for Microbiome Sciences, RIKEN Center for Integrative Medical Sciences, Yokohama 230-0045, Japan
| | - Masahira Hattori
- Cooperative Major in Advanced Health Science, Graduate School of Advanced Science and Engineering, Waseda University, Tokyo 169-8555, Japan
- Laboratory for Microbiome Sciences, RIKEN Center for Integrative Medical Sciences, Yokohama 230-0045, Japan
| | - Wataru Suda
- Laboratory for Microbiome Sciences, RIKEN Center for Integrative Medical Sciences, Yokohama 230-0045, Japan
| |
Collapse
|
9
|
Stephens Z, Milosevic D, Kipp B, Grebe S, Iyer RK, Kocher JPA. PB-Motif-A Method for Identifying Gene/Pseudogene Rearrangements With Long Reads: An Application to CYP21A2 Genotyping. Front Genet 2021; 12:716586. [PMID: 34394200 PMCID: PMC8355628 DOI: 10.3389/fgene.2021.716586] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2021] [Accepted: 07/05/2021] [Indexed: 12/30/2022] Open
Abstract
Long read sequencing technologies have the potential to accurately detect and phase variation in genomic regions that are difficult to fully characterize with conventional short read methods. These difficult to sequence regions include several clinically relevant genes with highly homologous pseudogenes, many of which are prone to gene conversions or other types of complex structural rearrangements. We present PB-Motif, a new method for identifying rearrangements between two highly homologous genomic regions using PacBio long reads. PB-Motif leverages clustering and filtering techniques to efficiently report rearrangements in the presence of sequencing errors and other systematic artifacts. Supporting reads for each high-confidence rearrangement can then be used for copy number estimation and phased variant calling. First, we demonstrate PB-Motif's accuracy with simulated sequence rearrangements of PMS2 and its pseudogene PMS2CL using simulated reads sweeping over a range of sequencing error rates. We then apply PB-Motif to 26 clinical samples, characterizing CYP21A2 and its pseudogene CYP21A1P as part of a diagnostic assay for congenital adrenal hyperplasia. We successfully identify damaging variation and patient carrier status concordant with clinical diagnosis obtained from multiplex ligation-dependent amplification (MLPA) and Sanger sequencing. The source code is available at: github.com/zstephens/pb-motif.
Collapse
Affiliation(s)
- Zachary Stephens
- Department of Electrical and Computer Engineering, University of Illinois Urbana-Champaign, Urbana, IL, United States
| | | | | | | | - Ravishankar K Iyer
- Department of Electrical and Computer Engineering, University of Illinois Urbana-Champaign, Urbana, IL, United States
| | | |
Collapse
|
10
|
Modlin SJ, Robinhold C, Morrissey C, Mitchell SN, Ramirez-Busby SM, Shmaya T, Valafar F. Exact mapping of Illumina blind spots in the Mycobacterium tuberculosis genome reveals platform-wide and workflow-specific biases. Microb Genom 2021; 7. [PMID: 33502304 PMCID: PMC8190613 DOI: 10.1099/mgen.0.000465] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Whole-genome sequencing (WGS) is fundamental to Mycobacterium tuberculosis basic research and many clinical applications. Coverage across Illumina-sequenced M. tuberculosis genomes is known to vary with sequence context, but this bias is poorly characterized. Here, through a novel application of phylogenomics that distinguishes genuine coverage bias from deletions, we discern Illumina ‘blind spots’ in the M. tuberculosis reference genome for seven sequencing workflows. We find blind spots to be widespread, affecting 529 genes, and provide their exact coordinates, enabling salvage of unaffected regions. Fifty-seven pe/ppe genes (the primary families assumed to exhibit Illumina bias) lack blind spots entirely, while the remaining pe/ppe genes account for 55.1 % of blind spots. Surprisingly, we find coverage bias persists in homopolymers as short as 6 bp, shorter tracts than previously reported. While G+C-rich regions challenge all Illumina sequencing workflows, a modified Nextera library preparation that amplifies DNA with a high-fidelity polymerase markedly attenuates coverage bias in G+C-rich and homopolymeric sequences, expanding the ‘Illumina-sequenceable’ genome. Through these findings, and by defining workflow-specific exclusion criteria, we spotlight effective strategies for handling bias in M. tuberculosis Illumina WGS. This empirical analysis framework may be used to systematically evaluate coverage bias in other species using existing sequencing data.
Collapse
Affiliation(s)
- Samuel J Modlin
- Laboratory for Pathogenesis of Clinical Drug Resistance and Persistence, School of Public Health, San Diego State University, San Diego, CA 92182, USA
| | - Cassidy Robinhold
- Laboratory for Pathogenesis of Clinical Drug Resistance and Persistence, School of Public Health, San Diego State University, San Diego, CA 92182, USA
| | - Christopher Morrissey
- Laboratory for Pathogenesis of Clinical Drug Resistance and Persistence, School of Public Health, San Diego State University, San Diego, CA 92182, USA
| | - Scott N Mitchell
- Laboratory for Pathogenesis of Clinical Drug Resistance and Persistence, School of Public Health, San Diego State University, San Diego, CA 92182, USA
| | - Sarah M Ramirez-Busby
- Laboratory for Pathogenesis of Clinical Drug Resistance and Persistence, School of Public Health, San Diego State University, San Diego, CA 92182, USA
| | - Tal Shmaya
- Laboratory for Pathogenesis of Clinical Drug Resistance and Persistence, School of Public Health, San Diego State University, San Diego, CA 92182, USA
| | - Faramarz Valafar
- Laboratory for Pathogenesis of Clinical Drug Resistance and Persistence, School of Public Health, San Diego State University, San Diego, CA 92182, USA
| |
Collapse
|
11
|
Abstract
The male-specific Y chromosome harbors genes important for sperm production. Because Y is repetitive, its DNA sequence was deciphered for only a few species, and its evolution remains elusive. Here we compared the Y chromosomes of great apes (human, chimpanzee, bonobo, gorilla, and orangutan) and found that many of their repetitive sequences and multicopy genes were likely already present in their common ancestor. Y repeats had increased intrachromosomal contacts, which might facilitate preservation of genes and gene regulatory elements. Chimpanzee and bonobo, experiencing high sperm competition, underwent many DNA changes and gene losses on the Y. Our research is significant for understanding the role of the Y chromosome in reproduction of nonhuman great apes, all of which are endangered. The mammalian male-specific Y chromosome plays a critical role in sex determination and male fertility. However, because of its repetitive and haploid nature, it is frequently absent from genome assemblies and remains enigmatic. The Y chromosomes of great apes represent a particular puzzle: their gene content is more similar between human and gorilla than between human and chimpanzee, even though human and chimpanzee share a more recent common ancestor. To solve this puzzle, here we constructed a dataset including Ys from all extant great ape genera. We generated assemblies of bonobo and orangutan Ys from short and long sequencing reads and aligned them with the publicly available human, chimpanzee, and gorilla Y assemblies. Analyzing this dataset, we found that the genus Pan, which includes chimpanzee and bonobo, experienced accelerated substitution rates. Pan also exhibited elevated gene death rates. These observations are consistent with high levels of sperm competition in Pan. Furthermore, we inferred that the great ape common ancestor already possessed multicopy sequences homologous to most human and chimpanzee palindromes. Nonetheless, each species also acquired distinct ampliconic sequences. We also detected increased chromatin contacts between and within palindromes (from Hi-C data), likely facilitating gene conversion and structural rearrangements. Our results highlight the dynamic mode of Y chromosome evolution and open avenues for studies of male-specific dispersal in endangered great ape species.
Collapse
|
12
|
Scheunert A, Dorfner M, Lingl T, Oberprieler C. Can we use it? On the utility of de novo and reference-based assembly of Nanopore data for plant plastome sequencing. PLoS One 2020; 15:e0226234. [PMID: 32208422 PMCID: PMC7092973 DOI: 10.1371/journal.pone.0226234] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2019] [Accepted: 02/28/2020] [Indexed: 12/13/2022] Open
Abstract
The chloroplast genome harbors plenty of valuable information for phylogenetic research. Illumina short-read data is generally used for de novo assembly of whole plastomes. PacBio or Oxford Nanopore long reads are additionally employed in hybrid approaches to enable assembly across the highly similar inverted repeats of a chloroplast genome. Unlike for PacBio, plastome assemblies based solely on Nanopore reads are rarely found, due to their high error rate and non-random error profile. However, the actual quality decline connected to their use has rarely been quantified. Furthermore, no study has employed reference-based assembly using Nanopore reads, which is common with Illumina data. Using Leucanthemum Mill. as an example, we compared the sequence quality of seven chloroplast genome assemblies of the same species, using combinations of two sequencing platforms and three analysis pipelines. In addition, we assessed the factors which might influence Nanopore assembly quality during sequence generation and bioinformatic processing. The consensus sequence derived from de novo assembly of Nanopore data had a sequence identity of 99.59% compared to Illumina short-read de novo assembly. Most of the errors detected were indels (81.5%), and a large majority of them is part of homopolymer regions. The quality of reference-based assembly is heavily dependent upon the choice of a close-enough reference. When using a reference with 0.83% sequence divergence from the studied species, mapping of Nanopore reads results in a consensus comparable to that from Nanopore de novo assembly, and of only slightly inferior quality compared to a reference-based assembly with Illumina data. For optimal de novo assembly of Nanopore data, appropriate filtering of contaminants and chimeric sequences, as well as employing moderate read coverage, is essential. Based on these results, we conclude that Nanopore long reads are a suitable alternative to Illumina short reads in plastome phylogenomics. Few errors remain in the finalized assembly, which can be easily masked in phylogenetic analyses without loss in analytical accuracy. The easily applicable and cost-effective technology might warrant more attention by researchers dealing with plant chloroplast genomes.
Collapse
Affiliation(s)
- Agnes Scheunert
- Evolutionary and Systematic Botany Group, Institute of Plant Sciences, University of Regensburg, Regensburg, Germany
| | - Marco Dorfner
- Evolutionary and Systematic Botany Group, Institute of Plant Sciences, University of Regensburg, Regensburg, Germany
| | - Thomas Lingl
- Evolutionary and Systematic Botany Group, Institute of Plant Sciences, University of Regensburg, Regensburg, Germany
| | - Christoph Oberprieler
- Evolutionary and Systematic Botany Group, Institute of Plant Sciences, University of Regensburg, Regensburg, Germany
| |
Collapse
|
13
|
Complete Sequence and Annotation of the Mycoplasma phocicerebrale Strain 1049 T Genome. Microbiol Resour Announc 2019; 8:8/25/e00514-19. [PMID: 31221652 PMCID: PMC6588373 DOI: 10.1128/mra.00514-19] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
The Mycoplasma phocicerebrale genome was analyzed to better understand this opportunistic pathogen. Amplification with ϕ29 polymerase was used to generate enough genomic DNA for large-insert library construction. Like other mycoplasmas from seals, M. phocicerebrale encodes an immunosuppressor that may predispose susceptibility to infection or influence intercurrent diseases of affected hosts. The Mycoplasma phocicerebrale genome was analyzed to better understand this opportunistic pathogen. Amplification with ϕ29 polymerase was used to generate enough genomic DNA for large-insert library construction. Like other mycoplasmas from seals, M. phocicerebrale encodes an immunosuppressor that may predispose susceptibility to infection or influence intercurrent diseases of affected hosts.
Collapse
|