Microarray-based genomic selection for high-throughput resequencing

Mini-review of DNA Methylation Detection Techniques and Their Potential Applications in Disease Diagnosis, Prognosis, and Treatment

DNA methylation is the dominant epigenetic mechanism for regulating gene expression in mammals, playing crucial roles in development, differentiation, and tissue homeostasis. Aberrations in DNA methylation are closely associated with the potential onset of various diseases. Consequently, numerous DNA methylation detection techniques have been successively developed. These methods not only facilitate the exploration of disease mechanisms but also hold significant promise for the development of diagnostic and prognostic strategies. In this Perspective, we present a comprehensive overview of commonly employed DNA methylation detection techniques as well as biosensing based on their underlying analytical techniques. For its medical applications, we begin by examining the pathogenesis of different diseases and then proceed to discuss how relevant technologies are applied in the context of these specific medical conditions. Additionally, we briefly discuss the current limitations of these techniques and highlight future challenges in advancing methylation detection and analysis methodologies.

Target Enrichment Approaches for Next-Generation Sequencing Applications in Oncology

Screening for genomic sequence variants in genes of predictive and prognostic significance is an integral part of precision medicine. Next-generation sequencing (NGS) technologies are progressively becoming platforms of choice to facilitate this, owing to their massively parallel sequencing capability, which can be used to simultaneously screen multiple markers in multiple samples for a variety of variants (single nucleotide and multi nucleotide variants, insertions and deletions, gene copy number variations, and fusions). A crucial step in the workflow of targeted NGS is the enrichment of the genomic regions of interest to be sequenced, against the whole genomic background. This ensures that the NGS effort is focused to predominantly screen target regions of interest with minimal off-target sequencing, making it more accurate and economical. Polymerase chain reaction-based (PCR, or amplicon-based) and hybridization capture-based methodologies are the two prominent approaches employed for target enrichment. This review summarizes the basic principles of target enrichment utilized by these methods, their multiple variations that have evolved over time, automation approaches, overall comparison of their advantages and drawbacks, and commercially available choices for these methodologies.

Progress in forensic bone DNA analysis: Lessons learned from ancient DNA

Research on ancient and forensic DNA is related in many ways, and the two fields must deal with similar obstacles. Therefore, communication between these two communities has the potential to improve results in both research fields. Here, we present the insights gained in the ancient DNA community with regard to analyzing DNA from aged skeletal material and the potential use of the developed protocols in forensic work. We discuss the various steps, from choosing samples for DNA extraction to deciding between classical PCR amplification and massively parallel sequencing approaches. Based on the progress made in ancient DNA analyses combined with the requirements of forensic work, we suggest that there is substantial potential for incorporating ancient DNA approaches into forensic protocols, a process that has already begun to a considerable extent. However, taking full advantage of the experiences gained from ancient DNA work will require comparative studies by the forensic DNA community to tailor the methods developed for ancient samples to the specific needs of forensic studies and case work. If successful, in our view, the benefits for both communities would be considerable.

Comparative Analysis of SNP Discovery and Genotyping in Fagus sylvatica L. and Quercus robur L. Using RADseq, GBS, and ddRAD Methods

Next-generation sequencing of reduced representation genomic libraries (RRL) is capable of providing large numbers of genetic markers for population genetic studies at relatively low costs. However, one major concern of these types of markers is the precision of genotyping, which is related to the common problem of missing data, which appears to be particularly important in association and genomic selection studies. We evaluated three RRL approaches (GBS, RADseq, ddRAD) and different SNP identification methods (de novo or based on a reference genome) to find the best solutions for future population genomics studies in two economically and ecologically important broadleaved tree species, namely F. sylvatica and Q. robur. We found that the use of ddRAD method coupled with SNP calling based on reference genomes provided the largest numbers of markers (28 k and 36 k for beech and oak, respectively), given standard filtering criteria. Using technical replicates of samples, we demonstrated that more than 80% of SNP loci should be considered as reliable markers in GBS and ddRAD, but not in RADseq data. According to the reference genomes’ annotations, more than 30% of the identified ddRAD loci appeared to be related to genes. Our findings provide a solid support for using ddRAD-based SNPs for future population genomics studies in beech and oak.

CRISPR-based enrichment strategies for targeted sequencing

The ability to easily produce or procure sequencing data has expanded to be within the reach of most clinics and research laboratories, but the complexity of sequence analysis remains a hurdle for many scientists, and a decline in sequencing cost means that the generation of gratuitous information in a given experiment is a challenge that is more and more often being encountered. To address this issue, methods have been present, some dating to the advent of nucleic acid sequencing, for capturing, targeting, or otherwise enriching specific nucleic acids in order to obtain greater depth of reads from a small portion of sequences within a complex sample. However, many of these methods have been complicated and laborious, relying on the design of hundreds to thousands of oligonucleotide probes, fabrication of microarray chips, and long hybridization times. Here, we review these methods, their benefits and uses, and catalog and discuss the implications of a recent development that has enabled a more efficient and expanded set of tools for enriching nucleic acids - the application of CRISPR technology. This introduction and analysis of the capabilities of new CRISPR-based enrichment strategies shows that it has the potential to expand the scope of enrichment to new possibilities, including the coupling of DNA and RNA targeting with long-read, portable sequencing platforms. Moreover, there are several areas where CRISPR-enrichment is a logical next step to more powerful and simplified sequencing for applications such as diagnostics and environmental monitoring.

EYS is a major gene involved in retinitis pigmentosa in Japan: genetic landscapes revealed by stepwise genetic screening

