Toon Rosseel

DNase SISPA-Next Generation Sequencing Confirms Schmallenberg Virus in Belgian Field Samples and Identifies Genetic Variation in Europe

In 2011, a novel Orthobunyavirus was identified in cattle and sheep in Germany and the Netherlands. This virus was named Schmallenberg virus (SBV). Later, presence of the virus was confirmed using real time RT-PCR in cases of congenital malformations of bovines and ovines in several European countries, including Belgium. In the absence of specific sequencing protocols for this novel virus we confirmed its presence in RT-qPCR positive field samples using DNase SISPA-next generation sequencing (NGS), a virus discovery method based on random amplification and next generation sequencing. An in vitro transcribed RNA was used to construct a standard curve allowing the quantification of viral RNA in the field samples. Two field samples of aborted lambs containing 7.66 and 7.64 log10 RNA copies per µL total RNA allowed unambiguous identification of SBV. One sample yielded 192 SBV reads covering about 81% of the L segment, 56% of the M segment and 13% of the S segment. The other sample resulted in 8 reads distributed over the L and M segments. Three weak positive field samples (one from an aborted calf, two from aborted lambs) containing virus quantities equivalent to 4.27–4.89 log10 RNA copies per µL did not allow identification using DNase SISPA-NGS. This partial sequence information was compared to the whole genome sequence of SBV isolated from bovines in Germany, identifying several sequence differences. The applied viral discovery method allowed the confirmation of SBV in RT-qPCR positive brain samples. However, the failure to confirm SBV in weak PCR-positive samples illustrates the importance of the selection of properly targeted and fresh field samples in any virus discovery method. The partial sequences derived from the field samples showed several differences compared to the sequences from bovines in Germany, indicating sequence divergence within the epidemic.

The Origin of Biased Sequence Depth in Sequence-Independent Nucleic Acid Amplification and Optimization for Efficient Massive Parallel Sequencing

Sequence Independent Single Primer Amplification is one of the most widely used random amplification approaches in virology for sequencing template preparation. This technique relies on oligonucleotides consisting of a 3′ random part used to prime complementary DNA synthesis and a 5′ defined tag sequence for subsequent amplification. Recently, this amplification method was combined with next generation sequencing to obtain viral sequences. However, these studies showed a biased distribution of the resulting sequence reads over the analyzed genomes. The aim of this study was to elucidate the mechanisms that lead to biased sequence depth when using random amplification. Avian paramyxovirus type 8 was used as a model RNA virus to investigate these mechanisms. We showed, based on in silico analysis of the sequence depth in relation to GC-content, predicted RNA secondary structure and sequence complementarity to the 3′ part of the tag sequence, that the tag sequence has the main contribution to the observed bias in sequence depth. We confirmed this finding experimentally using both fragmented and non-fragmented viral RNAs as well as primers differing in random oligomer length (6 or 12 nucleotides) and in the sequence of the amplification tag. The observed oligonucleotide annealing bias can be reduced by extending the random oligomer sequence and by in silico combining sequence data from SISPA experiments using different 5′ defined tag sequences. These findings contribute to the optimization of random nucleic acid amplification protocols that are currently required for downstream applications such as viral metagenomics and microarray analysis.

Identification and complete genome sequencing of paramyxoviruses in mallard ducks (Anas platyrhynchos) using random access amplification and next generation sequencing technologies

BackgroundDuring a wildlife screening program for avian influenza A viruses (AIV) and avian paramyxoviruses (APMV) in Belgium, we isolated two hemagglutinating agents from pools of cloacal swabs of wild mallards (Anas platyrhynchos) caught in a single sampling site at two different times. AIV and APMV1 were excluded using hemagglutination inhibition (HI) testing and specific real-time RT-PCR tests.MethodsTo refine the virological identification of APMV2-10 realized by HI subtyping tests and in lack of validated molecular tests for APMV2-10, random access amplification was used in combination with next generation sequencing for the sequence independent identification of the viruses and the determination of their genomes.ResultsThree different APMVs were identified. From one pooled sample, the complete genome sequence (15054 nucleotides) of an APMV4 was assembled from the random sequences. From the second pooled sample, the nearly complete genome sequence of an APMV6 (genome size of 16236 nucleotides) was determined, as well as a partial sequence for an APMV4. This APMV4 was closely related but not identical to the APMV4 isolated from the first sample. Although a cross-reactivity with other APMV subtypes did not allow formal identification, the HI subtyping revealed APMV4 and APMV6 in the respective pooled samples but failed to identify the co-infecting APMV4 in the APMV6 infected pool.ConclusionsThese data further contribute to the knowledge about the genetic diversity within the serotypes APMV4 and 6, and confirm the limited sensitivity of the HI subtyping test. Moreover, this study demonstrates the value of a random access nucleic acid amplification method in combination with massive parallel sequencing. Using only a moderate and economical sequencing effort, the characterization and full genome sequencing of APMVs can be obtained, including the identification of viruses in mixed infections.

