Measuring sequencer size bias using REcount: a novel method for highly accurate Illumina sequencing-based quantification

Spontaneously regenerative corticospinal neurons in mice

The spinal cord receives inputs from the cortex via corticospinal neurons (CSNs). While predominantly a contralateral projection, a less-investigated minority of its axons terminate in the ipsilateral spinal cord. We analyzed the spatial and molecular properties of these ipsilateral axons and their post-synaptic targets in mice and found they project primarily to the ventral horn, including directly to motor neurons. Barcode-based reconstruction of the ipsilateral axons revealed a class of primarily bilaterally-projecting CSNs with a distinct cortical distribution. The molecular properties of these ipsilaterally-projecting CSNs (IP-CSNs) are strikingly similar to the previously described molecular signature of embryonic-like regenerating CSNs. Finally, we show that IP-CSNs are spontaneously regenerative after spinal cord injury. The discovery of a class of spontaneously regenerative CSNs may prove valuable to the study of spinal cord injury. Additionally, this work suggests that the retention of juvenile-like characteristics may be a widespread phenomenon in adult nervous systems.

Combining TSS-MPRA and sensitive TSS profile dissimilarity scoring to study the sequence determinants of transcription initiation

Cis-regulatory elements (CREs) can be classified by the shapes of their transcription start site (TSS) profiles, which are indicative of distinct regulatory mechanisms. Massively parallel reporter assays (MPRAs) are increasingly being used to study CRE regulatory mechanisms, yet the degree to which MPRAs replicate individual endogenous TSS profiles has not been determined. Here, we present a new low-input MPRA protocol (TSS-MPRA) that enables measuring TSS profiles of episomal reporters as well as after lentiviral reporter chromatinization. To sensitively compare MPRA and endogenous TSS profiles, we developed a novel dissimilarity scoring algorithm (WIP score) that outperforms the frequently used earth mover's distance on experimental data. Using TSS-MPRA and WIP scoring on 500 unique reporter inserts, we found that short (153 bp) MPRA promoter inserts replicate the endogenous TSS patterns of ∼60% of promoters. Lentiviral reporter chromatinization did not improve fidelity of TSS-MPRA initiation patterns, and increasing insert size frequently led to activation of extraneous TSS in the MPRA that are not active in vivo. We discuss the implications of our findings, which highlight important caveats when using MPRAs to study transcription mechanisms. Finally, we illustrate how TSS-MPRA and WIP scoring can provide novel insights into the impact of transcription factor motif mutations and genetic variants on TSS patterns and transcription levels.

Methods in Nutrition & Gut Microbiome Research: An American Society for Nutrition Satellite Session [13 October 2022]

The microbial cells colonizing the human body form an ecosystem that is integral to the regulation and maintenance of human health. Elucidation of specific associations between the human microbiome and health outcomes is facilitating the development of microbiome-targeted recommendations and treatments (e.g., fecal microbiota transplant; pre-, pro-, and post-biotics) to help prevent and treat disease. However, the potential of such recommendations and treatments to improve human health has yet to be fully realized. Technological advances have led to the development and proliferation of a wide range of tools and methods to collect, store, sequence, and analyze microbiome samples. However, differences in methodology at each step in these analytic processes can lead to variability in results due to the unique biases and limitations of each component. This technical variability hampers the detection and validation of associations with small to medium effect sizes. Therefore, the American Society for Nutrition (ASN) Nutritional Microbiology Group Engaging Members (GEM), sponsored by the Institute for the Advancement of Food and Nutrition Sciences (IAFNS), hosted a satellite session on methods in nutrition and gut microbiome research to review currently available methods for microbiome research, best practices, as well as tools and standards to aid in comparability of methods and results. This manuscript summarizes the topics and research discussed at the session. Consideration of the guidelines and principles reviewed in this session will increase the accuracy, precision, and comparability of microbiome research and ultimately the understanding of the associations between the human microbiome and health.

Clonally resolved single-cell multi-omics identifies routes of cellular differentiation in acute myeloid leukemia

