High-Throughput Sequencing and Copy Number Variation Detection Using Formalin Fixed Embedded Tissue in Metastatic Gastric Cancer

In the era of targeted therapy, mutation profiling of cancer is a crucial aspect of making therapeutic decisions. To characterize cancer at a molecular level, the use of formalin-fixed paraffin-embedded tissue is important. We tested the Ion AmpliSeq Cancer Hotspot Panel v2 and nCounter Copy Number Variation Assay in 89 formalin-fixed paraffin-embedded gastric cancer samples to determine whether they are applicable in archival clinical samples for personalized targeted therapies. We validated the results with Sanger sequencing, real-time quantitative PCR, fluorescence in situ hybridization and immunohistochemistry. Frequently detected somatic mutations included TP53 (28.17%), APC (10.1%), PIK3CA (5.6%), KRAS (4.5%), SMO (3.4%), STK11 (3.4%), CDKN2A (3.4%) and SMAD4 (3.4%). Amplifications of HER2, CCNE1, MYC, KRAS and EGFR genes were observed in 8 (8.9%), 4 (4.5%), 2 (2.2%), 1 (1.1%) and 1 (1.1%) cases, respectively. In the cases with amplification, fluorescence in situ hybridization for HER2 verified gene amplification and immunohistochemistry for HER2, EGFR and CCNE1 verified the overexpression of proteins in tumor cells. In conclusion, we successfully performed semiconductor-based sequencing and nCounter copy number variation analyses in formalin-fixed paraffin-embedded gastric cancer samples. High-throughput screening in archival clinical samples enables faster, more accurate and cost-effective detection of hotspot mutations or amplification in genes.     

SLAF-seq: An Efficient Method of Large-Scale De Novo SNP Discovery and Genotyping Using High-Throughput Sequencing

Large-scale genotyping plays an important role in genetic association studies. It has provided new opportunities for gene discovery, especially when combined with high-throughput sequencing technologies. Here, we report an efficient solution for large-scale genotyping. We call it specific-locus amplified fragment sequencing (SLAF-seq). SLAF-seq technology has several distinguishing characteristics: i) deep sequencing to ensure genotyping accuracy; ii) reduced representation strategy to reduce sequencing costs; iii) pre-designed reduced representation scheme to optimize marker efficiency; and iv) double barcode system for large populations. In this study, we tested the efficiency of SLAF-seq on rice and soybean data. Both sets of results showed strong consistency between predicted and practical SLAFs and considerable genotyping accuracy. We also report the highest density genetic map yet created for any organism without a reference genome sequence, common carp in this case, using SLAF-seq data. We detected 50,530 high-quality SLAFs with 13,291 SNPs genotyped in 211 individual carp. The genetic map contained 5,885 markers with 0.68 cM intervals on average. A comparative genomics study between common carp genetic map and zebrafish genome sequence map showed high-quality SLAF-seq genotyping results. SLAF-seq provides a high-resolution strategy for large-scale genotyping and can be generally applicable to various species and populations.     

Direct Metagenomic Detection of Viral Pathogens in Nasal and Fecal Specimens Using an Unbiased High-Throughput Sequencing Approach

With the severe acute respiratory syndrome epidemic of 2003 and renewed attention on avian influenza viral pandemics, new surveillance systems are needed for the earlier detection of emerging infectious diseases. We applied a “next-generation” parallel sequencing platform for viral detection in nasopharyngeal and fecal samples collected during seasonal influenza virus (Flu) infections and norovirus outbreaks from 2005 to 2007 in Osaka, Japan. Random RT-PCR was performed to amplify RNA extracted from 0.1–0.25 ml of nasopharyngeal aspirates (N = 3) and fecal specimens (N = 5), and more than 10 µg of cDNA was synthesized. Unbiased high-throughput sequencing of these 8 samples yielded 15,298–32,335 (average 24,738) reads in a single 7.5 h run. In nasopharyngeal samples, although whole genome analysis was not available because the majority (>90%) of reads were host genome–derived, 20–460 Flu-reads were detected, which was sufficient for subtype identification. In fecal samples, bacteria and host cells were removed by centrifugation, resulting in gain of 484–15,260 reads of norovirus sequence (78–98% of the whole genome was covered), except for one specimen that was under-detectable by RT-PCR. These results suggest that our unbiased high-throughput sequencing approach is useful for directly detecting pathogenic viruses without advance genetic information. Although its cost and technological availability make it unlikely that this system will very soon be the diagnostic standard worldwide, this system could be useful for the earlier discovery of novel emerging viruses and bioterrorism, which are difficult to detect with conventional procedures.     