Next-generation sequencing (NGS) has greatly advanced the studies of causative genes and variants of inherited diseases. While it is sometimes challenging to determine the pathogenicity of identified variants in NGS, the American College of Medical Genetics and Genomics established the guidelines to help the interpretation. However, as to the genetic screenings for patients with retinitis pigmentosa (RP) in Japan, none of the previous studies utilized the guidelines. Considering that EYS is the major causative gene of RP in Japan, we conducted stepwise genetic screening of 220 Japanese patients with RP utilizing the guidelines. Step 1-4 comprised the following, in order: Sanger sequencing for two major EYS founder mutations; targeted sequencing of all coding regions of EYS; whole genome sequencing; Sanger sequencing for Alu element insertion in RP1, a recently determined founder mutation for RP. Among the detected variants, 2, 19, 173, and 1 variant(s) were considered pathogenic and 8, 41, 44, and 5 patients were genetically solved in step 1, 2, 3, and 4, respectively. Totally, 44.5% (98/220) of the patients were genetically solved, and 50 (51.0%) were EYS-associated and 5 (5.1%) were Alu element-associated. Among the unsolved 122 patients, 22 had at least one possible pathogenic variant.

AnthOligo: automating the design of oligonucleotides for capture/enrichment technologies

SUMMARY

A number of methods have been devised to address the need for targeted genomic resequencing. One of these methods, region-specific extraction (RSE) is characterized by the capture of long DNA fragments (15-20 kb) by magnetic beads, after enzymatic extension of oligonucleotides hybridized to selected genomic regions. Facilitating the selection of the most appropriate capture oligos for targeting a region of interest, satisfying the properties of temperature (Tm) and entropy (ΔG), while minimizing the formation of primer-dimers in a pooled experiment, is therefore necessary. Manual design and selection of oligos becomes very challenging, complicated by factors such as length of the target region and number of targeted regions. Here we describe, AnthOligo, a web-based application developed to optimally automate the process of generation of oligo sequences used to target and capture the continuum of large and complex genomic regions. Apart from generating oligos for RSE, this program may have wider applications in the design of customizable internal oligos to be used as baits for gene panel analysis or even probes for large-scale comparative genomic hybridization array processes. AnthOligo was tested by capturing the Major Histocompatibility Complex (MHC) of a random sample.The application provides users with a simple interface to upload an input file in BED format and customize parameters for each task. The task of probe design in AnthOligo commences when a user uploads an input file and concludes with the generation of a result-set containing an optimal set of region-specific oligos. AnthOligo is currently available as a public web application with URL: http://antholigo.chop.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The past, present and future of ancient bacterial DNA

Groundbreaking studies conducted in the mid-1980s demonstrated the possibility of sequencing ancient DNA (aDNA), which has allowed us to answer fundamental questions about the human past. Microbiologists were thus given a powerful tool to glimpse directly into inscrutable bacterial history, hitherto inaccessible due to a poor fossil record. Initially plagued by concerns regarding contamination, the field has grown alongside technical progress, with the advent of high-throughput sequencing being a breakthrough in sequence output and authentication. Albeit burdened with challenges unique to the analysis of bacteria, a growing number of viable sources for aDNA has opened multiple avenues of microbial research. Ancient pathogens have been extracted from bones, dental pulp, mummies and historical medical specimens and have answered focal historical questions such as identifying the aetiological agent of the black death as Yersinia pestis. Furthermore, ancient human microbiomes from fossilized faeces, mummies and dental plaque have shown shifts in human commensals through the Neolithic demographic transition and industrial revolution, whereas environmental isolates stemming from permafrost samples have revealed signs of ancient antimicrobial resistance. Culminating in an ever-growing repertoire of ancient genomes, the quickly expanding body of bacterial aDNA studies has also enabled comparisons of ancient genomes to their extant counterparts, illuminating the evolutionary history of bacteria. In this review we summarize the present avenues of research and contextualize them in the past of the field whilst also pointing towards questions still to be answered.

A Guide to Carrying Out a Phylogenomic Target Sequence Capture Project

High-throughput DNA sequencing techniques enable time- and cost-effective sequencing of large portions of the genome. Instead of sequencing and annotating whole genomes, many phylogenetic studies focus sequencing effort on large sets of pre-selected loci, which further reduces costs and bioinformatic challenges while increasing coverage. One common approach that enriches loci before sequencing is often referred to as target sequence capture. This technique has been shown to be applicable to phylogenetic studies of greatly varying evolutionary depth. Moreover, it has proven to produce powerful, large multi-locus DNA sequence datasets suitable for phylogenetic analyses. However, target capture requires careful considerations, which may greatly affect the success of experiments. Here we provide a simple flowchart for designing phylogenomic target capture experiments. We discuss necessary decisions from the identification of target loci to the final bioinformatic processing of sequence data. We outline challenges and solutions related to the taxonomic scope, sample quality, and available genomic resources of target capture projects. We hope this review will serve as a useful roadmap for designing and carrying out successful phylogenetic target capture studies.

Liquid biopsy for pediatric diffuse midline glioma: a review of circulating tumor DNA and cerebrospinal fluid tumor DNA

Diffuse midline glioma (DMG) is a highly malignant childhood tumor with an exceedingly poor prognosis and limited treatment options. The majority of these tumors harbor somatic mutations in genes encoding histone variants. These recurrent mutations correlate with treatment response and are forming the basis for molecularly guided clinical trials. The ability to detect these mutations, either in circulating tumor DNA (ctDNA) or cerebrospinal fluid tumor DNA (CSF-tDNA), may enable noninvasive molecular profiling and earlier prediction of treatment response. Here, the authors review ctDNA and CSF-tDNA detection methods, detail recent studies that have explored detection of ctDNA and CSF-tDNA in patients with DMG, and discuss the implications of liquid biopsies for patients with DMG.

Identifying MicroRNA Pathway Variants as Biomarkers of Patient Selection for Immune Therapy

