Long-read human genome sequencing and its applications

Straglr: discovering and genotyping tandem repeat expansions using whole genome long-read sequences

Tandem repeat (TR) expansion is the underlying cause of over 40 neurological disorders. Long-read sequencing offers an exciting avenue over conventional technologies for detecting TR expansions. Here, we present Straglr, a robust software tool for both targeted genotyping and novel expansion detection from long-read alignments. We benchmark Straglr using various simulations, targeted genotyping data of cell lines carrying expansions of known diseases, and whole genome sequencing data with chromosome-scale assembly. Our results suggest that Straglr may be useful for investigating disease-associated TR expansions using long-read sequencing.

Perspectives and Benefits of High-Throughput Long-Read Sequencing in Microbial Ecology

Short-read, high-throughput sequencing (HTS) methods have yielded numerous important insights into microbial ecology and function. Yet, in many instances short-read HTS techniques are suboptimal, for example, by providing insufficient phylogenetic resolution or low integrity of assembled genomes. Single-molecule and synthetic long-read (SLR) HTS methods have successfully ameliorated these limitations. In addition, nanopore sequencing has generated a number of unique analysis opportunities, such as rapid molecular diagnostics and direct RNA sequencing, and both Pacific Biosciences (PacBio) and nanopore sequencing support detection of epigenetic modifications. Although initially suffering from relatively low sequence quality, recent advances have greatly improved the accuracy of long-read sequencing technologies. In spite of great technological progress in recent years, the long-read HTS methods (PacBio and nanopore sequencing) are still relatively costly, require large amounts of high-quality starting material, and commonly need specific solutions in various analysis steps. Despite these challenges, long-read sequencing technologies offer high-quality, cutting-edge alternatives for testing hypotheses about microbiome structure and functioning as well as assembly of eukaryote genomes from complex environmental DNA samples.

Enzymatic approaches for profiling cytosine methylation and hydroxymethylation

Background

In mammals, modifications to cytosine bases, particularly in cytosine-guanine (CpG) dinucleotide contexts, play a major role in shaping the epigenome. The canonical epigenetic mark is 5-methylcytosine (5mC), but oxidized versions of 5mC, including 5-hydroxymethylcytosine (5hmC), are now known to be important players in epigenomic dynamics. Understanding the functional role of these modifications in gene regulation, normal development, and pathological conditions requires the ability to localize these modifications in genomic DNA. The classical approach for sequencing cytosine modifications has involved differential deamination via the chemical sodium bisulfite; however, bisulfite is destructive, limiting its utility in important biological or clinical settings where detection of low frequency populations is critical. Additionally, bisulfite fails to resolve 5mC from 5hmC.

Scope of review

To summarize how enzymatic rather than chemical approaches can be leveraged to localize and resolve different cytosine modifications in a non-destructive manner.

Major conclusions

Nature offers a suite of enzymes with biological roles in cytosine modification in organisms spanning from bacteriophages to mammals. These enzymatic activities include methylation by DNA methyltransferases, oxidation of 5mC by TET family enzymes, hypermodification of 5hmC by glucosyltransferases, and the generation of transition mutations from cytosine to uracil by DNA deaminases. Here, we describe how insights into the natural reactivities of these DNA-modifying enzymes can be leveraged to convert them into powerful biotechnological tools. Application of these enzymes in sequencing can be accomplished by relying on their natural activity, exploiting their ability to discriminate between cytosine modification states, reacting them with functionalized substrate analogs to introduce chemical handles, or engineering the DNA-modifying enzymes to take on new reactivities. We describe how these enzymatic reactions have been combined and permuted to localize DNA modifications with high specificity and without the destructive limitations posed by chemical methods for epigenetic sequencing.

•

Chemical sequencing methods damage DNA and can confound cytosine modifications.

•

DNA modifying enzymes offer non-destructive and selective biotechnological tools.

•

DNA deaminases, methyltransferases, oxygenases and glucosyltransferases can be used.

•

Permuting enzymes with various activities can reveal distinct cytosine states.

•

Engineered enzymes utilizing unnatural co-substrates expand sequencing scope.

The BS variant of C4 protects against age-related loss of white matter microstructural integrity

Age-related loss of white matter microstructural integrity is a major determinant of cognitive decline, dementia, and gait disorders. However, the mechanisms and molecular pathways that contribute to this loss of integrity remain elusive. We performed a GWAS of white matter microstructural integrity as quantified by diffusion MRI metrics (mean diffusivity, MD; and fractional anisotropy, FA) in up to 31,128 individuals from UK Biobank (age 45-81 years) based on a 2 degrees of freedom (2df) test of single nucleotide polymorphism (SNP) and SNP x age effects. We identified 18 loci that were associated at genome-wide significance with either MD (N = 16) or FA (N = 6). Among the top loci was a region on chromosome 6 encoding the human major histocompatibility complex (MHC). Variants in the MHC region were strongly associated with both MD (best SNP: 6:28866209_TTTTG_T, beta(SE)=-0.069(0.009); 2df p = 6.5x10-15) and FA (best SNP: rs3129787, beta(SE)=-0.056(0.008); 2df p = 3.5x10-12). Of the imputed HLA alleles and complement component 4 (C4) structural haplotype variants in the human MHC, the strongest association was with the C4-BS variant (for MD: beta(SE)=-0.070(0.010); p = 2.7x10-11; for FA: beta(SE)=-0.054(0.011); p = 1.6x10-7). After conditioning on C4-BS no associations with HLA alleles remained significant. The protective influence of C4-BS was stronger in older subjects (age ≥ 65; interaction p = 0.0019 (MD), p = 0.015 (FA)) and in subjects without a history of smoking (interaction p = 0.00093 (MD), p = 0.021 (FA)). Taken together, our findings demonstrate a role of the complement system and of gene-environment interactions in age-related loss of white matter microstructural integrity.

Targeted long-read sequencing identifies missing disease-causing variation