Identification of miRNAs and Their Target Genes in Peach (Prunus persica L.) Using High-Throughput Sequencing and Degradome Analysis

MicroRNAs play critical roles in various biological and metabolic processes. The function of miRNAs has been widely studied in model plants such as Arabidopsis and rice. However, the number of identified miRNAs and related miRNA targets in peach (Prunus persica) is limited. To understand further the relationship between miRNAs and their target genes during tissue development in peach, a small RNA library and three degradome libraries were constructed from three tissues for deep sequencing. We identified 117 conserved miRNAs and 186 novel miRNA candidates in peach by deep sequencing and 19 conserved miRNAs and 13 novel miRNAs were further evaluated for their expression by RT-qPCR. The number of gene targets that were identified for 26 conserved miRNA families and 38 novel miRNA candidates, were 172 and 87, respectively. Some of the identified miRNA targets were abundantly represented as conserved miRNA targets in plant. However, some of them were first identified and showed important roles in peach development. Our study provides information concerning the regulatory network of miRNAs in peach and advances our understanding of miRNA functions during tissue development.     

Ultra High Throughput Sequencing in Human DNA Variation Detection: A Comparative Study on the NDUFA3-PRPF31 Region

Background

Ultra high throughput sequencing (UHTS) technologies find an important application in targeted resequencing of candidate genes or of genomic intervals from genetic association studies. Despite the extraordinary power of these new methods, they are still rarely used in routine analysis of human genomic variants, in part because of the absence of specific standard procedures. The aim of this work is to provide human molecular geneticists with a tool to evaluate the best UHTS methodology for efficiently detecting DNA changes, from common SNPs to rare mutations.

Methodology/Principal Findings

We tested the three most widespread UHTS platforms (Roche/454 GS FLX Titanium, Illumina/Solexa Genome Analyzer II and Applied Biosystems/SOLiD System 3) on a well-studied region of the human genome containing many polymorphisms and a very rare heterozygous mutation located within an intronic repetitive DNA element. We identify the qualities and the limitations of each platform and describe some peculiarities of UHTS in resequencing projects.

Conclusions/Significance

When appropriate filtering and mapping procedures are applied UHTS technology can be safely and efficiently used as a tool for targeted human DNA variations detection. Unless particular and platform-dependent characteristics are needed for specific projects, the most relevant parameter to consider in mainstream human genome resequencing procedures is the cost per sequenced base-pair associated to each machine.     

Pilon: An Integrated Tool for Comprehensive Microbial Variant Detection and Genome Assembly Improvement

Advances in modern sequencing technologies allow us to generate sufficient data to analyze hundreds of bacterial genomes from a single machine in a single day. This potential for sequencing massive numbers of genomes calls for fully automated methods to produce high-quality assemblies and variant calls. We introduce Pilon, a fully automated, all-in-one tool for correcting draft assemblies and calling sequence variants of multiple sizes, including very large insertions and deletions. Pilon works with many types of sequence data, but is particularly strong when supplied with paired end data from two Illumina libraries with small e.g., 180 bp and large e.g., 3–5 Kb inserts. Pilon significantly improves draft genome assemblies by correcting bases, fixing mis-assemblies and filling gaps. For both haploid and diploid genomes, Pilon produces more contiguous genomes with fewer errors, enabling identification of more biologically relevant genes. Furthermore, Pilon identifies small variants with high accuracy as compared to state-of-the-art tools and is unique in its ability to accurately identify large sequence variants including duplications and resolve large insertions. Pilon is being used to improve the assemblies of thousands of new genomes and to identify variants from thousands of clinically relevant bacterial strains. Pilon is freely available as open source software.     

Detection and Typing of Human Rotavirus in Reference to Repeated Acute Gastroenteritis in Infants