In this chapter we discuss the discovery and validation of microRNA (miRNA) associated germline biomarkers, as well as their application on a cohort of patients treated with immune therapy to predict response and toxicity. MiRNAs are the first class of noncoding RNAs discovered, and these pathways have been shown to be important regulators of the systemic stress response, including that to cancer therapy. We detail the original discovery efforts identifying germline biomarkers that disrupt miRNA circuitry, and then the selection, application, and validation of these biomarkers and their potential to predict important outcomes to checkpoint therapy.

Phylogenetic relationships and systematics of the Amazonian poison frog genus Ameerega using ultraconserved genomic elements

The Amazonian poison frog genus Ameerega is one of the largest yet most understudied of the brightly colored genera in the anuran family Dendrobatidae, with 30 described species ranging throughout tropical South America. Phylogenetic analyses of Ameerega are highly discordant, lacking consistency due to variation in data types and methods, and often with limited coverage of species diversity in the genus. Here, we present a comprehensive phylogenomic reconstruction of Ameerega, utilizing state-of-the-art sequence capture techniques and phylogenetic methods. We sequenced thousands of ultraconserved elements from over 100 tissue samples, representing almost every described Ameerega species, as well as undescribed cryptic diversity. We generated topologies using maximum likelihood and coalescent methods and compared the use of maximum likelihood and Bayesian methods for estimating divergence times. Our phylogenetic inference diverged strongly from those of previous studies, and we recommend steps to bring Ameerega taxonomy in line with the new phylogeny. We place several species in a phylogeny for the first time, as well as provide evidence for six potential candidate species. We estimate that Ameerega experienced a rapid radiation approximately 7-11 million years ago and that the ancestor of all Ameerega was likely an aposematic, montane species. This study underscores the utility of phylogenomic data in improving our understanding of the phylogeny of understudied clades and making novel inferences about their evolution.

Nonfunctionalized PNAs as Beacons for Nucleic Acid Detection in a Nanopore System

In this work, single-channel current recordings were used to selectively detect individual ssDNA strands in the vestibule of the α-hemolysin (α-HL) protein nanopore. The sensing mechanism was based on the detection of the intrinsic topological change of target ssDNA molecules after the hybridization with complementary PNA fragments. The readily distinguishable current signatures of PNA-DNA duplexes reversible association with the α-HL's vestibule, in terms of blockade amplitudes and kinetic features, allows specific detection of nucleic acid hybridization.

Gene Saturation: An Approach to Assess Exploration Stage of Gene Interaction Networks

The gene interaction network is one of the most important biological networks and has been studied by many researchers. The gene interaction network provides information about whether the genes in the network can cause or heal diseases. As gene-gene interaction relations are constantly explored, gene interaction networks are evolving. To describe how much a gene has been studied, an approach based on a logistic model for each gene called gene saturation has been proposed, which in most cases, satisfies non-decreasing, correlation and robustness principles. The average saturation of a group of genes can be used to assess the network constructed by these genes. Saturation reflects the distance between known gene interaction networks and the real gene interaction network in a cell. Furthermore, the saturation values of 546 disease gene networks that belong to 15 categories of diseases have been calculated. The disease gene networks’ saturation for cancer is significantly higher than that of all other diseases, which means that the disease gene networks’ structure for cancer has been more deeply studied than other disease. Gene saturation provides guidance for selecting an experimental subject gene, which may have a large number of unknown interactions.

Comprehensive assay for the molecular profiling of cancer by target enrichment from formalin-fixed paraffin-embedded specimens

Tumor molecular profiling is becoming a standard of care for patients with cancer, but the optimal platform for cancer sequencing remains undetermined. We established a comprehensive assay, the Todai OncoPanel (TOP), which consists of DNA and RNA hybridization capture‐based next‐generation sequencing panels. A novel method for target enrichment, named the junction capture method, was developed for the RNA panel to accurately and cost‐effectively detect 365 fusion genes as well as aberrantly spliced transcripts. The TOP RNA panel can also measure the expression profiles of an additional 109 genes. The TOP DNA panel was developed to detect single nucleotide variants and insertions/deletions for 464 genes, to calculate tumor mutation burden and microsatellite instability status, and to infer chromosomal copy number. Clinically relevant somatic mutations were identified in 32.2% (59/183) of patients by prospective TOP testing, signifying the clinical utility of TOP for providing personalized medicine to cancer patients.

Hybrid Capture-Based Next Generation Sequencing and Its Application to Human Infectious Diseases

This review describes target-enrichment approaches followed by next generation sequencing and their recent application to the research and diagnostic field of modern and past infectious human diseases caused by viruses, bacteria, parasites and fungi.

Design and validation of a universal influenza virus enrichment probe set and its utility in deep sequence analysis of primary cloacal swab surveillance samples of wild birds

Influenza virus infections in humans and animals are major public health concerns. In the current study, a set of universal influenza enrichment probes was developed to increase the sensitivity of sequence-based virus detection and characterization for all influenza viruses. This universal influenza enrichment probe set contains 46,953 120nt RNA biotin-labeled probes designed based on all available influenza viral sequences and it can be used to enrich for influenza sequences without prior knowledge of type or subtype. Marked enrichment was demonstrated in influenza A/H1N1, influenza B, and H1-to-H16 hemagglutinin plasmids spiked into human DNA and in cultured influenza A/H2N1 virus. Furthermore, enrichment effects and mixed influenza A virus infections were revealed in wild bird cloacal swab samples. Therefore, this universal influenza virus enrichment probe system can capture and enrich influenza viral sequences selectively and effectively in different samples, especially ones with degraded RNA or containing low amount of influenza RNA.

