Automation enables high-throughput and reproducible single-cell transcriptomics library preparation.

Next-generation sequencing (NGS) has revolutionized genomics, decreasing sequencing costs and allowing researchers to draw correlations between diseases and DNA or RNA changes. Technical advances have enabled the analysis of RNA expression changes between single cells within a heterogeneous population, known as single-cell RNA-seq (scRNA-seq). Despite resolving transcriptomes of cellular subpopulations, scRNA-seq has not replaced RNA-seq, due to higher costs and longer hands-on time. Here, we developed an automated workflow to increase throughput (up to 48 reactions) and to reduce by 75% the hands-on time of scRNA-seq library preparation, using the 10X Genomics Single Cell 3' kit. After gel bead-in-emulsion (GEM) generation on the 10X Genomics Chromium Controller, cDNA amplification was performed, and the product was normalized and subjected to either the manual, standard library preparation method or a fully automated, walk-away method using a Biomek i7 Hybrid liquid handler. Control metrics showed that both quantity and quality of the single-cell gene expression libraries generated were equivalent in size and yield. Key scRNA-seq downstream quality metrics, such as unique molecular identifiers count, mitochondrial RNA content, and cell and gene counts, further showed high correlations between automated and manual workflows. Using the UMAP dimensionality reduction technique to visualize all cells, we were able to further correlate the results observed between the manual and automated methods (R=0.971). The method developed here allows for the fast, error-free, and reproducible multiplex generation of high-quality single-cell gene expression libraries.

Automation of high-throughput mRNA-seq library preparation: a robust, hands-free and time efficient methodology.

Over the last decade, whole transcriptome profiling, also known as RNA-sequencing (RNA-seq), has quickly gained traction as a reliable method for unbiased assessment of gene expression. Integration of RNA-seq expression data into other omics datasets (e.g., proteomics, metabolomics, or epigenetics) solidifies our understanding of cell-specific regulatory patterns, yielding pathways to investigate the key rules of gene regulation. A limitation to efficient, at-scale utilization of RNA-seq is the time-demanding library preparation workflows, which is a 2-day or longer endeavor per cohort/sample size. To tackle this bottleneck, we designed an automated workflow that increases throughput capacity, while minimizing human error to enhance reproducibility. To this end, we converted the manual protocol of the NEBNext Directional Ultra II RNA Library Prep Kit for Illumina on the Beckman Coulter liquid handler, Biomek i7 Hybrid workstation. A total of 84 RNA samples were isolated from two human cell lines and subjected to comparative manual and automated library preparation methods. Qualitative and quantitative results indicated a high degree of similarity between libraries generated manually or through automation. Yet, there was a significant reduction in both hands-on and assay time from a 2-day manual to a 9-hour automated workflow. Using linear regression analysis, we found the Pearson correlation coefficient between libraries generated manually or by automation to be almost identical to a sample being sequenced twice (R²= 0.985 vs 0.983). This demonstrates that high-throughput automated workflows can be of great benefit to genomic laboratories by enhancing efficiency of library preparation, reducing hands-on time and increasing throughput potential.

Optimization of Automated Sample Preparation for Vitamin D Determination on a Biomek i7 Workstation.

Vitamin D belongs to the fat-soluble vitamins and is an integral part of bone metabolism. In the human body, a decreased vitamin D level can be an additional risk factor for diseases like cancer, diabetes, and mental diseases. As a result, an enormous increase in the demand for vitamin D testing has been observed in recent years, increasing the demand for powerful methods for vitamin D determination at the same time.Automation is the key factor in increasing sample throughput. This study compares three fully automated sample preparation methods for the determination of 25(OH)D2 and 25(OH)D3 in plasma and serum samples. Starting from a semiautomated reference method, the method is tested manually and subsequently fully automated on the Biomek i7 Workstation by integrating a centrifuge and a positive pressure extractor into the workstation. Alternatively, the centrifugation for the separation of protein aggregates and supernatant is replaced by a filter plate. Finally, the sample throughput is further increased by using phospholipid removal cartridges. The results show that phospholipid removal significantly increases the recovery rates in liquid chromatography-mass spectrometry. With the phospholipid removal cartridges, recovery rates of 97.36% for 25(OH)D2 and 102.5% for 25(OH)D3 were achieved, whereas with the automated classic automated preparation method, the recovery rates were 83.31% for 25(OH)D2 and 86.54% for 25(OH)D3. In addition to the technical evaluation, the different methods were also examined with regard to their economic efficiency. Finally, the qualitative and quantitative performance of the developed methods is benchmarked with a selected semiautomatic reference method.

A Plasma Sample Preparation for Mass Spectrometry using an Automated Workstation.

