Gene set enrichment meta-learning analysis: next- generation sequencing versus microarrays

Background  

Reproducibility of results can have a significant impact on the acceptance of new technologies in gene expression analysis. With the recent introduction of the so-called next-generation sequencing (NGS) technology and established microarrays, one is able to choose between two completely different platforms for gene expression measurements. This study introduces a novel methodology for gene-ranking stability analysis that is applied to the evaluation of gene-ranking reproducibility on NGS and microarray data.     

Combinatorial analysis and algorithms for quasispecies reconstruction using next-generation sequencing

Background  

Next-generation sequencing (NGS) offers a unique opportunity for high-throughput genomics and has potential to replace Sanger sequencing in many fields, including de-novo sequencing, re-sequencing, meta-genomics, and characterisation of infectious pathogens, such as viral quasispecies. Although methodologies and software for whole genome assembly and genome variation analysis have been developed and refined for NGS data, reconstructing a viral quasispecies using NGS data remains a challenge. This application would be useful for analysing intra-host evolutionary pathways in relation to immune responses and antiretroviral therapy exposures. Here we introduce a set of formulae for the combinatorial analysis of a quasispecies, given a NGS re-sequencing experiment and an algorithm for quasispecies reconstruction. We require that sequenced fragments are aligned against a reference genome, and that the reference genome is partitioned into a set of sliding windows (amplicons). The reconstruction algorithm is based on combinations of multinomial distributions and is designed to minimise the reconstruction of false variants, called in-silico recombinants.     

SAMMate: a GUI tool for processing short read alignments in SAM/BAM format

Background  

Next Generation Sequencing (NGS) technology generates tens of millions of short reads for each DNA/RNA sample. A key step in NGS data analysis is the short read alignment of the generated sequences to a reference genome. Although storing alignment information in the Sequence Alignment/Map (SAM) or Binary SAM (BAM) format is now standard, biomedical researchers still have difficulty accessing this information.     

UreB immunodominant epitope-specific CD8+ T-cell responses were beneficial in reducing gastric symptoms in Helicobacter pylori-infected individuals

Background and aims

Although Helicobacter pylori is recognized as an extracellular infection bacterium, it can lead to an increase in the number of CD8⁺ T cells after infection. At present, the characteristics of H. pylori antigen-specific CD8⁺ T cells and the epitope response have not been elucidated. This study was focused on putative protective antigen UreB to detect specific CD8⁺ T-cell responses in vitro and screen for predominant response epitopes.

Methods

The PBMCs collected from H. pylori-infected individuals were stimulated by UreB peptide pools in vitro to identify the immunodominant CD8⁺ T-cell epitopes. Furthermore, their HLA restriction characteristics were detected accordingly by NGS. Finally, the relationship between immunodominant responses and appearance of gastric symptoms after H. pylori infection was conducted.

Results

UreB-specific CD8⁺ T-cell responses were detected in H. pylori-infected individuals. Three of UreB dominant epitopes (A-2 (UreB_443–451: GVKPNMIIK), B-4 (UreB_420–428: SEYVGSVEV), and C-1 (UreB_5–13: SRKEYVSMY)) were firstly identified and mainly presented by HLA-A*1101, HLA-B*4001 and HLA-C*0702 alleles, respectively. C-1 responses were mostly occurred in H. pylori-infected subjects without gastric symptoms and may alleviate the degree of gastric inflammation.

Conclusions

The UreB dominant epitope-specific CD8⁺ T-cell response was closely related to the gastric symptoms after H. pylori infection, and the C-1 (UreB_5-13) dominant peptides may be protective epitopes.     

Targeted next generation sequencing identifies novel NOTCH3 gene mutations in CADASIL diagnostics patients

Background

Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is a monogenic, hereditary, small vessel disease of the brain causing stroke and vascular dementia in adults. CADASIL has previously been shown to be caused by varying mutations in the NOTCH3 gene. The disorder is often misdiagnosed due to its significant clinical heterogeneic manifestation with familial hemiplegic migraine and several ataxia disorders as well as the location of the currently identified causative mutations. The aim of this study was to develop a new, comprehensive and efficient single assay strategy for complete molecular diagnosis of NOTCH3 mutations through the use of a custom next-generation sequencing (NGS) panel for improved routine clinical molecular diagnostic testing.

Results

Our custom NGS panel identified nine genetic variants in NOTCH3 (p.D139V, p.C183R, p.R332C, p.Y465C, p.C597W, p.R607H, p.E813E, p.C977G and p.Y1106C). Six mutations were stereotypical CADASIL mutations leading to an odd number of cysteine residues in one of the 34 NOTCH3 gene epidermal growth factor (EGF)-like repeats, including three new typical cysteine mutations identified in exon 11 (p.C597W; c.1791C>G); exon 18 (p.C977G; c.2929T>G) and exon 20 (p.Y1106C; c.3317A>G). Interestingly, a novel missense mutation in the CACNA1A gene was also identified in one CADASIL patient. All variants identified (novel and known) were further investigated using in silico bioinformatic analyses and confirmed through Sanger sequencing.

