Christophe Van Neste

Forensic STR analysis using massive parallel sequencing

We explore the applicability of second generation sequencing (SGS) to sequence multiplexed forensic STR amplicons, both in a single contributor sample as in multiple-person mixtures with different ratios. We compare the results of a commercial STR profiling kit (Applied Biosystems AmpFlSTR(®) Profiler Plus(®)), analyzed both with capillary electrophoresis and with Roche GS FLX sequencing. An easy to use open-source software pipeline is provided, chaining together the different steps needed to start the analysis from a GS FLX FASTA file, resulting in a FASTA file containing the called and quantified alleles present in the data. Sequencing of multiplexed STR amplicons using Roche GS FLX titanium technology is technically feasible but the technology is not ideal for this purpose. The fraction of full length reads is small and the homopolymer sequencing error rate is high. The pipeline compresses the homopolymers to a single base to avoid false results caused by these homopolymers. The qualitative and quantitative results from the SGS STR analysis pipeline are comparable to the electrophoresis method. Additionally, the SGS method provides extra information and is able to call allele subtypes based on STR sequences in a database. In mixed samples, all alleles were reported from individuals that contributed at least 10% to the mixture.

Massively parallel sequencing of forensic STRs: Considerations of the DNA commission of the International Society for Forensic Genetics (ISFG) on minimal nomenclature requirements

The DNA Commission of the International Society for Forensic Genetics (ISFG) is reviewing factors that need to be considered ahead of the adoption by the forensic community of short tandem repeat (STR) genotyping by massively parallel sequencing (MPS) technologies. MPS produces sequence data that provide a precise description of the repeat allele structure of a STR marker and variants that may reside in the flanking areas of the repeat region. When a STR contains a complex arrangement of repeat motifs, the level of genetic polymorphism revealed by the sequence data can increase substantially. As repeat structures can be complex and include substitutions, insertions, deletions, variable tandem repeat arrangements of multiple nucleotide motifs, and flanking region SNPs, established capillary electrophoresis (CE) allele descriptions must be supplemented by a new system of STR allele nomenclature, which retains backward compatibility with the CE data that currently populate national DNA databases and that will continue to be produced for the coming years. Thus, there is a pressing need to produce a standardized framework for describing complex sequences that enable comparison with currently used repeat allele nomenclature derived from conventional CE systems. It is important to discern three levels of information in hierarchical order (i) the sequence, (ii) the alignment, and (iii) the nomenclature of STR sequence data. We propose a sequence (text) string format the minimal requirement of data storage that laboratories should follow when adopting MPS of STRs. We further discuss the variant annotation and sequence comparison framework necessary to maintain compatibility among established and future data. This system must be easy to use and interpret by the DNA specialist, based on a universally accessible genome assembly, and in place before the uptake of MPS by the general forensic community starts to generate sequence data on a large scale. While the established nomenclature for CE-based STR analysis will remain unchanged in the future, the nomenclature of sequence-based STR genotypes will need to follow updated rules and be generated by expert systems that translate MPS sequences to match CE conventions in order to guarantee compatibility between the different generations of STR data.

My-Forensic-Loci-queries (MyFLq) framework for analysis of forensic STR data generated by massive parallel sequencing

Forensic scientists are currently investigating how to transition from capillary electrophoresis (CE) to massive parallel sequencing (MPS) for analysis of forensic DNA profiles. MPS offers several advantages over CE such as virtually unlimited multiplexy of loci, combining both short tandem repeat (STR) and single nucleotide polymorphism (SNP) loci, small amplicons without constraints of size separation, more discrimination power, deep mixture resolution and sample multiplexing. We present our bioinformatic framework My-Forensic-Loci-queries (MyFLq) for analysis of MPS forensic data. For allele calling, the framework uses a MySQL reference allele database with automatically determined regions of interest (ROIs) by a generic maximal flanking algorithm which makes it possible to use any STR or SNP forensic locus. Python scripts were designed to automatically make allele calls starting from raw MPS data. We also present a method to assess the usefulness and overall performance of a forensic locus with respect to MPS, as well as methods to estimate whether an unknown allele, which sequence is not present in the MySQL database, is in fact a new allele or a sequencing error. The MyFLq framework was applied to an Illumina MiSeq dataset of a forensic Illumina amplicon library, generated from multilocus STR polymerase chain reaction (PCR) on both single contributor samples and multiple person DNA mixtures. Although the multilocus PCR was not yet optimized for MPS in terms of amplicon length or locus selection, the results show excellent results for most loci. The results show a high signal-to-noise ratio, correct allele calls, and a low limit of detection for minor DNA contributors in mixed DNA samples. Technically, forensic MPS affords great promise for routine implementation in forensic genomics. The method is also applicable to adjacent disciplines such as molecular autopsy in legal medicine and in mitochondrial DNA research.