Sample preparation for mass spectrometry analysis in proteomics requires enzymatic cleavage of proteins into a peptide mixture. This process involves numerous incubation and liquid transfer steps in order to achieve denaturation, reduction, alkylation, and cleavage. Adapting this workflow onto an automated workstation can increase efficiency and reduce coefficients of variance, thereby providing more reliable data for statistical comparisons between sample types. We previously described an automated proteomic sample preparation workflow1. Here, we report the development of a more efficient and better controlled workflow with the following advantages: 1) The number of liquid transfer steps is reduced from nine to six by combining reagents; 2) Pipetting time is reduced by selective tip pipetting using a 96-position pipetting head with multiple channels; 3) Potential throughput is increased by the availability of up to 45 deck positions; 4) Complete enclosure of the system provides improved temperature and environmental control and reduces the potential for contamination of samples or reagents; and 5) The addition of stable isotope labeled peptides, as well as ß-galactosidase protein, to each sample makes monitoring and quality control possible throughout the entire process. These hardware and process improvements provide good reproducibility and improve intra-assay and inter-assay precision (CV of less than 20%) for LC-MS based protein and peptide quantification. The entire workflow for digesting 96 samples in a 96-well plate can be completed in approximately 5 hours.

The influence of trematode parasite burden on gene expression in a mammalian host.

BACKGROUND: Parasites can profoundly impact their hosts and are responsible for a plethora of debilitating diseases. To identify global changes in host gene expression related to parasite infection, we sequenced, assembled, and annotated the liver transcriptomes of Balb/cj mice infected with the trematode parasite Schistosoma mansoni and compared the results to uninfected mice. We used two different methodologies (i.e. de novo and reference guided) to evaluate the influence of parasite sequences on host transcriptome assembly. RESULTS: Our results demonstrate that the choice of assembly methodology significantly impacted the proportion of parasitic reads detected from the host library, yet the presence of non-target (xenobiotic) sequences did not create significant structural errors in the assembly. After removing parasite sequences from the mouse transcriptomes, we analyzed host gene expression under different parasite infection levels and observed significant differences in the associated immunologic and metabolic responses based on infection level. In particular, genes associated with T-helper type 1 (Th-1) and T-helper type 2 (Th-2) were up-regulated in infected mice whereas genes related to amino acid and carbohydrate metabolism were down-regulated in infected mice. These changes in gene expression scale with infection status and likely impact the evolutionary fitness of hosts. CONCLUSIONS: Overall, our data indicate that a) infected mice reduce the expression of key metabolic genes in direct proportion to their infection level; b) infected mice similarly increase the expression of key immune genes in response to infection; c) patterns of gene expression correspond to the pathological symptoms of schistosomiasis; and d) identifying and filtering out non-target sequences (xenobiotics) improves differential expression prediction. Our findings identify parasite targets for RNAi or other therapies and provide a better understanding of the pathology and host immune repertoire involved in response to S. mansoni infections.

The importance of including imperfect detection models in eDNA experimental design.

Environmental DNA (eDNA) is DNA that has been isolated from field samples, and it is increasingly used to infer the presence or absence of particular species in an ecosystem. However, the combination of sampling procedures and subsequent molecular amplification of eDNA can lead to spurious results. As such, it is imperative that eDNA studies include a statistical framework for interpreting eDNA presence/absence data. We reviewed published literature for studies that utilized eDNA where the species density was known and compared the probability of detecting the focal species to the sampling and analysis protocols. Although biomass of the target species and the volume per sample did not impact detectability, the number of field replicates and number of samples from each replicate were positively related to detection. Additionally, increased number of PCR replicates and increased primer specificity significantly increased detectability. Accordingly, we advocate for increased use of occupancy modelling as a method to incorporate effects of sampling effort and PCR sensitivity in eDNA study design. Based on simulation results and the hierarchical nature of occupancy models, we suggest that field replicates, as opposed to molecular replicates, result in better detection probabilities of target species.

Formalizing the definition of meta-analysis in Molecular Ecology.

Meta-analysis, the statistical synthesis of pertinent literature to develop evidence-based conclusions, is relatively new to the field of molecular ecology, with the first meta-analysis published in the journal Molecular Ecology in 2003 (Slate & Phua 2003). The goal of this article is to formalize the definition of meta-analysis for the authors, editors, reviewers and readers of Molecular Ecology by completing a review of the meta-analyses previously published in this journal. We also provide a brief overview of the many components required for meta-analysis with a more specific discussion of the issues related to the field of molecular ecology, including the use and statistical considerations of Wright's FST and its related analogues as effect sizes in meta-analysis. We performed a literature review to identify articles published as 'meta-analyses' in Molecular Ecology, which were then evaluated by at least two reviewers. We specifically targeted Molecular Ecology publications because as a flagship journal in this field, meta-analyses published in Molecular Ecology have the potential to set the standard for meta-analyses in other journals. We found that while many of these reviewed articles were strong meta-analyses, others failed to follow standard meta-analytical techniques. One of these unsatisfactory meta-analyses was in fact a secondary analysis. Other studies attempted meta-analyses but lacked the fundamental statistics that are considered necessary for an effective and powerful meta-analysis. By drawing attention to the inconsistency of studies labelled as meta-analyses, we emphasize the importance of understanding the components of traditional meta-analyses to fully embrace the strengths of quantitative data synthesis in the field of molecular ecology.