Despite widespread clinical genetic testing, many individuals with suspected genetic conditions lack a precise diagnosis, limiting their opportunity to take advantage of state-of-the-art treatments. In some cases, testing reveals difficult-to-evaluate structural differences, candidate variants that do not fully explain the phenotype, single pathogenic variants in recessive disorders, or no variants in genes of interest. Thus, there is a need for better tools to identify a precise genetic diagnosis in individuals when conventional testing approaches have been exhausted. We performed targeted long-read sequencing (T-LRS) using adaptive sampling on the Oxford Nanopore platform on 40 individuals, 10 of whom lacked a complete molecular diagnosis. We computationally targeted up to 151 Mbp of sequence per individual and searched for pathogenic substitutions, structural variants, and methylation differences using a single data source. We detected all genomic aberrations-including single-nucleotide variants, copy number changes, repeat expansions, and methylation differences-identified by prior clinical testing. In 8/8 individuals with complex structural rearrangements, T-LRS enabled more precise resolution of the mutation, leading to changes in clinical management in one case. In ten individuals with suspected Mendelian conditions lacking a precise genetic diagnosis, T-LRS identified pathogenic or likely pathogenic variants in six and variants of uncertain significance in two others. T-LRS accurately identifies pathogenic structural variants, resolves complex rearrangements, and identifies Mendelian variants not detected by other technologies. T-LRS represents an efficient and cost-effective strategy to evaluate high-priority genes and regions or complex clinical testing results.

Transcription Factor-Based Genetic Engineering in Microalgae

Sequence-specific DNA-binding transcription factors (TFs) are key components of gene regulatory networks. Advances in high-throughput sequencing have facilitated the rapid acquisition of whole genome assembly and TF repertoires in microalgal species. In this review, we summarize recent advances in gene discovery and functional analyses, especially for transcription factors in microalgal species. Specifically, we provide examples of the genome-scale identification of transcription factors in genome-sequenced microalgal species and showcase their application in the discovery of regulators involved in various cellular functions. Herein, we highlight TF-based genetic engineering as a promising framework for designing microalgal strains for microalgal-based bioproduction.

Association of assisted reproductive technology, germline de novo mutations and congenital heart defects in a prospective birth cohort study

Emerging evidence suggests that children conceived through assisted reproductive technology (ART) have a higher risk of congenital heart defects (CHDs) even when there is no family history. De novo mutation (DNM) is a well-known cause of sporadic congenital diseases; however, whether ART procedures increase the number of germline DNM (gDNM) has not yet been well studied. Here, we performed whole-genome sequencing of 1137 individuals from 160 families conceived through ART and 205 families conceived spontaneously. Children conceived via ART carried 4.59 more gDNMs than children conceived spontaneously, including 3.32 paternal and 1.26 maternal DNMs, after correcting for parental age at conception, cigarette smoking, alcohol drinking, and exercise behaviors. Paternal DNMs in offspring conceived via ART are characterized by C>T substitutions at CpG sites, which potentially affect protein-coding genes and are significantly associated with the increased risk of CHD. In addition, the accumulation of non-coding functional mutations was independently associated with CHD and 87.9% of the mutations were originated from the father. Among ART offspring, infertility of the father was associated with elevated paternal DNMs; usage of both recombinant and urinary follicle-stimulating hormone and high-dosage human chorionic gonadotropin trigger was associated with an increase of maternal DNMs. In sum, the increased gDNMs in offspring conceived by ART were primarily originated from fathers, indicating that ART itself may not be a major reason for the accumulation of gDNMs. Our findings emphasize the importance of evaluating the germline status of the fathers in families with the use of ART.

Characterizing nucleotide variation and expansion dynamics in human-specific variable number tandem repeats

There are more than 55,000 variable number tandem repeats (VNTRs) in the human genome, notable for both their striking polymorphism and mutability. Despite their role in human evolution and genomic variation, they have yet to be studied collectively and in detail, partially owing to their large size, variability, and predominant location in noncoding regions. Here, we examine 467 VNTRs that are human-specific expansions, unique to one location in the genome, and not associated with retrotransposons. We leverage publicly available long-read genomes, including from the Human Genome Structural Variant Consortium, to ascertain the exact nucleotide composition of these VNTRs and compare their composition of alleles. We then confirm repeat unit composition in more than 3000 short-read samples from the 1000 Genomes Project. Our analysis reveals that these VNTRs contain highly structured repeat motif organization, modified by frequent deletion and duplication events. Although overall VNTR compositions tend to remain similar between 1000 Genomes Project superpopulations, we describe a notable exception with substantial differences in repeat composition (in PCBP3), as well as several VNTRs that are significantly different in length between superpopulations (in ART1, PROP1, DYNC2I1, and LOC102723906). We also observe that most of these VNTRs are expanded in archaic human genomes, yet remain stable in length between single generations. Collectively, our findings indicate that repeat motif variability, repeat composition, and repeat length are all informative modalities to consider when characterizing VNTRs and their contribution to genomic variation.

Application of long-read sequencing to the detection of structural variants in human cancer genomes

In recent years, the so-called long-read sequencing technology has had a substantial impact on various aspects of genome sciences. Here, we introduce recent studies of cancerous structural variants (SVs) using long-read sequencing technologies, namely Pacific Biosciences (PacBio) sequencers, Oxford Nanopore Technologies (ONT) sequencers, and linked-read methods. By taking advantage of long-read lengths, these technologies have enabled the precise detection of SVs, including long insertions by transposable elements, such as LINE-1. In addition to SV detection, the epigenome status (including DNA methylation and haplotype information) surrounding SV loci has also been unveiled by long-read sequencing technologies, to identify the effects of SVs. Among the various research fields in which long-read sequencing has been applied, cancer genomics has shown the most remarkable advances. In fact, many studies are beginning to shed light on the detection of SVs and the elucidation of their complex structures in various types of cancer. In the particular case of cancers, we summarize the technical limitations of the application of this technology to the analysis of clinical samples. We will introduce recent achievements from this viewpoint. However, a similar approach will be started for other applications in the near future. Therefore, by complementing the current short-read sequencing analysis, long-read sequencing should reveal the complex nature of human genomes in their healthy and disease states, which will open a new opportunity for a better understanding of disease development and for a novel strategy for drug development.