NR5A1 is a novel disease gene for 46,XX testicular and ovotesticular disorders of sex development

Purpose:We aimed to identify the genetic cause in a cohort of 11 unrelated cases and two sisters with 46,XX SRY-negative (ovo)testicular disorders of sex development (DSD).Methods:Whole-exome sequencing (n = 9), targeted resequencing (n = 4), and haplotyping were performed. Immunohistochemistry of sex-specific markers was performed on patients’ gonads. The consequences of mutation were investigated using luciferase assays, localization studies, and RNA-seq.Results:We identified a novel heterozygous NR5A1 mutation, c.274C>T p.(Arg92Trp), in three unrelated patients. The Arg92 residue is highly conserved and located in the Ftz-F1 region, probably involved in DNA-binding specificity and stability. There were no consistent changes in transcriptional activation or subcellular localization. Transcriptomics in patient-derived lymphocytes showed upregulation of MAMLD1, a direct NR5A1 target previously associated with 46,XY DSD. In gonads of affected individuals, ovarian FOXL2 and testicular SRY-independent SOX9 expression observed.Conclusions:We propose NR5A1, previously associated with 46,XY DSD and 46,XX primary ovarian insufficiency, as a novel gene for 46,XX (ovo)testicular DSD. We hypothesize that p.(Arg92Trp) results in decreased inhibition of the male developmental pathway through downregulation of female antitestis genes, thereby tipping the balance toward testicular differentiation in 46,XX individuals. In conclusion, our study supports a role for NR5A1 in testis differentiation in the XX gonad.Genet Med 19 4, 367–376.

Evaluation of convenient pretreatment protocols for RNA virus metagenomics in serum and tissue samples

Viral metagenomic approaches are increasingly being used for viral discovery. Various strategies are applied to enrich viral sequences, but there is often a lack of knowledge about their effective influence on the viral discovery sensitivity. We evaluate some convenient and widely used approaches for RNA virus discovery in clinical samples in order to reveal their sensitivity and potential bias introduced by the enrichment or amplifications steps. An RNA virus was artificially spiked at a fixed titer in serum and lung tissue, respectively, low and high nucleic acid content matrices. For serum, a simple DNase treatment on the RNA extract gave the maximum gain in proportion of viral sequences (83×), and a subsequent ribosomal RNA removal nearly doubled once more the proportion of viral sequences. For lung tissue, a ribosomal RNA depletion step on the RNA extract had the biggest gain in proportion of viral sequences (32×). We show also that direct sequencing of cDNA is recommended above an extra random PCR amplification step, and a that the virion enrichment strategy (filtration and nuclease treatment) has a beneficial effect for sequencing-based virus discovery. Our findings provide sample-dependent guidelines for targeted virus discovery strategies.

Phylogeographic analysis of avian influenza viruses isolated from Charadriiformes in Belgium confirms intercontinental reassortment in gulls

Nine influenza viruses isolated from gulls and shorebirds in Belgium (2008-2010), including H3N8, H5N2, H6N1, H11N9, H13N6, H13N8, and H16N3 subtypes, were targeted using random amplification and next-generation sequencing. The gene segments of these viruses segregated into three phylogeographic lineage types: (1) segments circulating in waterfowl in Eurasia with sporadic introduction in other species and in the Americas (“Eurasian avian”), (2) segments circulating in American waterfowl with sporadic introduction to other species and regions (“American avian”), and (3) segments circulating exclusively in gulls and shorebirds and having increased connectivity between the two hemispheres (“Charadriiformes specific”). Notably, an H6N1 and an H5N2 isolated from L. argentatus had mainly Eurasian avian genes but shared a matrix segment of American avian origin (first documentation in European gulls of transhemispheric reassortment). These data support the growing evidence of an important role of Charadriiformes birds in the dynamic nature of avian influenza ecology.

Next-generation sequencing in veterinary medicine : how can the massive amount of information arising from high-throughput technologies improve diagnosis, control, and management of infectious diseases?