Targeted resequencing of phosphorus metabolism‑related genes in 86 patients with hypophosphatemic rickets/osteomalacia

Hypophosphatemic rickets/osteomalacia is characterized by defective renal phosphate reabsorption and abnormal bone mineralization. Hypophosphatemic rickets/osteomalacia consists of inherited and acquired forms, many of which have unknown aetiology. In the present study, next‑generation sequencing‑based resequencing was used on samples from Chinese subjects with hypophosphatemic rickets/osteomalacia, aiming to detect the spectrum of pathogenic genes in these patients. A total of 86 hypophosphatemic rickets/osteomalacia patients (ranging from 3 to 70 years old) were recruited. Patients with tumour‑induced osteomalacia (TIO), renal tubular acidosis, renal osteodystrophy, and adefovir‑induced Fanconi syndrome were excluded. Targeted massively parallel resequencing of 196 candidate genes for hypophosphatemic rickets/osteomalacia was performed in the 86 affected unrelated individuals (cases) and in 100 unrelated healthy controls to identify new genes and mutations in known genes that cause hypophosphatemic rickets/osteomalacia. The results identified seven phosphate‑regulating gene with homologies to endopeptidases on the X chromosome (PHEX) mutations (of which two were novel) and one novel dentin matrix protein 1 (DMP1) mutation in eight patients. Following targeted exome sequencing data analysis, 14 candidate disease‑related gene loci were selected, two of which were of most concern regarding disease severity. Further validation of the present results is warranted, with additional sequencing projects and functional tests. To our knowledge, the present study is the largest cohort of cases with hypophosphatemic rickets/osteomalacia to undergo targeted resequencing. The diagnosis and understanding of the molecular aetiologies of these disorders will be improved by this fast and efficient approach.

Comprehensive viral enrichment enables sensitive respiratory virus genomic identification and analysis by next generation sequencing

Next generation sequencing (NGS) technologies have revolutionized the genomics field and are becoming more commonplace for identification of human infectious diseases. However, due to the low abundance of viral nucleic acids (NAs) in relation to host, viral identification using direct NGS technologies often lacks sufficient sensitivity. Here, we describe an approach based on two complementary enrichment strategies that significantly improves the sensitivity of NGS-based virus identification. To start, we developed two sets of DNA probes to enrich virus NAs associated with respiratory diseases. The first set of probes spans the genomes, allowing for identification of known viruses and full genome sequencing, while the second set targets regions conserved among viral families or genera, providing the ability to detect both known and potentially novel members of those virus groups. Efficiency of enrichment was assessed by NGS testing reference virus and clinical samples with known infection. We show significant improvement in viral identification using enriched NGS compared to unenriched NGS. Without enrichment, we observed an average of 0.3% targeted viral reads per sample. However, after enrichment, 50%–99% of the reads per sample were the targeted viral reads for both the reference isolates and clinical specimens using both probe sets. Importantly, dramatic improvements on genome coverage were also observed following virus-specific probe enrichment. The methods described here provide improved sensitivity for virus identification by NGS, allowing for a more comprehensive analysis of disease etiology.

Methylation-based enrichment facilitates low-cost, noninvasive genomic scale sequencing of populations from feces

Obtaining high-quality samples from wild animals is a major obstacle for genomic studies of many taxa, particularly at the population level, as collection methods for such samples are typically invasive. DNA from feces is easy to obtain noninvasively, but is dominated by bacterial and other non-host DNA. The high proportion of non-host DNA drastically reduces the efficiency of high-throughput sequencing for host animal genomics. To address this issue, we developed an inexpensive capture method for enriching host DNA from noninvasive fecal samples. Our method exploits natural differences in CpG-methylation density between vertebrate and bacterial genomes to preferentially bind and isolate host DNA from majority-bacterial samples. We demonstrate that the enrichment is robust, efficient, and compatible with downstream library preparation methods useful for population studies (e.g., RADseq). Compared to other enrichment strategies, our method is quick and inexpensive, adding only a negligible cost to sample preparation. In combination with downstream methods such as RADseq, our approach allows for cost-effective and customizable genomic-scale genotyping that was previously feasible in practice only with invasive samples. Because feces are widely available and convenient to collect, our method empowers researchers to explore genomic-scale population-level questions in organisms for which invasive sampling is challenging or undesirable.

High-sensitivity HLA typing by Saturated Tiling Capture Sequencing (STC-Seq)

Background

Highly polymorphic human leukocyte antigen (HLA) genes are responsible for fine-tuning the adaptive immune system. High-resolution HLA typing is important for the treatment of autoimmune and infectious diseases. Additionally, it is routinely performed for identifying matched donors in transplantation medicine. Although many HLA typing approaches have been developed, the complexity, low-efficiency and high-cost of current HLA-typing assays limit their application in population-based high-throughput HLA typing for donors, which is required for creating large-scale databases for transplantation and precision medicine.

Results

Here, we present a cost-efficient Saturated Tiling Capture Sequencing (STC-Seq) approach to capturing 14 HLA class I and II genes. The highly efficient capture (an approximately 23,000-fold enrichment) of these genes allows for simplified allele calling. Tests on five genes (HLA-A/B/C/DRB1/DQB1) from 31 human samples and 351 datasets using STC-Seq showed results that were 98% consistent with the known two sets of digitals (field1 and field2) genotypes. Additionally, STC can capture genomic DNA fragments longer than 3 kb from HLA loci, making the library compatible with the third-generation sequencing.

Conclusions

STC-Seq is a highly accurate and cost-efficient method for HLA typing which can be used to facilitate the establishment of population-based HLA databases for the precision and transplantation medicine.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4431-5) contains supplementary material, which is available to authorized users.

Protocol for M3P: A Comprehensive and Clinical Oriented Targeted Sequencing Panel for Routine Molecular Analysis in Multiple Myeloma