Conclusions

NGS provides an improved and effective methodology for the diagnosis of CADASIL. The NGS approach reduced time and cost for comprehensive genetic diagnosis, placing genetic diagnostic testing within reach of more patients.

     

Alterations in cell growth and signaling in ErbB3 binding protein-1 (Ebp1) deficient mice

Background  

The ErbB3 binding protein-1 (Ebp1) belongs to a family of DNA/RNA binding proteins implicated in cell growth, apoptosis and differentiation. However, the physiological role of Ebp1 in the whole organism is not known. Therefore, we generated Ebp1-deficient mice carrying a gene trap insertion in intron 2 of the Ebp1 (pa2g4) gene.     

Molecular and functional analyses of COPT/Ctr-type copper transporter-like gene family in rice

Background  

The copper (Cu) transporter (COPT/Ctr) gene family has an important role in the maintenance of Cu homeostasis in different species. The rice COPT-type gene family consists of seven members (COPT1 to COPT7). However, only two, COPT1 and COPT5, have been characterized for their functions in Cu transport.     

Over expression of the selectable marker blasticidin S deaminase gene is toxic to human keratinocytes and murine BALB/MK cells

Background  

The blasticidin S resistance gene (bsr) is a selectable marker used for gene transfer experiments. The bsr gene encodes for blasticidin S (BS) deaminase, which has a specific activity upon BS. Therefore, its expression is supposed to be harmless in cells. The work reported on herein consisted of experiments to verify a possible toxicity of bsr on mammalian cells, which include several cell lines and primary cultures.     

Detection of Helicobacter pylori in stool samples of young children using real‐time polymerase chain reaction

Background

The aims of this study were to develop and validate a multiplex real‐time polymerase chain reaction (q‐PCR) assay of Helicobacter pylori in stool samples of healthy children. Additionally, we determined the prevalence of clarithromycin resistance and cagA gene in H. pylori‐positive samples.

Materials and methods

Archived stool samples from 188 children aged 6‐9 years and 272 samples of 92 infants aged 2‐18 months were tested for H. pylori antigens using enzyme immunoassay (EIA). A multiplex q‐PCR assay was designed to detect H. pylori 16S rRNA and urease and the human RNase P gene as an internal control. Kappa coefficient was calculated to assess the agreement between q‐PCR and EIA.

Results

Laboratory validation of the q‐PCR assay using quantitated H. pylori ATCC 43504 extracted DNA showed S‐shaped amplification curves for all genes; the limit of detection was 1 CFU/reaction. No cross‐reactivity with other bacterial pathogens was noted. Applying the multiplex q‐PCR to DNA extracted from fecal samples showed clear amplification curves for urease gene, but not for 16S rRNA. The prevalence of H. pylori infection was 50% (95% CI 43%‐57%) by q‐PCR (urease cycle threshold <44) vs 59% (95% CI 52%‐66%) by EIA. Kappa coefficient was .80 (P < .001) and .44 (P < .001) for children aged 6‐9 years and 2‐18 months, respectively. Sixteen samples were positive for cagA and three were positive for clarithromycin resistance mutation (A2143G) as confirmed by sequencing.

Conclusions

The developed q‐PCR can be used as a cotechnique to enhance the accuracy of H. pylori detection in epidemiological studies and in clinical settings.     

Assessment of FAE1 polymorphisms in three Brassica species using EcoTILLING and their association with differences in seed erucic acid contents

Background  

FAE1 (fatty acid elongase1) is the key gene in the control of erucic acid synthesis in seeds of Brassica species. Due to oil with low erucic acid (LEA) content is essential for human health and not enough LEA resource could be available, thus new LEA genetic resources are being sought for Brassica breeding. EcoTILLING, a powerful genotyping method, can readily be used to identify polymorphisms in Brassica.     

PileLine: a toolbox to handle genome position information in next-generation sequencing studies

Background  

Genomic position (GP) files currently used in next-generation sequencing (NGS) studies are always difficult to manipulate due to their huge size and the lack of appropriate tools to properly manage them. The structure of these flat files is based on representing one line per position that has been covered by at least one aligned read, imposing significant restrictions from a computational performance perspective.     

A fault-tolerant method for HLA typing with PacBio data

Background

Human leukocyte antigen (HLA) genes are critical genes involved in important biomedical aspects, including organ transplantation, autoimmune diseases and infectious diseases. The gene family contains the most polymorphic genes in humans and the difference between two alleles is only a single base pair substitution in many cases. The next generation sequencing (NGS) technologies could be used for high throughput HLA typing but in silico methods are still needed to correctly assign the alleles of a sample. Computer scientists have developed such methods for various NGS platforms, such as Illumina, Roche 454 and Ion Torrent, based on the characteristics of the reads they generate. However, the method for PacBio reads was less addressed, probably owing to its high error rates. The PacBio system has the longest read length among available NGS platforms, and therefore is the only platform capable of having exon 2 and exon 3 of HLA genes on the same read to unequivocally solve the ambiguity problem caused by the “phasing” issue.