Inter-patient variability and the similarity of healthy and leukemic stem cells (LSCs) have impeded the characterization of LSCs in acute myeloid leukemia (AML) and their differentiation landscape. Here, we introduce CloneTracer, a novel method that adds clonal resolution to single-cell RNA-seq datasets. Applied to samples from 19 AML patients, CloneTracer revealed routes of leukemic differentiation. Although residual healthy and preleukemic cells dominated the dormant stem cell compartment, active LSCs resembled their healthy counterpart and retained erythroid capacity. By contrast, downstream myeloid progenitors constituted a highly aberrant, disease-defining compartment: their gene expression and differentiation state affected both the chemotherapy response and leukemia's ability to differentiate into transcriptomically normal monocytes. Finally, we demonstrated the potential of CloneTracer to identify surface markers misregulated specifically in leukemic cells. Taken together, CloneTracer reveals a differentiation landscape that mimics its healthy counterpart and may determine biology and therapy response in AML.

Sequencing Illumina libraries at high accuracy on the ONT MinION using R2C2

High-throughput short-read sequencing has taken on a central role in research and diagnostics. Hundreds of different assays take advantage of Illumina short-read sequencers, the predominant short-read sequencing technology available today. Although other short-read sequencing technologies exist, the ubiquity of Illumina sequencers in sequencing core facilities and the high capital costs of these technologies have limited their adoption. Among a new generation of sequencing technologies, Oxford Nanopore Technologies (ONT) holds a unique position because the ONT MinION, an error-prone long-read sequencer, is associated with little to no capital cost. Here we show that we can make short-read Illumina libraries compatible with the ONT MinION by using the rolling circle to concatemeric consensus (R2C2) method to circularize and amplify the short library molecules. This results in longer DNA molecules containing tandem repeats of the original short library molecules. This longer DNA is ideally suited for the ONT MinION, and after sequencing, the tandem repeats in the resulting raw reads can be converted into high-accuracy consensus reads with similar error rates to that of the Illumina MiSeq. We highlight this capability by producing and benchmarking RNA-seq, ChIP-seq, and regular and target-enriched Tn5 libraries. We also explore the use of this approach for rapid evaluation of sequencing library metrics by implementing a real-time analysis workflow.

Genome Sequence Variations of Infectious Bronchitis Virus Serotypes From Commercial Chickens in Mexico

New variants of infectious bronchitis viruses (IBVs; Coronaviridae) continuously emerge despite routine vaccinations. Here, we report genome sequence variations of IBVs identified by random non-targeted next generation sequencing (NGS) of vaccine and field samples collected on FTA cards from commercial flocks in Mexico in 2019–2021. Paired-ended sequencing libraries prepared from rRNA-depleted RNAs were sequenced using Illumina MiSeq. IBV RNA was detected in 60.07% (n = 167) of the analyzed samples, from which 33 complete genome sequences were de novo assembled. The genomes are organized as 5'UTR-[Rep1a-Rep1b-S-3a-3b-E-M-4b-4c-5a-5b-N-6b]-3'UTR, except in eight sequences lacking non-structural protein genes (accessory genes) 4b, 4c, and 6b. Seventeen sequences have auxiliary S2' cleavage site located 153 residues downstream the canonically conserved primary furin-specific S1/S2 cleavage site. The sequences distinctly cluster into lineages GI-1 (Mass-type; n = 8), GI-3 (Holte/Iowa-97; n = 2), GI-9 (Arkansas-like; n = 8), GI-13 (793B; n = 14), and GI-17 (California variant; CAV; n = 1), with regional distribution in Mexico; this is the first report of the presence of 793B- and CAV-like strains in the country. Various point mutations, substitutions, insertions and deletions are present in the S1 hypervariable regions (HVRs I-III) across all 5 lineages, including in residues 38, 43, 56, 63, 66, and 69 that are critical in viral attachment to respiratory tract tissues. Nine intra-/inter-lineage recombination events are present in the S proteins of three Mass-type sequences, two each of Holte/Iowa-97 and Ark-like sequence, and one each of 793B-like and CAV-like sequences. This study demonstrates the feasibility of FTA cards as an attractive, adoptable low-cost sampling option for untargeted discovery of avian viral agents in field-collected clinical samples. Collectively, our data points to co-circulation of multiple distinct IBVs in Mexican commercial flocks, underscoring the need for active surveillance and a review of IBV vaccines currently used in Mexico and the larger Latin America region.