Stool specimens from infants who visited a clinic because of acute gastroenteritis were tested for the presence of human rotavirus. Among the samples obtained were specimens taken from seven patients who had visited the clinic at two different times. In six of these seven children, human rotavirus (HRV) was detected in only one of the specimens taken (i.e. during only one of the two visits). One patient was shown to have excreted HRV twice; in both cases the HRV was serotyped to be type 1. The present results indicate that the symptomatic reinfection of HRV was not a widely occurring phenomenon in the group of infants tested.     

Combined Targeted DNA Sequencing in Non-Small Cell Lung Cancer (NSCLC) Using UNCseq and NGScopy,and RNA Sequencing Using UNCqeR for the Detection of Genetic Aberrations in NSCLC

The recent FDA approval of the MiSeqDx platform provides a unique opportunity to develop targeted next generation sequencing (NGS) panels for human disease, including cancer. We have developed a scalable, targeted panel-based assay termed UNCseq, which involves a NGS panel of over 200 cancer-associated genes and a standardized downstream bioinformatics pipeline for detection of single nucleotide variations (SNV) as well as small insertions and deletions (indel). In addition, we developed a novel algorithm, NGScopy, designed for samples with sparse sequencing coverage to detect large-scale copy number variations (CNV), similar to human SNP Array 6.0 as well as small-scale intragenic CNV. Overall, we applied this assay to 100 snap-frozen lung cancer specimens lacking same-patient germline DNA (07–0120 tissue cohort) and validated our results against Sanger sequencing, SNP Array, and our recently published integrated DNA-seq/RNA-seq assay, UNCqeR, where RNA-seq of same-patient tumor specimens confirmed SNV detected by DNA-seq, if RNA-seq coverage depth was adequate. In addition, we applied the UNCseq assay on an independent lung cancer tumor tissue collection with available same-patient germline DNA (11–1115 tissue cohort) and confirmed mutations using assays performed in a CLIA-certified laboratory. We conclude that UNCseq can identify SNV, indel, and CNV in tumor specimens lacking germline DNA in a cost-efficient fashion.     

Genus- and Species-Specific PCR-Based Detection of Dairy Propionibacteria in Environmental Samples by Using Primers Targeted to the Genes Encoding 16S rRNA

PCR assays with primers targeted to the genes encoding 16S rRNA were developed for detection of dairy propionibacteria. Propionibacterium thoenii specific oligonucleotide PT3 was selected after partial resequencing. Tests allowed the detection of less than 10 cells per reaction from milk and cheese and 10² cells per reaction from forage and soil.     

Unveiling of the Diversity of Prasinoviruses (Phycodnaviridae) in Marine Samples by Using High-Throughput Sequencing Analyses of PCR-Amplified DNA Polymerase and Major Capsid Protein Genes

Viruses strongly influence the ecology and evolution of their eukaryotic hosts in the marine environment, but little is known about their diversity and distribution. Prasinoviruses infect an abundant and widespread class of phytoplankton, the Mamiellophyceae, and thereby exert a specific and important role in microbial ecosystems. However, molecular tools to specifically identify this viral genus in environmental samples are still lacking. We developed two primer sets, designed for use with polymerase chain reactions and 454 pyrosequencing technologies, to target two conserved genes, encoding the DNA polymerase (PolB gene) and the major capsid protein (MCP gene). While only one copy of the PolB gene is present in Prasinovirus genomes, there are at least seven paralogs for MCP, the copy we named number 6 being shared with other eukaryotic alga-infecting viruses. Primer sets for PolB and MCP6 were thus designed and tested on 6 samples from the Tara Oceans project. The results suggest that the MCP6 amplicons show greater richness but that PolB gave a wider coverage of Prasinovirus diversity. As a consequence, we recommend use of the PolB primer set, which will certainly reveal exciting new insights about the diversity and distribution of prasinoviruses at the community scale.     

Training in High-Throughput Sequencing: Common Guidelines to Enable Material Sharing,Dissemination, and Reusability