Results

We proposed a new method BayesTyping1 to assign HLA alleles for PacBio circular consensus sequencing reads using Bayes’ theorem. The method was applied to simulated data of the three loci HLA-A, HLA-B and HLA-DRB1. The experimental results showed its capability to tolerate the disturbance of sequencing errors and external noise reads.

Conclusions

The BayesTyping1 method could overcome the problems of HLA typing using PacBio reads, which mostly arise from sequencing errors of PacBio reads and the divergence of HLA genes, to some extent.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-296) contains supplementary material, which is available to authorized users.     

Development and application of a next-generation-sequencing (NGS) approach to detect known and novel gene defects underlying retinal diseases

Background

Inherited retinal disorders are clinically and genetically heterogeneous with more than 150 gene defects accounting for the diversity of disease phenotypes. So far, mutation detection was mainly performed by APEX technology and direct Sanger sequencing of known genes. However, these methods are time consuming, expensive and unable to provide a result if the patient carries a new gene mutation. In addition, multiplicity of phenotypes associated with the same gene defect may be overlooked.

Methods

To overcome these challenges, we designed an exon sequencing array to target 254 known and candidate genes using Agilent capture. Subsequently, 20 DNA samples from 17 different families, including four patients with known mutations were sequenced using Illumina Genome Analyzer IIx next-generation-sequencing (NGS) platform. Different filtering approaches were applied to identify the genetic defect. The most likely disease causing variants were analyzed by Sanger sequencing. Co-segregation and sequencing analysis of control samples validated the pathogenicity of the observed variants.

Results

The phenotype of the patients included retinitis pigmentosa, congenital stationary night blindness, Best disease, early-onset cone dystrophy and Stargardt disease. In three of four control samples with known genotypes NGS detected the expected mutations. Three known and five novel mutations were identified in NR2E3, PRPF3, EYS, PRPF8, CRB1, TRPM1 and CACNA1F. One of the control samples with a known genotype belongs to a family with two clinical phenotypes (Best and CSNB), where a novel mutation was identified for CSNB. In six families the disease associated mutations were not found, indicating that novel gene defects remain to be identified.

Conclusions

In summary, this unbiased and time-efficient NGS approach allowed mutation detection in 75% of control cases and in 57% of test cases. Furthermore, it has the possibility of associating known gene defects with novel phenotypes and mode of inheritance.     

Characterization of a novel non-specific nuclease from thermophilic bacteriophage GBSV1

Background  

Thermostable enzymes from thermophiles have attracted extensive studies. In this investigation, a nuclease-encoding gene (designated as GBSV1-NSN) was obtained from a thermophilic bacteriophage GBSV1 for the first time.     

Characterization of a heat resistant ß-glucosidase as a new reporter in cells and mice

Background   

Reporter genes are widely used in biology and only a limited number are available. We present a new reporter gene for the localization of mammalian cells and transgenic tissues based on detection of the bglA (SYNbglA) gene of Caldocellum saccharolyticum that encodes a thermophilic β-glucosidase.     

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.     

Targeted Next Generation Sequencing as a Reliable Diagnostic Assay for the Detection of Somatic Mutations in Tumours Using Minimal DNA Amounts from Formalin Fixed Paraffin Embedded Material

Background

Targeted Next Generation Sequencing (NGS) offers a way to implement testing of multiple genetic aberrations in diagnostic pathology practice, which is necessary for personalized cancer treatment. However, no standards regarding input material have been defined. This study therefore aimed to determine the effect of the type of input material (e.g. formalin fixed paraffin embedded (FFPE) versus fresh frozen (FF) tissue) on NGS derived results. Moreover, this study aimed to explore a standardized analysis pipeline to support consistent clinical decision-making.

Method

We used the Ion Torrent PGM sequencing platform in combination with the Ion AmpliSeq Cancer Hotspot Panel v2 to sequence frequently mutated regions in 50 cancer related genes, and validated the NGS detected variants in 250 FFPE samples using standard diagnostic assays. Next, 386 tumour samples were sequenced to explore the effect of input material on variant detection variables. For variant calling, Ion Torrent analysis software was supplemented with additional variant annotation and filtering.

Results

Both FFPE and FF tissue could be sequenced reliably with a sensitivity of 99.1%. Validation showed a 98.5% concordance between NGS and conventional sequencing techniques, where NGS provided both the advantage of low input DNA concentration and the detection of low-frequency variants. The reliability of mutation analysis could be further improved with manual inspection of sequence data.

Conclusion

Targeted NGS can be reliably implemented in cancer diagnostics using both FFPE and FF tissue when using appropriate analysis settings, even with low input DNA.