Applications of Raman Spectroscopy in Bacterial Infections: Principles, Advantages, and Shortcomings

Infectious diseases caused by bacterial pathogens are important public issues. In addition, due to the overuse of antibiotics, many multidrug-resistant bacterial pathogens have been widely encountered in clinical settings. Thus, the fast identification of bacteria pathogens and profiling of antibiotic resistance could greatly facilitate the precise treatment strategy of infectious diseases. So far, many conventional and molecular methods, both manual or automatized, have been developed for in vitro diagnostics, which have been proven to be accurate, reliable, and time efficient. Although Raman spectroscopy (RS) is an established technique in various fields such as geochemistry and material science, it is still considered as an emerging tool in research and diagnosis of infectious diseases. Based on current studies, it is too early to claim that RS may provide practical guidelines for microbiologists and clinicians because there is still a gap between basic research and clinical implementation. However, due to the promising prospects of label-free detection and noninvasive identification of bacterial infections and antibiotic resistance in several single steps, it is necessary to have an overview of the technique in terms of its strong points and shortcomings. Thus, in this review, we went through recent studies of RS in the field of infectious diseases, highlighting the application potentials of the technique and also current challenges that prevent its real-world applications.

Complete Genome Assembly of Myxococcus xanthus Strain DZ2 Using Long High-Fidelity (HiFi) Reads Generated with PacBio Technology

Myxococcus xanthus is a Gram-negative social bacterium belonging to the order Myxococcales of the class Deltaproteobacteria. It is a facultative social predator found in soils across the globe and is thought to be crucial for the microbial ecosystem. Here, we report a complete high-quality reference genome of the M. xanthus strain DZ2.

High-molecular weight DNA extraction, clean-up and size selection for long-read sequencing

Rapid advancements in long-read sequencing technologies have transformed read lengths from bps to Mbps, which has enabled chromosome-scale genome assemblies. However, read lengths are now becoming limited by the extraction of pure high-molecular weight DNA suitable for long-read sequencing, which is particularly challenging in plants and fungi. To overcome this, we present a protocol collection; high-molecular weight DNA extraction, clean-up and size selection for long-read sequencing. We optimised a gentle magnetic bead based high-molecular weight DNA extraction, which is presented here in detail. The protocol circumvents spin columns and high-centrifugation, to limit DNA fragmentation. The protocol is scalable based on tissue input, which can be used on many species of plants, fungi, reptiles and bacteria. It is also cost effective compared to kit-based protocols and hence applicable at scale in low resource settings. An optional sorbitol wash is listed and is highly recommended for plant and fungal tissues. To further remove any remaining contaminants such as phenols and polysaccharides, optional DNA clean-up and size selection strategies are given. This protocol collection is suitable for all common long-read sequencing platforms, such as technologies offered by PacBio and Oxford Nanopore. Using these protocols, sequencing on the Oxford Nanopore MinION can achieve read length N50 values of 30-50 kb, with reads exceeding 200 kb and outputs ranging from 15-30 Gbp. This has been routinely achieved with various plant, fungi, animal and bacteria samples.

A Curriculum for Genomic Education of Molecular Genetic Pathology Fellows: A Report of the Association for Molecular Pathology Training and Education Committee

Molecular genetic pathology (MGP) is a subspecialty of pathology and medical genetics and genomics. Genomic testing, which we define as that which generates large data sets and interrogates large segments of the genome in a single assay, is increasingly recognized as essential for optimal patient care through precision medicine. The most common genomic testing technologies in clinical laboratories are next-generation sequencing and microarray. It is essential to train in these methods and to consider the data generated in the context of the diagnosis, medical history, and other clinical findings of individual patients. Accordingly, updating the MGP fellowship curriculum to include genomics is timely, important, and challenging. At the completion of training, an MGP fellow should be capable of independently interpreting and signing out results of a wide range of genomic assays and, given the appropriate context and institutional support, of developing and validating new assays in compliance with applicable regulations. The Genomics Task Force of the MGP Program Directors, a working group of the Association for Molecular Pathology Training and Education Committee, has developed a genomics curriculum framework and recommendations specific to the MGP fellowship. These recommendations are presented for consideration and implementation by MGP fellowship programs with the understanding that MGP programs exist in a diversity of clinical practice environments with a spectrum of available resources.

A simple approach for effective shearing and reliable concentration measurement of ultra-high-molecular-weight DNA

Pipetting and concentration measurement of viscous ultra-high-molecular-weight (UHMW) DNA samples is challenging and often highly imprecise. Effective guidelines for handling UHMW samples are missing in the field. Herein, a simple and low-cost workflow is presented that enables accurate pipetting and reliable concentration measurement. Central to the workflow is the shearing of representative small aliquots of UHMW DNA samples to a fragment size <150 kb by vortexing them for 1 min with a glass bead in a round-bottomed 2-ml tube. Additionally, a solution is provided for accurate quantitation of high-molecular-weight DNA with fluorometric (Qubit [Thermo Fisher Scientific, MA, USA]) methods by using an appropriate genomic DNA standard, resulting in values that match spectrophotometric (Nanodrop [Thermo Fisher Scientific]) optical density readings.

Brain Regional Differences in Hexanucleotide Repeat Length in X-Linked Dystonia-Parkinsonism Using Nanopore Sequencing

Objective

Our study investigated the presence of regional differences in hexanucleotide repeat number in postmortem brain tissues of 2 patients with X-linked dystonia-parkinsonism (XDP), a combined dystonia-parkinsonism syndrome modified by a (CCCTCT)_n repeat within the causal SINE-VNTR-Alu retrotransposon insertion in the TAF1 gene.

Methods

Genomic DNA was extracted from blood and postmortem brain samples, including the basal ganglia and cortex from both patients and from the cerebellum, midbrain, and pituitary gland from 1 patient. Repeat sizing was performed using fragment analysis, small-pool PCR-based Southern blotting, and Oxford nanopore sequencing.