The advancement of high-throughput sequencing (HTS) technologies and the rapid development of numerous analysis algorithms and pipelines in this field has resulted in an unprecedentedly high demand for training scientists in HTS data analysis. Embarking on developing new training materials is challenging for many reasons. Trainers often do not have prior experience in preparing or delivering such materials and struggle to keep them up to date. A repository of curated HTS training materials would support trainers in materials preparation, reduce the duplication of effort by increasing the usage of existing materials, and allow for the sharing of teaching experience among the HTS trainers’ community. To achieve this, we have developed a strategy for materials’ curation and dissemination. Standards for describing training materials have been proposed and applied to the curation of existing materials. A Git repository has been set up for sharing annotated materials that can now be reused, modified, or incorporated into new courses. This repository uses Git; hence, it is decentralized and self-managed by the community and can be forked/built-upon by all users. The repository is accessible at http://bioinformatics.upsc.se/htmr.     

Using High-Throughput Sequencing to Leverage Surveillance of Genetic Diversity and Oseltamivir Resistance: A Pilot Study during the 2009 Influenza A(H1N1) Pandemic

Background

Influenza viruses display a high mutation rate and complex evolutionary patterns. Next-generation sequencing (NGS) has been widely used for qualitative and semi-quantitative assessment of genetic diversity in complex biological samples. The “deep sequencing” approach, enabled by the enormous throughput of current NGS platforms, allows the identification of rare genetic viral variants in targeted genetic regions, but is usually limited to a small number of samples.

Methodology and Principal Findings

We designed a proof-of-principle study to test whether redistributing sequencing throughput from a high depth-small sample number towards a low depth-large sample number approach is feasible and contributes to influenza epidemiological surveillance. Using 454-Roche sequencing, we sequenced at a rather low depth, a 307 bp amplicon of the neuraminidase gene of the Influenza A(H1N1) pandemic (A(H1N1)pdm) virus from cDNA amplicons pooled in 48 barcoded libraries obtained from nasal swab samples of infected patients (n  =  299) taken from May to November, 2009 pandemic period in Mexico. This approach revealed that during the transition from the first (May-July) to second wave (September-November) of the pandemic, the initial genetic variants were replaced by the N248D mutation in the NA gene, and enabled the establishment of temporal and geographic associations with genetic diversity and the identification of mutations associated with oseltamivir resistance.

Conclusions

NGS sequencing of a short amplicon from the NA gene at low sequencing depth allowed genetic screening of a large number of samples, providing insights to viral genetic diversity dynamics and the identification of genetic variants associated with oseltamivir resistance. Further research is needed to explain the observed replacement of the genetic variants seen during the second wave. As sequencing throughput rises and library multiplexing and automation improves, we foresee that the approach presented here can be scaled up for global genetic surveillance of influenza and other infectious diseases.     

Unraveling Multiple MHC Gene Associations with Systemic Lupus Erythematosus: Model Choice Indicates a Role for HLA Alleles and Non-HLA Genes in Europeans

We have performed a meta-analysis of the major-histocompatibility-complex (MHC) region in systemic lupus erythematosus (SLE) to determine the association with both SNPs and classical human-leukocyte-antigen (HLA) alleles. More specifically, we combined results from six studies and well-known out-of-study control data sets, providing us with 3,701 independent SLE cases and 12,110 independent controls of European ancestry. This study used genotypes for 7,199 SNPs within the MHC region and for classical HLA alleles (typed and imputed). Our results from conditional analysis and model choice with the use of the Bayesian information criterion show that the best model for SLE association includes both classical loci (HLA-DRB1^∗03:01, HLA-DRB1^∗08:01, and HLA-DQA1^∗01:02) and two SNPs, rs8192591 (in class III and upstream of NOTCH4) and rs2246618 (MICB in class I). Our approach was to perform a stepwise search from multiple baseline models deduced from a priori evidence on HLA-DRB1 lupus-associated alleles, a stepwise regression on SNPs alone, and a stepwise regression on HLA alleles. With this approach, we were able to identify a model that was an overwhelmingly better fit to the data than one identified by simple stepwise regression either on SNPs alone (Bayes factor [BF] > 50) or on classical HLA alleles alone (BF > 1,000).     

RADtyping: An Integrated Package for Accurate De Novo Codominant and Dominant RAD Genotyping in Mapping Populations

