Filip Van Nieuwerburgh

Library construction for next-generation sequencing: Overviews and challenges

High-throughput sequencing, also known as next-generation sequencing (NGS), has revolutionized genomic research. In recent years, NGS technology has steadily improved, with costs dropping and the number and range of sequencing applications increasing exponentially. Here, we examine the critical role of sequencing library quality and consider important challenges when preparing NGS libraries from DNA and RNA sources. Factors such as the quantity and physical characteristics of the RNA or DNA source material as well as the desired application (i.e., genome sequencing, targeted sequencing, RNA-seq, ChIP-seq, RIP-seq, and methylation) are addressed in the context of preparing high quality sequencing libraries. In addition, the current methods for preparing NGS libraries from single cells are also discussed.

Ancient Dispersal of the Human Fungal Pathogen Cryptococcus gattii from the Amazon Rainforest

Over the past two decades, several fungal outbreaks have occurred, including the high-profile ‘Vancouver Island’ and ‘Pacific Northwest’ outbreaks, caused by Cryptococcus gattii, which has affected hundreds of otherwise healthy humans and animals. Over the same time period, C. gattii was the cause of several additional case clusters at localities outside of the tropical and subtropical climate zones where the species normally occurs. In every case, the causative agent belongs to a previously rare genotype of C. gattii called AFLP6/VGII, but the origin of the outbreak clades remains enigmatic. Here we used phylogenetic and recombination analyses, based on AFLP and multiple MLST datasets, and coalescence gene genealogy to demonstrate that these outbreaks have arisen from a highly-recombining C. gattii population in the native rainforest of Northern Brazil. Thus the modern virulent C. gattii AFLP6/VGII outbreak lineages derived from mating events in South America and then dispersed to temperate regions where they cause serious infections in humans and animals.

Dissection of the phytohormonal regulation of trichome formation and biosynthesis of the antimalarial compound artemisinin in Artemisia annua plants

• Biosynthesis of the sesquiterpene lactone and potent antimalarial drug artemisinin occurs in glandular trichomes of Artemisia annua plants and is subjected to a strict network of developmental and other regulatory cues. • The effects of three hormones, jasmonate, gibberellin and cytokinin, were studied at the structural and molecular levels in two different A. annua chemotypes by microscopic analysis of gland development, and by targeted metabolite and transcript profiling. Furthermore, a genome-wide cDNA-amplified fragment length polymorphism (AFLP)-based transcriptome profiling was carried out of jasmonate-elicited leaves at different developmental stages. • Although cytokinin and gibberellin positively affected at least one aspect of gland formation, these two hormones did not stimulate artemisinin biosynthesis. Only jasmonate simultaneously promoted gland formation and coordinated transcriptional activation of biosynthetic gene expression, which ultimately led to increased sesquiterpenoid accumulation with chemotype-dependent effects on the distinct pathway branches. Transcriptome profiling revealed a trichome-specific fatty acyl- coenzyme A reductase, trichome-specific fatty acyl-CoA reductase 1 (TFAR1), the expression of which correlates with trichome development and sesquiterpenoid biosynthesis. • TFAR1 is potentially involved in cuticular wax formation during glandular trichome expansion in leaves and flowers of A. annua plants. Analysis of phytohormone-modulated transcriptional regulons provides clues to dissect the concerted regulation of metabolism and development of plant trichomes.

Association between the dopamine D2 receptor TaqI A2 allele and low activity COMT allele with obsessive–compulsive disorder in males

BACKGROUND Mounting evidence suggests the involvement of the dopamine system in the pathophysiology of obsessive-compulsive disorder. METHOD The relationship of the dopamine D(2) receptor (DRD2) TaqI A, and catechol-O-methyl-transferase (COMT) NlaIII High/Low activity polymorphism to obsessive-compulsive disorder (OCD) was examined in a sample of 150 patients and 150 controls. RESULTS OCD patients did not show significant differences in genotype distribution and allele frequency for polymorphisms investigated relative to controls. However, when the sample was stratified by gender, there was a trend to a significant predominance of the DRD2 A2A2 genotype (p=0.049), and a higher frequency of the DRD2 A2 allele (p=0.020) and low-activity COMT allele (p=0.035) in male OCD patients compared to male controls. In addition, we observed an association of the DRD2 A2A2 genotype in patients with an early onset of disease (<or=15 years) (p=0.033). CONCLUSIONS Our findings replicate previous reports and provide support for a potential role of the COMT and DRD2 locus in subgroup of male, early onset patients with OCD.

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.

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.