Results

The basal ganglia (p < 0.001) and cerebellum (p < 0.001) showed higher median repeat numbers and higher degrees of repeat instability compared with blood.

Conclusions

Somatic repeat instability may predominate in brain regions selectively affected in XDP, thereby hinting at its potential role in disease manifestation and modification.

Building pan-genome infrastructures for crop plants and their use in association genetics

Pan-genomic studies aim at representing the entire sequence diversity within a species to provide useful resources for evolutionary studies, functional genomics and breeding of cultivated plants. Cost reductions in high-throughput sequencing and advances in sequence assembly algorithms have made it possible to create multiple reference genomes along with a catalogue of all forms of genetic variations in plant species with large and complex or polyploid genomes. In this review, we summarize the current approaches to building pan-genomes as an in silico representation of plant sequence diversity and outline relevant methods for their effective utilization in linking structural with phenotypic variation. We propose as future research avenues (i) transcriptomic and epigenomic studies across multiple reference genomes and (ii) the development of user-friendly and feature-rich pan-genome browsers.

Exome variant discrepancies due to reference-genome differences

Despite release of the GRCh38 human reference genome more than seven years ago, GRCh37 remains more widely used by most research and clinical laboratories. To date, no study has quantified the impact of utilizing different reference assemblies for the identification of variants associated with rare and common diseases from large-scale exome-sequencing data. By calling variants on both the GRCh37 and GRCh38 references, we identified single-nucleotide variants (SNVs) and insertion-deletions (indels) in 1,572 exomes from participants with Mendelian diseases and their family members. We found that a total of 1.5% of SNVs and 2.0% of indels were discordant when different references were used. Notably, 76.6% of the discordant variants were clustered within discrete discordant reference patches (DISCREPs) comprising only 0.9% of loci targeted by exome sequencing. These DISCREPs were enriched for genomic elements including segmental duplications, fix patch sequences, and loci known to contain alternate haplotypes. We identified 206 genes significantly enriched for discordant variants, most of which were in DISCREPs and caused by multi-mapped reads on the reference assembly that lacked the variant call. Among these 206 genes, eight are implicated in known Mendelian diseases and 53 are associated with common phenotypes from genome-wide association studies. In addition, variant interpretations could also be influenced by the reference after lifting-over variant loci to another assembly. Overall, we identified genes and genomic loci affected by reference assembly choice, including genes associated with Mendelian disorders and complex human diseases that require careful evaluation in both research and clinical applications.

Nanopore Technology and Its Applications in Gene Sequencing

In recent years, nanopore technology has become increasingly important in the field of life science and biomedical research. By embedding a nano-scale hole in a thin membrane and measuring the electrochemical signal, nanopore technology can be used to investigate the nucleic acids and other biomacromolecules. One of the most successful applications of nanopore technology, the Oxford Nanopore Technology, marks the beginning of the fourth generation of gene sequencing technology. In this review, the operational principle and the technology for signal processing of the nanopore gene sequencing are documented. Moreover, this review focuses on the applications using nanopore gene sequencing technology, including the diagnosis of cancer, detection of viruses and other microbes, and the assembly of genomes. These applications show that nanopore technology is promising in the field of biological and biomedical sensing.

Comparison of long-read sequencing technologies in interrogating bacteria and fly genomes

The newest generation of DNA sequencing technology is highlighted by the ability to generate sequence reads hundreds of kilobases in length. Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) have pioneered competitive long read platforms, with more recent work focused on improving sequencing throughput and per-base accuracy. We used whole-genome sequencing data produced by three PacBio protocols (Sequel II CLR, Sequel II HiFi, RS II) and two ONT protocols (Rapid Sequencing and Ligation Sequencing) to compare assemblies of the bacteria Escherichia coli and the fruit fly Drosophila ananassae. In both organisms tested, Sequel II assemblies had the highest consensus accuracy, even after accounting for differences in sequencing throughput. ONT and PacBio CLR had the longest reads sequenced compared to PacBio RS II and HiFi, and genome contiguity was highest when assembling these datasets. ONT Rapid Sequencing libraries had the fewest chimeric reads in addition to superior quantification of E. coli plasmids versus ligation-based libraries. The quality of assemblies can be enhanced by adopting hybrid approaches using Illumina libraries for bacterial genome assembly or polishing eukaryotic genome assemblies, and an ONT-Illumina hybrid approach would be more cost-effective for many users. Genome-wide DNA methylation could be detected using both technologies, however ONT libraries enabled the identification of a broader range of known E. coli methyltransferase recognition motifs in addition to undocumented D. ananassae motifs. The ideal choice of long read technology may depend on several factors including the question or hypothesis under examination. No single technology outperformed others in all metrics examined.

Comparative Analysis for the Performance of Long-Read-Based Structural Variation Detection Pipelines in Tandem Repeat Regions

There has been growing recognition of the vital links between structural variations (SVs) and diverse diseases. Research suggests that, with much longer DNA fragments and abundant contextual information, long-read technologies have advantages in SV detection even in complex repetitive regions. So far, several pipelines for calling SVs from long-read sequencing data have been proposed and used in human genome research. However, the performance of these pipelines is still lack of deep exploration and adequate comparison. In this study, we comprehensively evaluated the performance of three commonly used long-read SV detection pipelines, namely PBSV, Sniffles and PBHoney, especially the performance on detecting the SVs in tandem repeat regions (TRRs). Evaluated by using a robust benchmark for germline SV detection as the gold standard, we thoroughly estimated the precision, recall and F1 score of insertions and deletions detected by the pipelines. Our results revealed that all these pipelines clearly exhibited better performance outside TRRs than that in TRRs. The F1 scores of Sniffles in and outside TRRs were 0.60 and 0.76, respectively. The performance of PBSV was similar to that of Sniffles, and was generally higher than that of PBHoney. In conclusion, our findings can be benefit for choosing the appropriate pipelines in real practice and are good complementary to the application of long-read sequencing technologies in the research of rare diseases.