Horizontal gene transfer in schistosomes: A critical assessment.

Horizontal gene transfer (HGT), the movement of genetic material between distinct evolutionary lineages, has long been known as a principal force of diversification and adaptation of prokaryotes. More recently, genomic and transcriptomic datasets have suggested gene transfers among various eukaryotic taxa (e.g., Porifera, Cnidaria, Nematoda, Arthropoda, Rotifera, Craniata, and Plantae). Although the exact mechanism of eukaryotic HGT is often unknown, host-parasite interactions may provide ample opportunities for HGT. Schistosomes are trematode blood parasites with complex life cycles that have been repeatedly implicated in HGT. We employed molecular, bioinformatic and phylogenetic approaches to critically analyze 13 published reports of direct HGTs between schistosomes and their hosts to better understand host-parasite co-evolution. Our research suggests that reported cases of schistosome-associated HGT may be due to technical artifacts as opposed to biological reality as we were unable to substantiate them. HGT clearly occurs in eukaryotic organisms, but the burden of proof is high and we emphasize the importance of multiple lines of evidence to conclusively document HGT.

The genomic proliferation of transposable elements in colonizing populations: Schistosoma mansoni in the new world.

Transposable elements (TEs) are mobile genes with an inherent ability to move within and among genomes. Theory predicts that TEs proliferate extensively during physiological stress due to the breakdown of TE repression systems. We tested this hypothesis in Schistosoma mansoni, a widespread trematode parasite that causes the human disease schistosomiasis. According to phylogenetic analysis, S. mansoni invaded the new world during the last 500 years. We hypothesized that new world strains of S. mansoni would have more copies of TEs than old world strains due to the physiological stress associated with invasion of the new world. We quantified the copy number of six TEs (Saci-1, Saci-2 and Saci-3, Perere-1, Merlin-sm1, and SmTRC1) in the genome and the transcriptome of old world and new world strains of S. mansoni, using qPCR relative quantification. As predicted, the genomes of new world parasites contain significantly more copies of class I and class II TEs in both laboratory and field strains. However, such differences are not observed in the transcriptome suggesting that either TE silencing mechanisms have reactivated to control the expression of these elements or the presence of inactive truncated copies of TEs.

The genome sequence of a widespread apex predator, the golden eagle (Aquila chrysaetos).

Biologists routinely use molecular markers to identify conservation units, to quantify genetic connectivity, to estimate population sizes, and to identify targets of selection. Many imperiled eagle populations require such efforts and would benefit from enhanced genomic resources. We sequenced, assembled, and annotated the first eagle genome using DNA from a male golden eagle (Aquila chrysaetos) captured in western North America. We constructed genomic libraries that were sequenced using Illumina technology and assembled the high-quality data to a depth of ∼40x coverage. The genome assembly includes 2,552 scaffolds >10 Kb and 415 scaffolds >1.2 Mb. We annotated 16,571 genes that are involved in myriad biological processes, including such disparate traits as beak formation and color vision. We also identified repetitive regions spanning 92 Mb (∼6% of the assembly), including LINES, SINES, LTR-RTs and DNA transposons. The mitochondrial genome encompasses 17,332 bp and is ∼91% identical to the Mountain Hawk-Eagle (Nisaetus nipalensis). Finally, the data reveal that several anonymous microsatellites commonly used for population studies are embedded within protein-coding genes and thus may not have evolved in a neutral fashion. Because the genome sequence includes ∼800,000 novel polymorphisms, markers can now be chosen based on their proximity to functional genes involved in migration, carnivory, and other biological processes.

Hosts, parasites, and horizontal gene transfer.

Mendelian inheritance transfers genes vertically within lineages, whereas horizontal gene transfer (HGT) moves genetic material between or among lineages. Herein, we explore possible mechanisms of HGT between parasites and their hosts, as their intimate contact affords substantial opportunities for HGT. We review studies of host-parasite HGT, discussing their merits, their shortcomings, and the multiple lines of evidence needed to conclusively document HGT while avoiding false positives. We focus mainly on schistosomes and other parasites with complex life cycles as they provide numerous opportunities for HGT among the parasite, intermediate, and definitive host genomes. We also highlight future research directions that could prove illuminating with regard to the occurrence, prevalence, and overall importance of HGT in host-parasite coevolution.

Of contigs and quagmires: next-generation sequencing pitfalls associated with transcriptomic studies.

Molecular ecologists have good reasons to be excited about the newest DNA/RNA sequencing technologies. However, this exuberance should be tempered with a hefty dose of reality: new sequencing technologies come with significant new challenges. Herein, we offer a brief overview of some practical problems encountered during transcriptomics studies conducted in our laboratory, and of nontrivial issues that prospective practitioners should consider. These include template contamination (e.g. from xenobiotics) and the cutting-room floor problem, whereby most of the data are often unassembled, unannotated and unused. We also highlight computational requirements, including hardware, personnel time and associated skill sets. We are very optimistic about the future of molecular ecology, but we hope this cautionary overview will help neophytes better recognize some key challenges associated with new technologies.