Optimization of enzymatic fragmentation is crucial to maximize genome coverage: a comparison of library preparation methods for Illumina sequencing

Background

Novel commercial kits for whole genome library preparation for next-generation sequencing on Illumina platforms promise shorter workflows, lower inputs and cost savings. Time savings are achieved by employing enzymatic DNA fragmentation and by combining end-repair and tailing reactions. Fewer cleanup steps also allow greater DNA input flexibility (1 ng-1 μg), PCR-free options from 100 ng DNA, and lower price as compared to the well-established sonication and tagmentation-based DNA library preparation kits.

Results

We compared the performance of four enzymatic fragmentation-based DNA library preparation kits (from New England Biolabs, Roche, Swift Biosciences and Quantabio) to a tagmentation-based kit (Illumina) using low input DNA amounts (10 ng) and PCR-free reactions with 100 ng DNA. With four technical replicates of each input amount and kit, we compared the kits’ fragmentation sequence-bias as well as performance parameters such as sequence coverage and the clinically relevant detection of single nucleotide and indel variants. While all kits produced high quality sequence data and demonstrated similar performance, several enzymatic fragmentation methods produced library insert sizes which deviated from those intended. Libraries with longer insert lengths performed better in terms of coverage, SNV and indel detection. Lower performance of shorter-insert libraries could be explained by loss of sequence coverage to overlapping paired-end reads, exacerbated by the preferential sequencing of shorter fragments on Illumina sequencers. We also observed that libraries prepared with minimal or no PCR performed best with regard to indel detection.

Conclusions

The enzymatic fragmentation-based DNA library preparation kits from NEB, Roche, Swift and Quantabio are good alternatives to the tagmentation based Nextera DNA flex kit from Illumina, offering reproducible results using flexible DNA inputs, quick workflows and lower prices. Libraries with insert DNA fragments longer than the cumulative sum of both read lengths avoid read overlap, thus produce more informative data that leads to strongly improved genome coverage and consequently also increased sensitivity and precision of SNP and indel detection. In order to best utilize such enzymatic fragmentation reagents, researchers should be prepared to invest time to optimize fragmentation conditions for their particular samples.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08316-y.

Measurement Biases Distort Cell-Free DNA Fragmentation Profiles and Define the Sensitivity of Metagenomic Cell-Free DNA Sequencing Assays

BACKGROUND

Metagenomic sequencing of microbial cell-free DNA (cfDNA) in blood and urine is increasingly used as a tool for unbiased infection screening. The sensitivity of metagenomic cfDNA sequencing assays is determined by the efficiency by which the assay recovers microbial cfDNA vs host-specific cfDNA. We hypothesized that the choice of methods used for DNA isolation, DNA sequencing library preparation, and sequencing would affect the sensitivity of metagenomic cfDNA sequencing.

METHODS

We characterized the fragment length biases inherent to select DNA isolation and library preparation procedures and developed a model to correct for these biases. We analyzed 305 cfDNA sequencing data sets, including publicly available data sets and 124 newly generated data sets, to evaluate the dependence of the sensitivity of metagenomic cfDNA sequencing on pre-analytical variables.

RESULTS

Length bias correction of fragment length distributions measured from different experimental procedures revealed the ultrashort (<100 bp) nature of microbial-, mitochondrial-, and host-specific urinary cfDNA. The sensitivity of metagenomic sequencing assays to detect the clinically reported microorganism differed by more than 5-fold depending on the combination of DNA isolation and library preparation used.

CONCLUSIONS

Substantial gains in the sensitivity of microbial and other short fragment recovery can be achieved by easy-to-implement changes in the sample preparation protocol, which highlights the need for standardization in the liquid biopsy field.

Customization of a dada2-based pipeline for fungal internal transcribed spacer 1 (ITS 1) amplicon datasets