Comparative genome analysis using sample-specific string detection in accurate long reads

Motivation

Comparative genome analysis of two or more whole-genome sequenced (WGS) samples is at the core of most applications in genomics. These include the discovery of genomic differences segregating in populations, case-control analysis in common diseases and diagnosing rare disorders. With the current progress of accurate long-read sequencing technologies (e.g. circular consensus sequencing from PacBio sequencers), we can dive into studying repeat regions of the genome (e.g. segmental duplications) and hard-to-detect variants (e.g. complex structural variants).

Results

We propose a novel framework for comparative genome analysis through the discovery of strings that are specific to one genome ('samples-specific' strings). We have developed a novel, accurate and efficient computational method for the discovery of sample-specific strings between two groups of WGS samples. The proposed approach will give us the ability to perform comparative genome analysis without the need to map the reads and is not hindered by shortcomings of the reference genome and mapping algorithms. We show that the proposed approach is capable of accurately finding sample-specific strings representing nearly all variation (>98%) reported across pairs or trios of WGS samples using accurate long reads (e.g. PacBio HiFi data).

Availability and implementation

Data, code and instructions for reproducing the results presented in this manuscript are publicly available at https://github.com/Parsoa/PingPong.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Integrating Viral Metagenomics into an Ecological Framework

Viral metagenomics has expanded our knowledge of the ecology of uncultured viruses, within both environmental (e.g., terrestrial and aquatic) and host-associated (e.g., plants and animals, including humans) contexts. Here, we emphasize the implementation of an ecological framework in viral metagenomic studies to address questions in virology rarely considered ecological, which can change our perception of viruses and how they interact with their surroundings. An ecological framework explicitly considers diverse variants of viruses in populations that make up communities of interacting viruses, with ecosystem-level effects. It provides a structure for the study of the diversity, distributions, dynamics, and interactions of viruses with one another, hosts, and the ecosystem, including interactions with abiotic factors. An ecological framework in viral metagenomics stands poised to broadly expand our knowledge in basic and applied virology. We highlight specific fundamental research needs to capitalize on its potential and advance the field. Expected final online publication date for the Annual Review of Virology, Volume 8 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Long-read whole-genome methylation patterning using enzymatic base conversion and nanopore sequencing

Long-read whole-genome sequencing analysis of DNA methylation would provide useful information on the chromosomal context of gene expression regulation. Here we describe the development of a method that improves the read length generated by using the bisulfite-sequencing-based approach. In this method, we combined recently developed enzymatic base conversion, where an unmethylated cytosine (C) should be converted to thymine (T), with nanopore sequencing. After methylation-sensitive base conversion, the sequencing library was constructed using long-range polymerase chain reaction. This type of analysis is possible using a minimum of 1 ng genomic DNA, and an N50 read length of 3.4–7.6 kb is achieved. To analyze the produced data, which contained a substantial number of base mismatches due to sequence conversion and an inaccurate base read of the nanopore sequencing, a new analytical pipeline was constructed. To demonstrate the performance of long-read methylation sequencing, breast cancer cell lines and clinical specimens were subjected to analysis, which revealed the chromosomal methylation context of key cancer-related genes, allele-specific methylated genes, and repetitive or deletion regions. This method should convert the intractable specimens for which the amount of available genomic DNA is limited to the tractable targets.

Mining the gaps of chromosome 8

The first gapless, telomere-to-telomere sequence of a human autosome, chromosome 8, is complete. Sequencing and assembly of the corresponding centromere in the chimpanzee, orangutan and macaque reveals details of its rapid evolution over the past 25 million years.

Hi-C as a molecular rangefinder to examine genomic rearrangements

The mammalian genome is highly packed into the nucleus. Over the past decade, the development of Hi-C has contributed significantly to our understanding of the three-dimensional (3D) chromatin structure, uncovering the principles and functions of higher-order chromatin organizations. Recent studies have repositioned its property in spatial proximity measurement to address challenging problems in genome analyses including genome assembly, haplotype phasing, and the detection of genomic rearrangements. In particular, the power of Hi-C in detecting large-scale structural variations (SVs) in the cancer genome has been demonstrated, which is challenging to be addressed solely with short-read-based whole-genome sequencing analyses. In this review, we first provide a comprehensive view of Hi-C as an intuitive and effective SV detection tool. Then, we introduce recently developed bioinformatics tools utilizing Hi-C to investigate genomic rearrangements. Finally, we discuss the potential application of single-cell Hi-C to address the heterogeneity of genomic rearrangements and sub-population identification in the cancer genome.

Advancing crop genomics from lab to field

Crop genomics remains a key element in ensuring scientific progress to secure global food security. It has been two decades since the sequence of the first plant genome, that of Arabidopsis thaliana, was released, and soon after that the draft sequencing of the rice genome was completed. Since then, the genomes of more than 100 crops have been sequenced, plant genome research has expanded across multiple fronts and the next few years promise to bring further advances spurred by the advent of new technologies and approaches. We are likely to see continued innovations in crop genome sequencing, genetic mapping and the acquisition of multiple levels of biological data. There will be exciting opportunities to integrate genome-scale information across multiple scales of biological organization, leading to advances in our mechanistic understanding of crop biological processes, which will, in turn, provide greater impetus for translation of laboratory results to the field.

A high-quality bonobo genome refines the analysis of hominid evolution