Over the past 10 years next generation sequencing (NGS) approaches deciphered a large number of genomes from a wide variety of tumor types. However, despite most relevant findings, this technology has not yet been implemented into standard diagnostic workflows. Broad access to NGS technology is still limited, sequencing/analysis times exceed clinically relevant timeframes and despite huge cuts, costs remain significant. We proposed a custom-tailored gene panel, which focuses on a selected number of relevant genes and developed a clinically oriented NGS targeted sequencing approach for the molecular characterization of Multiple Myeloma (MM) tumors, allowing the description of the tumor genetic heterogeneity and its changes under selective pressure of antitumor therapy, in a more cost effective and faster turnaround timeframe.

Isotachophoresis applied to biomolecular reactions

This review discusses research developments and applications of isotachophoresis (ITP) to the initiation, control, and acceleration of chemical reactions, emphasizing reactions involving biomolecular reactants such as nucleic acids, proteins, and live cells. ITP is a versatile technique which requires no specific geometric design or material, and is compatible with a wide range of microfluidic and automated platforms. Though ITP has traditionally been used as a purification and separation technique, recent years have seen its emergence as a method to automate and speed up chemical reactions. ITP has been used to demonstrate up to 14 000-fold acceleration of nucleic acid assays, and has been used to enhance lateral flow and other immunoassays, and even whole bacterial cell detection assays. We here classify these studies into two categories: homogeneous (all reactants in solution) and heterogeneous (at least one reactant immobilized on a solid surface) assay configurations. For each category, we review and describe physical modeling and scaling of ITP-aided reaction assays, and elucidate key principles in ITP assay design. We summarize experimental advances, and identify common threads and approaches which researchers have used to optimize assay performance. Lastly, we propose unaddressed challenges and opportunities that could further improve these applications of ITP.

Solution-Phase Exome Capture

This protocol describes the construction of a paired-end library of genomic DNA (gDNA) and subsequent capture of specific regions of a genome using NimbleGen sequence capture probes and Illumina TruSeq oligos. The captured DNA, purified and quantitated, is appropriate for use as template in Illumina sequencing systems. A procedure is also provided for magnetic AMPure bead cleanup.

Liquid-phase sequence capture and targeted re-sequencing revealed novel polymorphisms in tomato genes belonging to the MEP carotenoid pathway

Tomato (Solanum lycopersicum L.) plants are characterized by having a variety of fruit colours that reflect the composition and accumulation of diverse carotenoids in the berries. Carotenoids are extensively studied for their health-promoting effects and this explains the great attention these pigments received by breeders and researchers worldwide. In this work we applied Agilent's SureSelect liquid-phase sequence capture and Illumina targeted re-sequencing of 34 tomato genes belonging to the methylerythritol phosphate (MEP) carotenoid pathway on a panel of 48 genotypes which differ for carotenoid content calculated as the sum of β-carotene, cis- and trans-lycopene. We targeted 230 kb of genomic regions including all exons and regulatory regions and observed ~40% of on-target capture. We found ample genetic variation among all the genotypes under study and generated an extensive catalog of SNPs/InDels located in both genic and regulatory regions. SNPs/InDels were also classified based on genomic location and putative biological effect. With our work we contributed to the identification of allelic variations possibly underpinning a key agronomic trait in tomato. Results from this study can be exploited for the promotion of novel studies on tomato bio-fortification as well as of breeding programs related to carotenoid accumulation in fruits.

ConcatSeq: A method for increasing throughput of single molecule sequencing by concatenating short DNA fragments

Single molecule sequencing (SMS) platforms enable base sequences to be read directly from individual strands of DNA in real-time. Though capable of long read lengths, SMS platforms currently suffer from low throughput compared to competing short-read sequencing technologies. Here, we present a novel strategy for sequencing library preparation, dubbed ConcatSeq, which increases the throughput of SMS platforms by generating long concatenated templates from pools of short DNA molecules. We demonstrate adaptation of this technique to two target enrichment workflows, commonly used for oncology applications, and feasibility using PacBio single molecule real-time (SMRT) technology. Our approach is capable of increasing the sequencing throughput of the PacBio RSII platform by more than five-fold, while maintaining the ability to correctly call allele frequencies of known single nucleotide variants. ConcatSeq provides a versatile new sample preparation tool for long-read sequencing technologies.

Targeted sequencing of both DNA strands barcoded and captured individually by RNA probes to identify genome-wide ultra-rare mutations

Next Generation Sequencing (NGS) has been widely implemented in biological research and has made a profound impact on patient care. One of the essential NGS applications is to identify disease-causing sequence variants, where high coverage and accuracy are needed. Here, we reported a novel NGS pipeline, termed a Sequencing System of Digitalized Barcode Encrypted Single-stranded Library from Extremely Low (quality and quantity) DNA Input with Probe-based DNA Enrichment by RNA probes targeting DNA duplex (DEEPER-Seq). This method combines an ultra-sensitive single-stranded library construction with barcoding error correction, termed DEEPER-Library; and a DNA capture approach using RNA probes targeting both DNA strands, termed DEEPER-Capture. DEEPER-Seq can create NGS libraries from as little as 20 pg DNA with PCR error correcting capabilities, and capture target sequences at an average ratio of 29.2% by targeting both DNA strands simultaneously with an over 98.6% coverage. Our method tags and sequences each of the two strands of a DNA duplex independently and only scores mutations that are found at the same position in both strands, which allows us to identify mutations with allelic fractions down to 0.03% in a whole exome sequencing (WES) study with a background error rate of one artificial error per 4.8 × 10⁹ nucleotides.

Exploiting the great potential of Sequence Capture data by a new tool, SUPER-CAP