Molecular diagnostics for congenital hearing loss including 15 deafness genes using a next generation sequencing platform.

BackgroundHereditary hearing loss (HL) can originate from mutations in one of many genes involved in the complex process of hearing. Identification of the genetic defects in patients is currently labor intensive and expensive. While screening with Sanger sequencing for GJB2 mutations is common, this is not the case for the other known deafness genes (> 60). Next generation sequencing technology (NGS) has the potential to be much more cost efficient. Published methods mainly use hybridization based target enrichment procedures that are time saving and efficient, but lead to loss in sensitivity. In this study we used a semi-automated PCR amplification and NGS in order to combine high sensitivity, speed and cost efficiency.ResultsIn this proof of concept study, we screened 15 autosomal recessive deafness genes in 5 patients with congenital genetic deafness. 646 specific primer pairs for all exons and most of the UTR of the 15 selected genes were designed using primerXL. Using patient specific identifiers, all amplicons were pooled and analyzed using the Roche 454 NGS technology. Three of these patients are members of families in which a region of interest has previously been characterized by linkage studies. In these, we were able to identify two new mutations in CDH23 and OTOF. For another patient, the etiology of deafness was unclear, and no causal mutation was found. In a fifth patient, included as a positive control, we could confirm a known mutation in TMC1.ConclusionsWe have developed an assay that holds great promise as a tool for screening patients with familial autosomal recessive nonsyndromal hearing loss (ARNSHL). For the first time, an efficient, reliable and cost effective genetic test, based on PCR enrichment, for newborns with undiagnosed deafness is available.

Assessment of Microbial Diversity in Biofilms Recovered from Endotracheal Tubes Using Culture Dependent and Independent Approaches

Ventilator-associated pneumonia (VAP) is a common nosocomial infection in mechanically ventilated patients. Biofilm formation is one of the mechanisms through which the endotracheal tube (ET) facilitates bacterial contamination of the lower airways. In the present study, we analyzed the composition of the ET biofilm flora by means of culture dependent and culture independent (16 S rRNA gene clone libraries and pyrosequencing) approaches. Overall, the microbial diversity was high and members of different phylogenetic lineages were detected (Actinobacteria, beta-Proteobacteria, Candida spp., Clostridia, epsilon-Proteobacteria, Firmicutes, Fusobacteria and gamma-Proteobacteria). Culture dependent analysis, based on the use of selective growth media and conventional microbiological tests, resulted in the identification of typical aerobic nosocomial pathogens which are known to play a role in the development of VAP, e.g. Staphylococcus aureus and Pseudomonas aeruginosa. Other opportunistic pathogens were also identified, including Staphylococcus epidermidis and Kocuria varians. In general, there was little correlation between the results obtained by sequencing 16 S rRNA gene clone libraries and by cultivation. Pyrosequencing of PCR amplified 16 S rRNA genes of four selected samples resulted in the identification of a much wider variety of bacteria. The results from the pyrosequencing analysis suggest that these four samples were dominated by members of the normal oral flora such as Prevotella spp., Peptostreptococcus spp. and lactic acid bacteria. A combination of methods is recommended to obtain a complete picture of the microbial diversity of the ET biofilm.

The Quorum Sensing Inhibitor Hamamelitannin Increases Antibiotic Susceptibility of Staphylococcus aureus Biofilms by Affecting Peptidoglycan Biosynthesis and eDNA Release.

Treatment of Staphylococcus aureus infections has become increasingly challenging due to the rapid emergence and dissemination of methicillin-resistant strains. In addition, S. aureus reside within biofilms at the site of infection. Few novel antibacterial agents have been developed in recent years and their bacteriostatic or bactericidal activity results in selective pressure, inevitably inducing antimicrobial resistance. Consequently, innovative antimicrobials with other modes of action are urgently needed. One alternative approach is targeting the bacterial quorum sensing (QS) system. Hamamelitannin (2′,5-di-O-galloyl-d-hamamelose; HAM) was previously suggested to block QS through the TraP QS system and was shown to increase S. aureus biofilm susceptibility towards vancomycin (VAN) although mechanistic insights are still lacking. In the present study we provide evidence that HAM specifically affects S. aureus biofilm susceptibility through the TraP receptor by affecting cell wall synthesis and extracellular DNA release of S. aureus. We further provide evidence that HAM can increase the susceptibility of S. aureus biofilms towards different classes of antibiotics in vitro. Finally, we show that HAM increases the susceptibility of S. aureus to antibiotic treatment in in vivo Caenorhabditis elegans and mouse mammary gland infection models.