The divergence of chimpanzee and bonobo provides one of the few examples of recent hominid speciation^1,2. Here we describe a fully annotated, high-quality bonobo genome assembly, which was constructed without guidance from reference genomes by applying a multiplatform genomics approach. We generate a bonobo genome assembly in which more than 98% of genes are completely annotated and 99% of the gaps are closed, including the resolution of about half of the segmental duplications and almost all of the full-length mobile elements. We compare the bonobo genome to those of other great apes^1,3–5 and identify more than 5,569 fixed structural variants that specifically distinguish the bonobo and chimpanzee lineages. We focus on genes that have been lost, changed in structure or expanded in the last few million years of bonobo evolution. We produce a high-resolution map of incomplete lineage sorting and estimate that around 5.1% of the human genome is genetically closer to chimpanzee or bonobo and that more than 36.5% of the genome shows incomplete lineage sorting if we consider a deeper phylogeny including gorilla and orangutan. We also show that 26% of the segments of incomplete lineage sorting between human and chimpanzee or human and bonobo are non-randomly distributed and that genes within these clustered segments show significant excess of amino acid replacement compared to the rest of the genome.

A high-quality bonobo genome assembly provides insights into incomplete lineage sorting in hominids and its relevance to gene evolution and the genetic relationship among living hominids.

Genomic Tackling of Human Satellite DNA: Breaking Barriers through Time

(Peri)centromeric repetitive sequences and, more specifically, satellite DNA (satDNA) sequences, constitute a major human genomic component. SatDNA sequences can vary on a large number of features, including nucleotide composition, complexity, and abundance. Several satDNA families have been identified and characterized in the human genome through time, albeit at different speeds. Human satDNA families present a high degree of sub-variability, leading to the definition of various subfamilies with different organization and clustered localization. Evolution of satDNA analysis has enabled the progressive characterization of satDNA features. Despite recent advances in the sequencing of centromeric arrays, comprehensive genomic studies to assess their variability are still required to provide accurate and proportional representation of satDNA (peri)centromeric/acrocentric short arm sequences. Approaches combining multiple techniques have been successfully applied and seem to be the path to follow for generating integrated knowledge in the promising field of human satDNA biology.

The Science and Art of Clinical Genetic Variant Classification and Its Impact on Test Accuracy

Clinical genetic variant classification science is a growing subspecialty of clinical genetics and genomics. The field's continued improvement is essential for the success of precision medicine in both germline (hereditary) and somatic (oncology) contexts. This review focuses on variant classification for DNA next-generation sequencing tests. We first summarize current limitations in variant discovery and definition, and then describe the current five- and four-tier classification systems outlined in dominant standards and guideline publications for germline and somatic tests, respectively. We then discuss measures of variant classification discordance and the field's bias for positive results, as well as considerations for panel size and population screening in the context of estimates of positive predictive value thatincorporate estimated variant classification imperfections. Finally, we share opinions on the current state of variant classification from some of the authors of the most widely used standards and guideline publications and from other domain experts.

Single-molecule optical genome mapping in nanochannels: multidisciplinarity at the nanoscale

The human genome contains multiple layers of information that extend beyond the genetic sequence. In fact, identical genetics do not necessarily yield identical phenotypes as evident for the case of two different cell types in the human body. The great variation in structure and function displayed by cells with identical genetic background is attributed to additional genomic information content. This includes large-scale genetic aberrations, as well as diverse epigenetic patterns that are crucial for regulating specific cell functions. These genetic and epigenetic patterns operate in concert in order to maintain specific cellular functions in health and disease. Single-molecule optical genome mapping is a high-throughput genome analysis method that is based on imaging long chromosomal fragments stretched in nanochannel arrays. The access to long DNA molecules coupled with fluorescent tagging of various genomic information presents a unique opportunity to study genetic and epigenetic patterns in the genome at a single-molecule level over large genomic distances. Optical mapping entwines synergistically chemical, physical, and computational advancements, to uncover invaluable biological insights, inaccessible by sequencing technologies. Here we describe the method's basic principles of operation, and review the various available mechanisms to fluorescently tag genomic information. We present some of the recent biological and clinical impact enabled by optical mapping and present recent approaches for increasing the method's resolution and accuracy. Finally, we discuss how multiple layers of genomic information may be mapped simultaneously on the same DNA molecule, thus paving the way for characterizing multiple genomic observables on individual DNA molecules.

Computational methods for chromosome-scale haplotype reconstruction

High-quality chromosome-scale haplotype sequences of diploid genomes, polyploid genomes, and metagenomes provide important insights into genetic variation associated with disease and biodiversity. However, whole-genome short read sequencing does not yield haplotype information spanning whole chromosomes directly. Computational assembly of shorter haplotype fragments is required for haplotype reconstruction, which can be challenging owing to limited fragment lengths and high haplotype and repeat variability across genomes. Recent advancements in long-read and chromosome-scale sequencing technologies, alongside computational innovations, are improving the reconstruction of haplotypes at the level of whole chromosomes. Here, we review recent and discuss methodological progress and perspectives in these areas.

Long-read genome sequencing for the molecular diagnosis of neurodevelopmental disorders

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. 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. We identified variants and created de novo assemblies in each trio, with global metrics indicating these datasets 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. The breadth and quality of variant detection, coupled to finding variants of clinical and research interest in two of six probands with unexplained NDDs, support the hypothesis that long-read genome sequencing can substantially improve rare disease genetic discovery rates.

Synteny-Based Genome Assembly for 16 Species of Heliconius Butterflies, and an Assessment of Structural Variation across the Genus

Heliconius butterflies (Lepidoptera: Nymphalidae) are a group of 48 neotropical species widely studied in evolutionary research. Despite the wealth of genomic data generated in past years, chromosomal level genome assemblies currently exist for only two species, Heliconius melpomene and Heliconius erato, each a representative of one of the two major clades of the genus. Here, we use these reference genomes to improve the contiguity of previously published draft genome assemblies of 16 Heliconius species. Using a reference-assisted scaffolding approach, we place and order the scaffolds of these genomes onto chromosomes, resulting in 95.7-99.9% of their genomes anchored to chromosomes. Genome sizes are somewhat variable among species (270-422 Mb) and in one small group of species (Heliconius hecale, Heliconius elevatus, and Heliconius pardalinus) expansions in genome size are driven mainly by repetitive sequences that map to four small regions in the H. melpomene reference genome. Genes from these repeat regions show an increase in exon copy number, an absence of internal stop codons, evidence of constraint on nonsynonymous changes, and increased expression, all of which suggest that at least some of the extra copies are functional. Finally, we conducted a systematic search for inversions and identified five moderately large inversions fixed between the two major Heliconius clades. We infer that one of these inversions was transferred by introgression between the lineages leading to the erato/sara and burneyi/doris clades. These reference-guided assemblies represent a major improvement in Heliconius genomic resources that enable further genetic and evolutionary discoveries in this genus.