The recent development of Sequence Capture methodology represents a powerful strategy for enhancing data generation to assess genetic variation of targeted genomic regions. Here, we present SUPER-CAP, a bioinformatics web tool aimed at handling Sequence Capture data, fine calculating the allele frequency of variations and building genotype-specific sequence of captured genes. The dataset used to develop this in silico strategy consists of 378 loci and related regulative regions in a collection of 44 tomato landraces. About 14,000 high-quality variants were identified. The high depth (>40×) of coverage and adopting the correct filtering criteria allowed identification of about 4,000 rare variants and 10 genes with a different copy number variation. We also show that the tool is capable to reconstruct genotype-specific sequences for each genotype by using the detected variants. This allows evaluating the combined effect of multiple variants in the same protein. The architecture and functionality of SUPER-CAP makes the software appropriate for a broad set of analyses including SNP discovery and mining. Its functionality, together with the capability to process large data sets and efficient detection of sequence variation, makes SUPER-CAP a valuable bioinformatics tool for genomics and breeding purposes.

Microfluidics for genome-wide studies involving next generation sequencing

Next-generation sequencing (NGS) has revolutionized how molecular biology studies are conducted. Its decreasing cost and increasing throughput permit profiling of genomic, transcriptomic, and epigenomic features for a wide range of applications. Microfluidics has been proven to be highly complementary to NGS technology with its unique capabilities for handling small volumes of samples and providing platforms for automation, integration, and multiplexing. In this article, we review recent progress on applying microfluidics to facilitate genome-wide studies. We emphasize on several technical aspects of NGS and how they benefit from coupling with microfluidic technology. We also summarize recent efforts on developing microfluidic technology for genomic, transcriptomic, and epigenomic studies, with emphasis on single cell analysis. We envision rapid growth in these directions, driven by the needs for testing scarce primary cell samples from patients in the context of precision medicine.

Targeted Capture of Complete Coding Regions across Divergent Species

Despite continued advances in sequencing technologies, there is a need for methods that can efficiently sequence large numbers of genes from diverse species. One approach to accomplish this is targeted capture (hybrid enrichment). While these methods are well established for genome resequencing projects, cross-species capture strategies are still being developed and generally focus on the capture of conserved regions, rather than complete coding regions from specific genes of interest. The resulting data is thus useful for phylogenetic studies, but the wealth of comparative data that could be used for evolutionary and functional studies is lost. Here, we design and implement a targeted capture method that enables recovery of complete coding regions across broad taxonomic scales. Capture probes were designed from multiple reference species and extensively tiled in order to facilitate cross-species capture. Using novel bioinformatics pipelines we were able to recover nearly all of the targeted genes with high completeness from species that were up to 200 myr divergent. Increased probe diversity and tiling for a subset of genes had a large positive effect on both recovery and completeness. The resulting data produced an accurate species tree, but importantly this same data can also be applied to studies of molecular evolution and function that will allow researchers to ask larger questions in broader phylogenetic contexts. Our method demonstrates the utility of cross-species approaches for the capture of full length coding sequences, and will substantially improve the ability for researchers to conduct large-scale comparative studies of molecular evolution and function.

From Genomes to GENE-omes: Exome Sequencing Concept and Applications in Crop Improvement

Exome sequencing represents targeted capture and sequencing of 1-2% of 'high-value genomic regions' (subset of the genome) which are enriched for functional variants and harbors low level of repetitive regions. We discuss here an overview of exome sequencing, ways to approach plant exomes, and advantages and applicability of this powerful approach in deciphering functional regions of genomes. Though initially this approach was developed as an alternative to whole genome sequencing (WGS), but the multitude of benefits conferred by sequence capture via hybridization approaches created a niche for itself to solve many of biological riddles, particularly for resolving phylogenetic distances. The technique has also proved to be successful in understanding the basis of natural and induced molecular variation, marker development and developing genomic resources for complex, wild and non-model species, which are still intractable for WGS efforts. Thus, with profound applications of this powerful sequencing strategy, near future is expected to witness a collective expansion of both techniques, i.e., sequence capture via hybridization for evolutionary and ecological research and WGS approaches for its universal accessibility.

Clinical Versus Research Sequencing

Historically, sequencing has been the key technology to assess variation in the genetic code, and has been widely accepted in clinical diagnostics of genetic disease. The advent of next-generation sequencing (NGS) methods increased the size of the analyzed target by several orders of magnitude, while at the same time drastically reducing the cost of sequencing. Current research allows sequencing of germline and tumor whole genomes. However, with the arrival of cutting-edge technology to the clinical diagnostic field, strict regulatory oversight is required to use the advances of the latest research when applied to routine clinical practice. We discuss the differences between sequencing in a research setting and sequencing in a clinical diagnostics setting, as applied to next-generation technology.

Capturing Darwin's dream

Evolutionary biologists from Darwin forward have dreamed of having data that would elucidate our understanding of evolutionary history and the diversity of life. Sequence capture is a relatively old DNA technology, but its use is growing rapidly due to advances in (i) massively parallel DNA sequencing approaches and instruments, (ii) massively parallel bait construction, (iii) methods to identify target regions and (iv) sample preparation. We give a little historical context to these developments, summarize some of the important advances reported in this special issue and point to further advances that can be made to help fulfill Darwin's dream.

Targeted sequencing for high-resolution evolutionary analyses following genome duplication in salmonid fish: Proof of concept for key components of the insulin-like growth factor axis