Identification and analysis of fungal communities commonly rely on internal transcribed spacer–based (ITS-based) amplicon sequencing. There is no gold standard used to infer and classify fungal constituents since methodologies have been adapted from analyses of bacterial communities. To achieve high-resolution inference of fungal constituents, we customized a DADA2-based pipeline using a mix of 11 medically relevant fungi. While DADA2 allowed the discrimination of ITS1 sequences differing by single nucleotides, quality filtering, sequencing bias, and database selection were identified as key variables determining the accuracy of sample inference. Due to species-specific differences in sequencing quality, default filtering settings removed most reads that originated from Aspergillus species, Saccharomyces cerevisiae, and Candida glabrata. By fine-tuning the quality filtering process, we achieved an improved representation of the fungal communities. By adapting a wobble nucleotide in the ITS1 forward primer region, we further increased the yield of S. cerevisiae and C. glabrata sequences. Finally, we showed that a BLAST-based algorithm based on the UNITE+INSD or the NCBI NT database achieved a higher reliability in species-level taxonomic annotation compared with the naive Bayesian classifier implemented in DADA2. These steps optimized a robust fungal ITS1 sequencing pipeline that, in most instances, enabled species-level assignment of community members.

Performance Assessment of the Devyser High-Throughput Sequencing-Based Assay for Chimerism Monitoring in Patients after Allogeneic Hematopoietic Stem Cell Transplantation

Chimerism analysis is widely used to aid in the clinical management of patients after allogeneic hematopoietic stem cell transplantation. Many laboratories currently use assays based on polymerase chain reaction followed by capillary electrophoresis, with a limit of quantification of 1% to 5%. Assays with a lower limit of quantification could allow for earlier relapse detection, resulting in improved patient care. This study investigated the analytical, clinical, technical, and practical performance of the Devyser NGS chimerism assay, a commercial high-throughput sequencing-based assay for chimerism analysis. Performance of this assay was compared with that of the Promega PowerPlex 16 HS assay, a commercial capillary electrophoresis-based assay. A limit of quantification of 0.1% was achievable with the Devyser assay. The repeatability, reproducibility, trueness, and linearity of the Devyser assay were acceptable. The Devyser assay showed potential for earlier relapse detection compared with the Promega assay. Conclusive analysis was not possible for 3% of donor-recipient pairs with the Devyser assay due to an insufficient number of informative markers; this factor was not an issue for the Promega assay. Further improvements in assay design or data analysis may allow the assay's applicability to be extended to all donor-recipient pairs studied. Technical performance criteria for chimerism analysis by high-throughput sequencing were suggested and evaluated. Both assays were found to be practical for use in a clinical diagnostics laboratory.

Dissecting and tuning primer editing by proofreading polymerases

Proofreading polymerases have 3′ to 5′ exonuclease activity that allows the excision and correction of mis-incorporated bases during DNA replication. In a previous study, we demonstrated that in addition to correcting substitution errors and lowering the error rate of DNA amplification, proofreading polymerases can also edit PCR primers to match template sequences. Primer editing is a feature that can be advantageous in certain experimental contexts, such as amplicon-based microbiome profiling. Here we develop a set of synthetic DNA standards to report on primer editing activity and use these standards to dissect this phenomenon. The primer editing standards allow next-generation sequencing-based enzymological measurements, reveal the extent of editing, and allow the comparison of different polymerases and cycling conditions. We demonstrate that proofreading polymerases edit PCR primers in a concentration-dependent manner, and we examine whether primer editing exhibits any sequence specificity. In addition, we use these standards to show that primer editing is tunable through the incorporation of phosphorothioate linkages. Finally, we demonstrate the ability of primer editing to robustly rescue the drop-out of taxa with 16S rRNA gene-targeting primer mismatches using mock communities and human skin microbiome samples.

Quantification, Dynamic Visualization, and Validation of Bias in ATAC-Seq Data with ataqv