The development of high-throughput molecular technologies and associated bioinformatics has dramatically changed the capacities of scientists to produce, handle, and analyze large amounts of genomic, transcriptomic, and proteomic data. A clear example of this step-change is represented by the amount of DNA sequence data that can be now produced using next-generation sequencing (NGS) platforms. Similarly, recent improvements in protein and peptide separation efficiencies and highly accurate mass spectrometry have promoted the identification and quantification of proteins in a given sample. These advancements in biotechnology have increasingly been applied to the study of animal infectious diseases and are beginning to revolutionize the way that biological and evolutionary processes can be studied at the molecular level. Studies have demonstrated the value of NGS technologies for molecular characterization, ranging from metagenomic characterization of unknown pathogens or microbial communities to molecular epidemiology and evolution of viral quasispecies. Moreover, high-throughput technologies now allow detailed studies of host-pathogen interactions at the level of their genomes (genomics), transcriptomes (transcriptomics), or proteomes (proteomics). Ultimately, the interaction between pathogen and host biological networks can be questioned by analytically integrating these levels (integrative OMICS and systems biology). The application of high-throughput biotechnology platforms in these fields and their typical low-cost per information content has revolutionized the resolution with which these processes can now be studied. The aim of this chapter is to provide a current and prospective view on the opportunities and challenges associated with the application of massive parallel sequencing technologies to veterinary medicine, with particular focus on applications that have a potential impact on disease control and management.

What's in a strain? Viral metagenomics identifies genetic variation and contaminating circoviruses in laboratory isolates of pigeon paramyxovirus type 1

We used next generation sequencing on random amplified viral nucleic acids to determine the genome sequence of 11 pigeon paramyxovirus type 1 (PPMV-1) isolates from Belgium (period 1998-2011). The PPMV-1 deep sequence data allowed identification of sequence variability in multiple PPMV-1 isolates, including one STOP codon in the Matrix gene which was present in 15% of the viral population of one isolate. Notably, mutations that were previously associated with pathogenicity in chickens were identified as minor sequence variants in one parent laboratory strain. A phylogenetic analysis of the consensus PPMV-1 genome sequences was performed. In addition to providing nearly complete paramyxovirus genome sequences, our sequence-independent approach identified the presence of pigeon circovirus (PiCV) sequences in four of these viral stocks. Real-time quantitative RT-PCR analysis specific for PMV-1 and PiCV showed that these contaminations were present in seven viral stocks consisting of allantoic fluids and was occasionally also detected in stocks passaged in embryonated chicken eggs. Phylogenetic analysis of the PiCV consensus genome sequences showed a circulation of PiCV covering the full genetic diversity of known PiCV. This study shows the value of novel sequence independent technologies for access to sequence information for the control of reference virus stocks and other biological materials, as co-infecting viruses or sequence variants from the original sample may persist in the stocks without being identified by the routine virus-specific diagnostic tools. The exact role of PiCV in pigeon disease - in particular Newcastle disease - and its potential interference with PPMV-1 diagnostics remains to be investigated.

Extensive Genetic Variability Linked to IS26 Insertions in the fljB Promoter Region of Atypical Monophasic Variants of Salmonella enterica Serovar Typhimurium

ABSTRACT Fifty-nine monophasic Salmonella enterica serovar Typhimurium isolates, collected in Belgium during the period from 2008 to 2011, have been serotyped as 4,[5]:i:− and shown to harbor an fljB coding sequence. The genetic differences between these strains and phenotypically biphasic Salmonella Typhimurium were analyzed through PCR and DNA sequencing. Genetic alterations in the fljB promoter region affecting expression of the phase 2 flagellin were observed in 53 isolates. Other genetic events in the invertible region carrying the fljB promoter were observed in 2 isolates. For the remaining 4 isolates, no molecular differences with a reference biphasic Salmonella Typhimurium strain could be observed. Next-generation sequencing of one representative isolate affected in the fljB promoter region revealed a 26-kb IS26 composite transposon insertion along with a local genomic rearrangement. Several other IS26 element-mediated alterations of this genomic region were observed. This group of monophasic Salmonella Typhimurium isolates was genetically heterogeneous, as revealed by multilocus variable-number tandem-repeat analysis (MLVA), PCR, and sequencing. Pigs and pork represented a major source of such monophasic isolates in Belgium, as reported in other countries. Three out of 5 isolates of human origin presented genetic profiles identical to those of food isolates, demonstrating the pathogenic potential of the newly characterized variants and potential dissemination along the food chain. This study highlighted the key role played by IS26 insertions in the loss of phase 2 flagellin expression and the subsequent generation of multiple monophasic variant lineages from biphasic Salmonella Typhimurium ancestors.

Complete Genome Sequence of Bovine Polyomavirus Type 1 from Aborted Cattle, Isolated in Belgium in 2014

ABSTRACT The complete and fully annotated genome sequence of a bovine polyomavirus type 1 (BPyV/BEL/1/2014) from aborted cattle was assembled from a metagenomics data set. The 4,697-bp circular dsDNA genome contains 6 protein-coding genes. Bovine polyomavirus is unlikely to be causally related to the abortion cases.