High-throughput sequencing has revolutionised comparative and evolutionary genome biology. It has now become relatively commonplace to generate multiple genomes and/or transcriptomes to characterize the evolution of large taxonomic groups of interest. Nevertheless, such efforts may be unsuited to some research questions or remain beyond the scope of some research groups. Here we show that targeted high-throughput sequencing offers a viable alternative to study genome evolution across a vertebrate family of great scientific interest. Specifically, we exploited sequence capture and Illumina sequencing to characterize the evolution of key components from the insulin-like growth (IGF) signalling axis of salmonid fish at unprecedented phylogenetic resolution. The IGF axis represents a central governor of vertebrate growth and its core components were expanded by whole genome duplication in the salmonid ancestor ~95Ma. Using RNA baits synthesised to genes encoding the complete family of IGF binding proteins (IGFBP) and an IGF hormone (IGF2), we captured, sequenced and assembled orthologous and paralogous exons from species representing all ten salmonid genera. This approach generated 299 novel sequences, most as complete or near-complete protein-coding sequences. Phylogenetic analyses confirmed congruent evolutionary histories for all nineteen recognized salmonid IGFBP family members and identified novel salmonid-specific IGF2 paralogues. Moreover, we reconstructed the evolution of duplicated IGF axis paralogues across a replete salmonid phylogeny, revealing complex historic selection regimes - both ancestral to salmonids and lineage-restricted - that frequently involved asymmetric paralogue divergence under positive and/or relaxed purifying selection. Our findings add to an emerging literature highlighting diverse applications for targeted sequencing in comparative-evolutionary genomics. We also set out a viable approach to obtain large sets of nuclear genes for any member of the salmonid family, which should enable insights into the evolutionary role of whole genome duplication before additional nuclear genome sequences become available.

Sequence capture by hybridization to explore modern and ancient genomic diversity in model and nonmodel organisms

The recent expansion of next-generation sequencing has significantly improved biological research. Nevertheless, deep exploration of genomes or metagenomic samples remains difficult because of the sequencing depth and the associated costs required. Therefore, different partitioning strategies have been developed to sequence informative subsets of studied genomes. Among these strategies, hybridization capture has proven to be an innovative and efficient tool for targeting and enriching specific biomarkers in complex DNA mixtures. It has been successfully applied in numerous areas of biology, such as exome resequencing for the identification of mutations underlying Mendelian or complex diseases and cancers, and its usefulness has been demonstrated in the agronomic field through the linking of genetic variants to agricultural phenotypic traits of interest. Moreover, hybridization capture has provided access to underexplored, but relevant fractions of genomes through its ability to enrich defined targets and their flanking regions. Finally, on the basis of restricted genomic information, this method has also allowed the expansion of knowledge of nonreference species and ancient genomes and provided a better understanding of metagenomic samples. In this review, we present the major advances and discoveries permitted by hybridization capture and highlight the potency of this approach in all areas of biology.

Optoelectrofluidic enhanced immunoreaction based on optically-induced dynamic AC electroosmosis

We report a novel optoelectrofluidic immunoreaction system based on electroosmotic flow for enhancing antibody-analyte binding efficiency on a surface-based sensing system. Two conventional indium tin oxide glass slides are assembled to provide a reaction chamber for a tiny volume of sample droplet (∼5 μL), in which the top layer is employed as an antibody-immobilized substrate and the bottom layer acts as a photoconductive layer of an optoelectrofluidic device. Under the application of an AC voltage, an illuminated light pattern on the photoconductive layer causes strong counter-rotating vortices to transport analytes from the bulk solution to the vicinity of the assay spot on the glass substrate. This configuration overcomes the slow immunoreaction problem of a diffusion-based sensing system, resulting in the enhancement of binding efficiency via an optoelectrofluidic method. Furthermore, we investigate the effect of optically-induced dynamic AC electroosmotic flow on optoelectrofluidic enhancement for surface-based immunoreaction with a mathematical simulation study and real experiments using immunoglobulin G (IgG) and anti-IgG. As a result, dynamic light patterns provided better immunoreaction efficiency than static light patterns due to effective mass transport of the target analyte, resulting in an achievement of 2.18-fold enhancement under a growing circular light pattern compared to the passive mode.

Advances in clinical next-generation sequencing: target enrichment and sequencing technologies

The huge parallel sequencing capabilities of next generation sequencing technologies have made them the tools of choice to characterize genomic aberrations for research and diagnostic purposes. For clinical applications, screening the whole genome or exome is challenging owing to the large genomic area to be sequenced, associated costs, complexity of data, and lack of known clinical significance of all genes. Consequently, routine screening involves limited markers with established clinical relevance. This process, referred to as targeted genome sequencing, requires selective enrichment of the genomic areas comprising these markers via one of several primer or probe-based enrichment strategies, followed by sequencing of the enriched genomic areas. Here, the authors review current target enrichment approaches and next generation sequencing platforms, focusing on the underlying principles, capabilities, and limitations of each technology along with validation and implementation for clinical testing.

DNA Methylation Analysis: Choosing the Right Method

In the burgeoning field of epigenetics, there are several methods available to determine the methylation status of DNA samples. However, choosing the method that is best suited to answering a particular biological question still proves to be a difficult task. This review aims to provide biologists, particularly those new to the field of epigenetics, with a simple algorithm to help guide them in the selection of the most appropriate assay to meet their research needs. First of all, we have separated all methods into two categories: those that are used for: (1) the discovery of unknown epigenetic changes; and (2) the assessment of DNA methylation within particular regulatory regions/genes of interest. The techniques are then scrutinized and ranked according to their robustness, high throughput capabilities and cost. This review includes the majority of methods available to date, but with a particular focus on commercially available kits or other simple and straightforward solutions that have proven to be useful.

Targeted capture in evolutionary and ecological genomics