The assay for transposase-accessible chromatin using sequencing (ATAC-seq) has become the preferred method for mapping chromatin accessibility due to its time and input material efficiency. However, it can be difficult to evaluate data quality and identify sources of technical bias across samples. Here, we present ataqv, a computational toolkit for efficiently measuring, visualizing, and comparing quality control (QC) results across samples and experiments. We use ataqv to analyze 2,009 public ATAC-seq datasets; their QC metrics display a 10-fold range. Tn5 dosage experiments and statistical modeling show that technical variation in the ratio of Tn5 transposase to nuclei and sequencing flowcell density induces systematic bias in ATAC-seq data by changing the enrichment of reads across functional genomic annotations including promoters, enhancers, and transcription-factor-bound regions, with the notable exception of CTCF. ataqv can be integrated into existing computational pipelines and is freely available at https://github.com/ParkerLab/ataqv/.

An optimized protocol for rapid, sensitive and robust on-bead ChIP-seq from primary cells

Integrative analysis of next-generation sequencing data can help understand disease mechanisms. Specifically, ChIP-seq can illuminate where transcription regulators bind to regulate transcription. A major obstacle to performing this assay on primary cells is the low numbers obtained from tissues. The extensively validated ChIP-seq protocol presented here uses small volumes and single-pot on-bead library preparation to generate diverse high-quality ChIP-seq data. This protocol allows for medium-to-high-throughput ChIP-seq of low-abundance cells and can also be applied to other mammalian cells. For complete details on the use and execution of this protocol, please refer to Brigidi et al. (2019), Carlin et al. (2018), Heinz et al. (2018), Nott et al. (2019), Sakai et al. (2019), and Seidman et al. (2020).

A rapid, cost-effective tailed amplicon method for sequencing SARS-CoV-2

Background

The global COVID-19 pandemic has led to an urgent need for scalable methods for clinical diagnostics and viral tracking. Next generation sequencing technologies have enabled large-scale genomic surveillance of SARS-CoV-2 as thousands of isolates are being sequenced around the world and deposited in public data repositories. A number of methods using both short- and long-read technologies are currently being applied for SARS-CoV-2 sequencing, including amplicon approaches, metagenomic methods, and sequence capture or enrichment methods. Given the small genome size, the ability to sequence SARS-CoV-2 at scale is limited by the cost and labor associated with making sequencing libraries.

Results

Here we describe a low-cost, streamlined, all amplicon-based method for sequencing SARS-CoV-2, which bypasses costly and time-consuming library preparation steps. We benchmark this tailed amplicon method against both the ARTIC amplicon protocol and sequence capture approaches and show that an optimized tailed amplicon approach achieves comparable amplicon balance, coverage metrics, and variant calls to the ARTIC v3 approach.

Conclusions

The tailed amplicon method we describe represents a cost-effective and highly scalable method for SARS-CoV-2 sequencing.

Plant DNA barcoding necessitates marker-specific efforts to establish more comprehensive reference databases

The problem of low species-level identification rates in plants by DNA barcoding is exacerbated by the fact that reference databases are far from being comprehensive. We investigate the impact of increased sampling depth on identification success by analyzing the efficacy of established plant barcode marker sequences (rbcL, matK, trnL-trnF, psbA-trnH, ITS). Adding sequences of the same species to the reference database led to an increase in correct species assignment of +10.9% for rbcL and +19.0% for ITS. Simultaneously, erroneous identification dropped from ∼40% to ∼12.5%. Despite its evolutionary constraints, ITS showed the highest identification rate and identification gain by increased sampling effort, which makes it a very suitable marker in the planning phase of a barcode study. The limited sequence availability of trnL-trnF is problematic for an otherwise very promising plastid plant barcoding marker. Future developments in machine learning algorithms have the potential to give new impetus to plant barcoding, but are dependent on extensive reference databases. We expect that our results will be incorporated into future plans for the development of DNA barcoding reference databases and will lead to these being developed with greater depth and taxonomic coverage.

Minority report: the intestinal mycobiota in systemic infections