Assessing the utility of long-read nanopore sequencing for rapid and efficient characterization of mobile element insertions

Short-read next generation sequencing (NGS) has become the predominant first-line technique used to diagnose patients with rare genetic conditions. Inherent limitations of short-read technology, notably for the detection and characterization of complex insertion-containing variants, are offset by the ability to concurrently screen many disease genes. "Third-generation" long-read sequencers are increasingly being deployed as an orthogonal adjunct technology, but their full potential for molecular genetic diagnosis has yet to be exploited. Here, we describe three diagnostic cases in which pathogenic mobile element insertions were refractory to characterization by short-read sequencing. To validate the accuracy of the long-read technology, we first used Sanger sequencing to confirm the integration sites and derive curated benchmark sequences of the variant-containing alleles. Long-read nanopore sequencing was then performed on locus-specific amplicons. Pairwise comparison between these data and the previously determined benchmark alleles revealed 100% identity of the variant-containing sequences. We demonstrate a number of technical advantages over existing wet-laboratory approaches, including in silico size selection of a mixed pool of amplification products, and the relative ease with which an automated informatics workflow can be established. Our findings add to a growing body of literature describing the diagnostic utility of long-read sequencing.

Opportunities and challenges for high-quality biodiversity tissue archives in the age of long-read sequencing

The technological ability to characterize genetic variation at a genome-wide scale provides an unprecedented opportunity to study the genetic underpinnings and evolutionary mechanisms that promote and sustain biodiversity. The transition from short- to long-read sequencing is particularly promising and allows a more holistic view on any changes in genetic diversity across time and space. Long-read sequencing has tremendous potential but sequencing success strongly depends on the long-range integrity of DNA molecules and therefore on the availability of high-quality tissue samples. With the scope of genomic experiments expanding and wild populations simultaneously disappearing at an unprecedented rate, access to high-quality samples may soon be a major concern for many projects. The need for high-quality biodiversity tissue archives is therefore urgent but sampling and preserving high-quality samples is not a trivial exercise. In this review, I will briefly outline how long-read sequencing can benefit the study of molecular ecology, how this will substantially increase the demand for high-quality tissues and why it is challenging to preserve DNA integrity. I will then provide an overview of preservation approaches and end with a call for support to acknowledge the efforts needed to assemble high-quality tissue archives. In doing so, I hope to simultaneously motivate field biologists to expand sampling practices and molecular biologists to develop (cost) efficient guidelines for the sampling and long-term storage of tissues. A concerted, interdisciplinary, effort is needed to catalogue the genetic variation underlying contemporary biodiversity and will eventually provide a critical resource for future studies.

Next-generation sequencing technologies: An overview

Since the days of Sanger sequencing, next-generation sequencing technologies have significantly evolved to provide increased data output, efficiencies, and applications. These next generations of technologies can be categorized based on read length. This review provides an overview of these technologies as two paradigms: short-read, or "second-generation," technologies, and long-read, or "third-generation," technologies. Herein, short-read sequencing approaches are represented by the most prevalent technologies, Illumina and Ion Torrent, and long-read sequencing approaches are represented by Pacific Biosciences and Oxford Nanopore technologies. All technologies are reviewed along with reported advantages and disadvantages. Until recently, short-read sequencing was thought to provide high accuracy limited by read-length, while long-read technologies afforded much longer read-lengths at the expense of accuracy. Emerging developments for third-generation technologies hold promise for the next wave of sequencing evolution, with the co-existence of longer read lengths and high accuracy.

Omics-based microbiome analysis in microbial ecology: from sequences to information

Microbial ecology is the study of microorganisms present in nature. It particularly focuses on microbial interactions with any biota and with surrounding environments. Microbial ecology is entering its golden age with innovative multi-omics methods triggered by next-generation sequencing technologies. However, the extraction of ecologically relevant information from ever-increasing omics data remains one of the most challenging tasks in microbial ecology. This special issue includes 11 review articles that provide an overview of the state of the art of omics-based approaches in the field of microbial ecology, with particular emphasis on the interpretation of omics data, environmental pollution tracking, interactions in microbiomes, and viral ecology.

Miniaturized DNA Sequencers for Personal Use: Unreachable Dreams or Achievable Goals

The appearance of next generation sequencing technology that features short read length with high measurement throughput and low cost has revolutionized the field of life science, medicine, and even computer science. The subsequent development of the third-generation sequencing technologies represented by nanopore and zero-mode waveguide techniques offers even higher speed and long read length with promising applications in portable and rapid genomic tests in field. Especially under the current circumstances, issues such as public health emergencies and global pandemics impose soaring demand on quick identification of origins and species of analytes through DNA sequences. In addition, future development of disease diagnosis, treatment, and tracking techniques may also require frequent DNA testing. As a result, DNA sequencers with miniaturized size and highly integrated components for personal and portable use to tackle increasing needs for disease prevention, personal medicine, and biohazard protection may become future trends. Just like many other biological and medical analytical systems that were originally bulky in sizes, collaborative work from various subjects in engineering and science eventually leads to the miniaturization of these systems. DNA sequencers that involve nanoprobes, detectors, microfluidics, microelectronics, and circuits as well as complex functional materials and structures are extremely complicated but may be miniaturized with technical advancement. This paper reviews the state-of-the-art technology in developing essential components in DNA sequencers and analyzes the feasibility to achieve miniaturized DNA sequencers for personal use. Future perspectives on the opportunities and associated challenges for compact DNA sequencers are also identified.

The extrachromosomal elements of the Naegleria genus: How little we know