Genetic linkage maps are indispensable tools in genetic, genomic and breeding studies. As one of genotyping-by-sequencing methods, RAD-Seq (restriction-site associated DNA sequencing) has gained particular popularity for construction of high-density linkage maps. Current RAD analytical tools are being predominantly used for typing codominant markers. However, no genotyping algorithm has been developed for dominant markers (resulting from recognition site disruption). Given their abundance in eukaryotic genomes, utilization of dominant markers would greatly diminish the extensive sequencing effort required for large-scale marker development. In this study, we established, for the first time, a novel statistical framework for de novo dominant genotyping in mapping populations. An integrated package called RADtyping was developed by incorporating both de novo codominant and dominant genotyping algorithms. We demonstrated the superb performance of RADtyping in achieving remarkably high genotyping accuracy based on simulated and real mapping datasets. The RADtyping package is freely available at http://www2.ouc.edu.cn/mollusk/ detailen.asp?id=727.     

Sequencing Biological and Physical Events Affects Specific Frequency Bands within the Human Premotor Cortex: An Intracerebral EEG Study

Evidence that the human premotor cortex (PMC) is activated by cognitive functions involving the motor domain is classically explained as the reactivation of a motor program decoupled from its executive functions, and exploited for different purposes by means of a motor simulation. In contrast, the evidence that PMC contributes to the sequencing of non-biological events cannot be explained by the simulationist theory. Here we investigated how motor simulation and event sequencing coexist within the PMC and how these mechanisms interact when both functions are executed. We asked patients with depth electrodes implanted in the PMC to passively observe a randomized arrangement of images depicting biological actions and physical events and, in a second block, to sequence them in the correct order. This task allowed us to disambiguate between the simple observation of actions, their sequencing (recruiting different motor simulation processes), as well as the sequencing of non-biological events (recruiting a sequencer mechanism non dependant on motor simulation). We analysed the response of the gamma, alpha and beta frequency bands to evaluate the contribution of each brain rhythm to the observation and sequencing of both biological and non-biological stimuli. We found that motor simulation (biological>physical) and event sequencing (sequencing>observation) differently affect the three investigated frequency bands: motor simulation was reflected on the gamma and, partially, in the beta, but not in the alpha band. In contrast, event sequencing was also reflected on the alpha band.     

Campylobacteriosis in Urban versus Rural Areas: A Case-Case Study Integrated with Molecular Typing to Validate Risk Factors and to Attribute Sources of Infection

Campylobacter infection is a leading cause of bacterial gastroenteritis worldwide, and most clinical cases appear as isolated, sporadic infections for which the source is rarely apparent. From July 2005 to December 2007 we conducted a prospective case-case study of sporadic, domestically-acquired Campylobacter enteritis in rural versus urban areas and a prevalence study of Campylobacter in animal and environmental sources in the Eastern Townships, Quebec. Isolates were typed using Multilocus Sequence Typing (MLST) to reinforce the case-case findings and to assign a source probability estimate for each human isolate. The risk of human campylobacteriosis was 1.89-fold higher in rural than urban areas. Unconditional multivariate logistic regression analysis identified two independent risk factors associated with human Campylobacter infections acquired in rural area: occupational exposure to animals (OR = 10.6, 95% CI: 1.2–91, p = 0.032), and household water coming from a private well (OR = 8.3, 95% CI: 3.4–20.4, p<0.0001). A total of 851 C. jejuni isolates (178 human, 257 chicken, 87 bovine, 266 water, 63 wild bird) were typed using MLST. Among human isolates, the incidence rates of clonal complexes (CC) CC-21, CC-45, and CC-61 were higher in rural than urban areas. MLST-based source attribution analysis indicated that 64.5% of human C. jejuni isolates were attributable to chicken, followed by cattle (25.8%), water (7.4%), and wild birds (2.3%). Chicken was the attributable source for the majority of cases, independent of residential area, sex and age. The increased incidence in rural compared to urban areas was associated with occupational exposure to animals, particularly cattle among those aged 15–34 years, and with consumption of private well water. Both bovine and water exposure appeared to contribute to the seasonal variation in campylobacteriosis. These results provide a basis for developing public education and preventive programs targeting the risk factors identified.