Compared to bacteria, fungi often exhibit a lower abundance and a higher temporal volatility in the intestinal microbiota. Analysis of fungi in the microbiota (mycobiota) faces technical limitations with tools that were originally developed for analyzing bacteria. Dysbiotic states of the intestinal mycobiota, often associated with disruption of the healthy bacterial microbiota, are characterized by overgrowth (domination) of specific fungal taxa and loss of diversity. Intestinal domination by Candida species has been shown to be a major source of Candida bloodstream infections. Fungal dysbiosis is also linked to the development and treatment response in non-fungal infections, for example Clostridioides difficile colitis and HIV. Further research is needed to define the contribution of intestinal mycobiota to human fungal and non-fungal infections.

Adapterama II: universal amplicon sequencing on Illumina platforms (TaggiMatrix)

Next-generation sequencing (NGS) of amplicons is used in a wide variety of contexts. In many cases, NGS amplicon sequencing remains overly expensive and inflexible, with library preparation strategies relying upon the fusion of locus-specific primers to full-length adapter sequences with a single identifying sequence or ligating adapters onto PCR products. In Adapterama I, we presented universal stubs and primers to produce thousands of unique index combinations and a modifiable system for incorporating them into Illumina libraries. Here, we describe multiple ways to use the Adapterama system and other approaches for amplicon sequencing on Illumina instruments. In the variant we use most frequently for large-scale projects, we fuse partial adapter sequences (TruSeq or Nextera) onto the 5' end of locus-specific PCR primers with variable-length tag sequences between the adapter and locus-specific sequences. These fusion primers can be used combinatorially to amplify samples within a 96-well plate (8 forward primers + 12 reverse primers yield 8 × 12 = 96 combinations), and the resulting amplicons can be pooled. The initial PCR products then serve as template for a second round of PCR with dual-indexed iTru or iNext primers (also used combinatorially) to make full-length libraries. The resulting quadruple-indexed amplicons have diversity at most base positions and can be pooled with any standard Illumina library for sequencing. The number of sequencing reads from the amplicon pools can be adjusted, facilitating deep sequencing when required or reducing sequencing costs per sample to an economically trivial amount when deep coverage is not needed. We demonstrate the utility and versatility of our approaches with results from six projects using different implementations of our protocols. Thus, we show that these methods facilitate amplicon library construction for Illumina instruments at reduced cost with increased flexibility. A simple web page to design fusion primers compatible with iTru primers is available at: http://baddna.uga.edu/tools-taggi.html. A fast and easy to use program to demultiplex amplicon pools with internal indexes is available at: https://github.com/lefeverde/Mr_Demuxy.

Deconvolution of nucleic-acid length distributions: a gel electrophoresis analysis tool and applications

Next-generation DNA-sequencing (NGS) technologies, which are designed to streamline the acquisition of massive amounts of sequencing data, are nonetheless dependent on various preparative steps to generate DNA fragments of required concentration, purity and average size (molecular weight). Current automated electrophoresis systems for DNA- and RNA-sample quality control, such as Agilent's Bioanalyzer® and TapeStation® products, are costly to acquire and use; they also provide limited information for samples having broad size distributions. Here, we describe a software tool that helps determine the size distribution of DNA fragments in an NGS library, or other DNA sample, based on gel-electrophoretic line profiles. The software, developed as an ImageJ plug-in, allows for straightforward processing of gel images, including lane selection and fitting of univariate functions to intensity distributions. The user selects the option of fitting either discrete profiles in cases where discrete gel bands are visible or continuous profiles, having multiple bands buried under a single broad peak. The method requires only modest imaging capabilities and is a cost-effective, rigorous alternative characterization method to augment existing techniques for library quality control.

Measuring sequencer size bias using REcount: a novel method for highly accurate Illumina sequencing-based quantification

Quantification of DNA sequence tags from engineered constructs such as plasmids, transposons, or other transgenes underlies many functional genomics measurements. Typically, such measurements rely on PCR followed by next-generation sequencing. However, PCR amplification can introduce significant quantitative error. We describe REcount, a novel PCR-free direct counting method. Comparing measurements of defined plasmid pools to droplet digital PCR data demonstrates that REcount is highly accurate and reproducible. We use REcount to provide new insights into clustering biases due to molecule length across different Illumina sequencers and illustrate the impacts on interpretation of next-generation sequencing data and the economics of data generation.