Differential transcriptome analysis of glandular and filamentous trichomes in Artemisia annua.

BackgroundThe medicinal plant Artemisia annua is covered with filamentous trichomes and glandular, artemisinin producing trichomes. A high artemisinin supply is needed at a reduced cost for treating malaria. Artemisinin production in bioreactors can be facilitated if a better insight is obtained in the biosynthesis of artemisinin and other metabolites. Therefore, metabolic activities of glandular and filamentous trichomes were investigated at the transcriptome level.ResultsBy laser pressure catapulting, glandular and filamentous trichomes as well as apical and sub-apical cells from glandular trichomes were collected and their transcriptome was sequenced using Illumina RNA-Seq. A de novo transcriptome was assembled (Trinity) and studied with a differential expression analysis (edgeR).A comparison of the transcriptome from glandular and filamentous trichomes shows that MEP, MVA, most terpene and lipid biosynthesis pathways are significantly upregulated in glandular trichomes. Conversely, some transcripts coding for specific sesquiterpenoid and triterpenoid enzymes such as 8-epi-cedrol synthase and an uncharacterized oxidosqualene cyclase were significantly upregulated in filamentous trichomes. All known artemisinin biosynthesis genes are upregulated in glandular trichomes and were detected in both the apical and sub-apical cells of the glandular trichomes. No significant differential expression could be observed between the apical and sub-apical cells.ConclusionsOur results underscore the vast metabolic capacities of A. annua glandular trichomes but nonetheless point to the existence of specific terpene metabolic pathways in the filamentous trichomes. Candidate genes that might be involved in artemisinin biosynthesis are proposed based on their putative function and their differential expression level.

TRAPID: an efficient online tool for the functional and comparative analysis of de novo RNA-Seq transcriptomes

Transcriptome analysis through next-generation sequencing technologies allows the generation of detailed gene catalogs for non-model species, at the cost of new challenges with regards to computational requirements and bioinformatics expertise. Here, we present TRAPID, an online tool for the fast and efficient processing of assembled RNA-Seq transcriptome data, developed to mitigate these challenges. TRAPID offers high-throughput open reading frame detection, frameshift correction and includes a functional, comparative and phylogenetic toolbox, making use of 175 reference proteomes. Benchmarking and comparison against state-of-the-art transcript analysis tools reveals the efficiency and unique features of the TRAPID system. TRAPID is freely available at http://bioinformatics.psb.ugent.be/webtools/trapid/.

Forensic massively parallel sequencing data analysis tool: Implementation of MyFLq as a standalone web- and Illumina BaseSpace(®)-application.

Routine use of massively parallel sequencing (MPS) for forensic genomics is on the horizon. The last few years, several algorithms and workflows have been developed to analyze forensic MPS data. However, none have yet been tailored to the needs of the forensic analyst who does not possess an extensive bioinformatics background. We developed our previously published forensic MPS data analysis framework MyFLq (My-Forensic-Loci-queries) into an open-source, user-friendly, web-based application. It can be installed as a standalone web application, or run directly from the Illumina BaseSpace environment. In the former, laboratories can keep their data on-site, while in the latter, data from forensic samples that are sequenced on an Illumina sequencer can be uploaded to Basespace during acquisition, and can subsequently be analyzed using the published MyFLq BaseSpace application. Additional features were implemented such as an interactive graphical report of the results, an interactive threshold selection bar, and an allele length-based analysis in addition to the sequenced-based analysis. Practical use of the application is demonstrated through the analysis of four 16-plex short tandem repeat (STR) samples, showing the complementarity between the sequence- and length-based analysis of the same MPS data.

Reference loci for RT-qPCR analysis of differentiating human embryonic stem cells