The rapid expansion of next-generation sequencing has yielded a powerful array of tools to address fundamental biological questions at a scale that was inconceivable just a few years ago. Various genome-partitioning strategies to sequence select subsets of the genome have emerged as powerful alternatives to whole-genome sequencing in ecological and evolutionary genomic studies. High-throughput targeted capture is one such strategy that involves the parallel enrichment of preselected genomic regions of interest. The growing use of targeted capture demonstrates its potential power to address a range of research questions, yet these approaches have yet to expand broadly across laboratories focused on evolutionary and ecological genomics. In part, the use of targeted capture has been hindered by the logistics of capture design and implementation in species without established reference genomes. Here we aim to (i) increase the accessibility of targeted capture to researchers working in nonmodel taxa by discussing capture methods that circumvent the need of a reference genome, (ii) highlight the evolutionary and ecological applications where this approach is emerging as a powerful sequencing strategy and (iii) discuss the future of targeted capture and other genome-partitioning approaches in the light of the increasing accessibility of whole-genome sequencing. Given the practical advantages and increasing feasibility of high-throughput targeted capture, we anticipate an ongoing expansion of capture-based approaches in evolutionary and ecological research, synergistic with an expansion of whole-genome sequencing.

Targeted amplification for enhanced detection of biothreat agents by next-generation sequencing

BACKGROUND

Historically, identification of causal agents of disease has relied heavily on the ability to culture the organism in the laboratory and/or the use of pathogen-specific antibodies or sequence-based probes. However, these methods can be limiting: Even highly sensitive PCR-based assays must be continually updated due to signature degradation as new target strains and near neighbors are sequenced. Thus, there has been a need for assays that do not suffer as greatly from these limitations and/or biases. Recent advances in library preparation technologies for Next-Generation Sequencing (NGS) are focusing on the use of targeted amplification and targeted enrichment/capture to ensure that the most highly discriminating regions of the genomes of known targets (organism-unique regions and/or regions containing functionally important genes or phylogenetically-discriminating SNPs) will be sequenced, regardless of the complex sample background.

RESULTS

In the present study, we have assessed the feasibility of targeted sequence enhancement via amplification to facilitate detection of a bacterial pathogen present in low copy numbers in a background of human genomic material. Our results indicate that the targeted amplification of signature regions can effectively identify pathogen genomic material present in as little as 10 copies per ml in a complex sample. Importantly, the correct species and strain calls could be made in amplified samples, while this was not possible in unamplified samples.

CONCLUSIONS

The results presented here demonstrate the efficacy of a targeted amplification approach to biothreat detection, using multiple highly-discriminative amplicons per biothreat organism that provide redundancy in case of variation in some primer regions. Importantly, strain level discrimination was possible at levels of 10 genome equivalents. Similar results could be obtained through use of panels focused on the identification of amplicons targeted for specific genes or SNPs instead of, or in addition to, those targeted for specific organisms (ongoing gene-targeting work to be reported later). Note that without some form of targeted enhancement, the enormous background present in complex clinical and environmental samples makes it highly unlikely that sufficient coverage of key pathogen(s) present in the sample will be achieved with current NGS technology to guarantee that the most highly discriminating regions will be sequenced.

Enhanced virome sequencing using targeted sequence capture

Metagenomic shotgun sequencing (MSS) is an important tool for characterizing viral populations. It is culture independent, requires no a priori knowledge of the viruses in the sample, and may provide useful genomic information. However, MSS can lack sensitivity and may yield insufficient data for detailed analysis. We have created a targeted sequence capture panel, ViroCap, designed to enrich nucleic acid from DNA and RNA viruses from 34 families that infect vertebrate hosts. A computational approach condensed ∼1 billion bp of viral reference sequence into <200 million bp of unique, representative sequence suitable for targeted sequence capture. We compared the effectiveness of detecting viruses in standard MSS versus MSS following targeted sequence capture. First, we analyzed two sets of samples, one derived from samples submitted to a diagnostic virology laboratory and one derived from samples collected in a study of fever in children. We detected 14 and 18 viruses in the two sets, comprising 19 genera from 10 families, with dramatic enhancement of genome representation following capture enrichment. The median fold-increases in percentage viral reads post-capture were 674 and 296. Median breadth of coverage increased from 2.1% to 83.2% post-capture in the first set and from 2.0% to 75.6% in the second set. Next, we analyzed samples containing a set of diverse anellovirus sequences and demonstrated that ViroCap could be used to detect viral sequences with up to 58% variation from the references used to select capture probes. ViroCap substantially enhances MSS for a comprehensive set of viruses and has utility for research and clinical applications.

Common and Rare Variant Association Study for Plasma Lipids and Coronary Artery Disease

Blood lipid levels are highly heritable and modifiable risk factors for coronary artery disease (CAD), and are the leading cause of death worldwide. These facts have motivated human genetic association studies that have the substantial potential to define the risk factors that are causal and to identify pathways and therapeutic targets for lipids and CAD.The success of the HapMap project that provided an extensive catalog of human genetic variations and the development of microarray based genotyping chips (typically containing variations with allele frequencies ＞ 5%) facilitated common variant association study (CVAS; formerly termed genome-wide association study, GWAS) identifying disease-associated variants in a genome-wide manner. To date, 157 loci associated with blood lipids and 46 loci with CAD have been successfully identified, accounting for approximately 12%-14% of heritability for lipids and 10% of heritability for CAD. However, there is yet a major challenge termed "missing heritability problem," namely the observation that loci detected by CVAS explain only a small fraction of the inferred genetic variations. To explain such missing portions, focuses in genetic association studies have shifted from common to rare variants. However, it is challenging to apply rare variant association study (RVAS) in an unbiased manner because such variants typically lack the sufficient number to be identified statistically.In this review, we provide a current understanding of the genetic architecture mostly derived from CVAS, and several updates on the progress and limitations of RVAS for lipids and CAD.