There are currently 47 characterized species in the Naegleria genus of free-living amoebae. Each amoeba has thousands of extrachromosomal elements that are closed circular structures comprised of a single ribosomal DNA (rDNA) copy and a large non-rDNA sequence. Despite the presence of putative open reading frames and introns, ribosomal RNA is the only established transcript. A single origin of DNA replication (ori) has been mapped within the non-rDNA sequence for one species (N. gruberi), a finding that strongly indicates that these episomes replicate independently of the cell's chromosomal DNA component. This article reviews that which has been published about these interesting DNA elements and by analyzing available sequence data, discusses the possibility that different phylogenetically related clusters of Naegleria species individually conserve ori structures and suggests where the rRNA promoter and termination sites may be located.

long-read-tools.org: an interactive catalogue of analysis methods for long-read sequencing data

Background

The data produced by long-read third-generation sequencers have unique characteristics compared to short-read sequencing data, often requiring tailored analysis tools for tasks ranging from quality control to downstream processing. The rapid growth in software that addresses these challenges for different genomics applications is difficult to keep track of, which makes it hard for users to choose the most appropriate tool for their analysis goal and for developers to identify areas of need and existing solutions to benchmark against.

Findings

We describe the implementation of long-read-tools.org, an open-source database that organizes the rapidly expanding collection of long-read data analysis tools and allows its exploration through interactive browsing and filtering. The current database release contains 478 tools across 32 categories. Most tools are developed in Python, and the most frequent analysis tasks include base calling, de novo assembly, error correction, quality checking/filtering, and isoform detection, while long-read single-cell data analysis and transcriptomics are areas with the fewest tools available.

Conclusion

Continued growth in the application of long-read sequencing in genomics research positions the long-read-tools.org database as an essential resource that allows researchers to keep abreast of both established and emerging software to help guide the selection of the most relevant tool for their analysis needs.

Potential of whole-genome sequencing-based pharmacogenetic profiling

Pharmacogenetics represents a major driver of precision medicine, promising individualized drug selection and dosing. Traditionally, pharmacogenetic profiling has been performed using targeted genotyping that focuses on common/known variants. Recently, whole-genome sequencing (WGS) is emerging as a more comprehensive short-read next-generation sequencing approach, enabling both gene diagnostics and pharmacogenetic profiling, including rare/novel variants, in a single assay. Using the example of the pharmacogene CYP2D6, we demonstrate the potential of WGS-based pharmacogenetic profiling as well as emphasize the limitations of short-read next-generation sequencing. In the near future, we envision a shift toward long-read sequencing as the predominant method for gene diagnostics and pharmacogenetic profiling, providing unprecedented data quality and improving patient care.

Ratatosk: hybrid error correction of long reads enables accurate variant calling and assembly

A major challenge to long read sequencing data is their high error rate of up to 15%. We present Ratatosk, a method to correct long reads with short read data. We demonstrate on 5 human genome trios that Ratatosk reduces the error rate of long reads 6-fold on average with a median error rate as low as 0.22 %. SNP calls in Ratatosk corrected reads are nearly 99 % accurate and indel calls accuracy is increased by up to 37 %. An assembly of Ratatosk corrected reads from an Ashkenazi individual yields a contig N50 of 45 Mbp and less misassemblies than a PacBio HiFi reads assembly.

The structure, function and evolution of a complete human chromosome 8

The complete assembly of each human chromosome is essential for understanding human biology and evolution^1,2. Here we use complementary long-read sequencing technologies to complete the linear assembly of human chromosome 8. Our assembly resolves the sequence of five previously long-standing gaps, including a 2.08-Mb centromeric α-satellite array, a 644-kb copy number polymorphism in the β-defensin gene cluster that is important for disease risk, and an 863-kb variable number tandem repeat at chromosome 8q21.2 that can function as a neocentromere. We show that the centromeric α-satellite array is generally methylated except for a 73-kb hypomethylated region of diverse higher-order α-satellites enriched with CENP-A nucleosomes, consistent with the location of the kinetochore. In addition, we confirm the overall organization and methylation pattern of the centromere in a diploid human genome. Using a dual long-read sequencing approach, we complete high-quality draft assemblies of the orthologous centromere from chromosome 8 in chimpanzee, orangutan and macaque to reconstruct its evolutionary history. Comparative and phylogenetic analyses show that the higher-order α-satellite structure evolved in the great ape ancestor with a layered symmetry, in which more ancient higher-order repeats locate peripherally to monomeric α-satellites. We estimate that the mutation rate of centromeric satellite DNA is accelerated by more than 2.2-fold compared to the unique portions of the genome, and this acceleration extends into the flanking sequence.

Comparison of long-read methods for sequencing and assembly of a plant genome

BACKGROUND

Sequencing technologies have advanced to the point where it is possible to generate high-accuracy, haplotype-resolved, chromosome-scale assemblies. Several long-read sequencing technologies are available, and a growing number of algorithms have been developed to assemble the reads generated by those technologies. When starting a new genome project, it is therefore challenging to select the most cost-effective sequencing technology, as well as the most appropriate software for assembly and polishing. It is thus important to benchmark different approaches applied to the same sample.

RESULTS

Here, we report a comparison of 3 long-read sequencing technologies applied to the de novo assembly of a plant genome, Macadamia jansenii. We have generated sequencing data using Pacific Biosciences (Sequel I), Oxford Nanopore Technologies (PromethION), and BGI (single-tube Long Fragment Read) technologies for the same sample. Several assemblers were benchmarked in the assembly of Pacific Biosciences and Nanopore reads. Results obtained from combining long-read technologies or short-read and long-read technologies are also presented. The assemblies were compared for contiguity, base accuracy, and completeness, as well as sequencing costs and DNA material requirements.

CONCLUSIONS

The 3 long-read technologies produced highly contiguous and complete genome assemblies of M. jansenii. At the time of sequencing, the cost associated with each method was significantly different, but continuous improvements in technologies have resulted in greater accuracy, increased throughput, and reduced costs. We propose updating this comparison regularly with reports on significant iterations of the sequencing technologies.