BackgroundSelecting stably expressed reference genes is essential for proper reverse transcription quantitative polymerase chain reaction gene expression analysis. However, this choice is not always straightforward. In the case of differentiating human embryonic stem (hES) cells, differentiation itself introduces changes whereby reference gene stability may be influenced.ResultsIn this study, we evaluated the stability of various references during retinoic acid-induced (2 microM) differentiation of hES cells. Out of 12 candidate references, beta-2-microglobulin, ribosomal protein L13A and Alu repeats are found to be the most stable for this experimental set-up.ConclusionsOur results show that some of the commonly used reference genes are actually not amongst the most stable loci during hES cell differentiation promoted by retinoic acid. Moreover, a novel normalization strategy based on expressed Alu repeats is validated for use in hES cell experiments.

Differential transcriptome analysis of the common shrimp Crangon crangon: Special focus on the nuclear receptors and RNAi-related genes

The decapod Crangon crangon is one of the most valuable European fisheries commodities. Despite its economic importance, little sequence data is available for this shrimp species. In this paper, we report the transcriptome sequencing for five different stages of C. crangon (early embryo, late embryo, larva, female adults and male adults) and the annotation and stage-specific expression analysis of nuclear receptors (NRs) and RNA interference (RNAi)-related genes. The NRs are transcription factors that play an essential role in growth, development, cell differentiation, molting/metamorphosis and reproduction, while the RNAi-related genes are very important for internal gene expression regulation and in antiviral defense. We discovered a NR in the female C. crangon which is either a very rapidly evolved homolog of HR10, or a novel NR altogether. This new NR could act as a biological marker for sex determination as it is not expressed in male adults. Most RNAi-related genes were present in C. crangon, proving that the requirements for successful RNAi is present in this decapod shrimp. RNAi-based applications in Crangon such as its use in functional genomics or as antiviral therapeutics could become very important in the near future.

Dynamic epigenetic changes to VHL occur with sunitinib in metastatic clear cell renal cancer

Background Genetic intratumoral heterogeneity (ITH) hinders biomarker development in metastatic clear cell renal cancer (mccRCC). Epigenetic relative to genetic ITH or the presence of consistent epigenetic changes following targeted therapy in mccRCC have not been evaluated. The aim of this study was to determine methylome/genetic ITH and to evaluate specific epigenetic and genetic changes associated with sunitinib therapy. Patients and methods Multi-region DNA sampling performed on sequential frozen pairs of primary tumor tissue from 14 metastatic ccRCC patients, in the Upfront Sunitinib (SU011248) Therapy Followed by Surgery in Patients with Metastatic Renal Cancer: a Pilot Phase II Study (SuMR; ClinicalTrials.gov identifier: NCT01024205), at presentation (biopsy) and after 3-cycles of 50mg sunitinib (nephrectomy). Untreated biopsy and nephrectomy samples before and after renal artery ligation were controls. Ion Proton sequencing of 48 key ccRCC genes, and MethylCap-seq DNA methylation analysis was performed, data was analysed using the statistical computing environment R. Results Unsupervised hierarchical clustering revealed complete methylome clustering of biopsy and three nephrectomy samples for each patient (14/14 patients). For mutational status, untreated biopsy and all treated nephrectomy samples clustered together in 8/13 (61.5%) patients. The only methylation target significantly altered following sunitinib therapy was VHL promoter region 7896829 which was hypermethylated with treatment (FDR=0.077, P<0.001) and consistent for all patients (pre-treatment 50% patients had VHL mutations, 14% patients VHL hypermethylation). Renal artery ligation did not affect this result. No significant differences in driver or private mutation count was found with sunitinib treatment. Conclusions Demonstration of relative methylome homogeneity and consistent VHL hypermethylation, after sunitinib, may overcome the hurdle of ITH present at other molecular levels for biomarker research.

Retrospective study of the impact of miniSTRs on forensic DNA profiling of touch DNA samples

The theoretical advantages of miniSTRs are undeniable. Several studies show that miniSTRs are more sensitive and robust in the analysis of low template and degraded DNA. In this study we want to show the overall benefit of using miniSTRs in real forensic casework samples and show the percentage of samples that benefit from analysis with additional miniSTR loci in terms of resulting in a useful profile. The considered samples were 3064 touch DNA samples, analyzed in our accredited routine forensic DNA profiling laboratory between mid 2009 and mid 2013. Of these 3064 samples, 618 samples were analyzed using 13 loci, 532 samples using 15 loci and 1914 samples using 20 loci of which 5 were the mini- and midi-STR loci that were added to the extended European Standard Set (ESS). The retrospective results show a small increased success rate after implementation of extra loci and an even smaller increase after the implementation of the mini- and midi-STR analysis. The percentage of touch DNA samples that benefit from the analysis of additional mini- and midi-STR loci is limited.