DNA secondary structures are associated with recombination in major Plasmodium falciparum variable surface antigen gene families.

Many bacterial, viral and parasitic pathogens undergo antigenic variation to counter host immune defense mechanisms. In Plasmodium falciparum, the most lethal of human malaria parasites, switching of var gene expression results in alternating expression of the adhesion proteins of the Plasmodium falciparum-erythrocyte membrane protein 1 class on the infected erythrocyte surface. Recombination clearly generates var diversity, but the nature and control of the genetic exchanges involved remain unclear. By experimental and bioinformatic identification of recombination events and genome-wide recombination hotspots in var genes, we show that during the parasite's sexual stages, ectopic recombination between isogenous var paralogs occurs near low folding free energy DNA 50-mers and that these sequences are heavily concentrated at the boundaries of regions encoding individual Plasmodium falciparum-erythrocyte membrane protein 1 structural domains. The recombinogenic potential of these 50-mers is not parasite-specific because these sequences also induce recombination when transferred to the yeast Saccharomyces cerevisiae. Genetic cross data suggest that DNA secondary structures (DSS) act as inducers of recombination during DNA replication in P. falciparum sexual stages, and that these DSS-regulated genetic exchanges generate functional and diverse P. falciparum adhesion antigens. DSS-induced recombination may represent a common mechanism for optimizing the evolvability of virulence gene families in pathogens.

One hundred twenty years of koala retrovirus evolution determined from museum skins.

Although endogenous retroviruses are common across vertebrate genomes, the koala retrovirus (KoRV) is the only retrovirus known to be currently invading the germ line of its host. KoRV is believed to have first infected koalas in northern Australia less than two centuries ago. We examined KoRV in 28 koala museum skins collected in the late 19th and 20th centuries and deep sequenced the complete proviral envelope region from five northern Australian specimens. Strikingly, KoRV env sequences were conserved among koalas collected over the span of a century, and two functional motifs that affect viral infectivity were fixed across the museum koala specimens. We detected only 20 env polymorphisms among the koalas, likely representing derived mutations subject to purifying selection. Among northern Australian koalas, KoRV was already ubiquitous by the late 19th century, suggesting that KoRV evolved and spread among koala populations more slowly than previously believed. Given that museum and modern koalas share nearly identical KoRV sequences, it is likely that koala populations, for more than a century, have experienced increased susceptibility to diseases caused by viral pathogenesis.

Genetic diversity among pandemic 2009 influenza viruses isolated from a transmission chain.

BACKGROUND: Influenza viruses such as swine-origin influenza A(H1N1) virus (A(H1N1)pdm09) generate genetic diversity due to the high error rate of their RNA polymerase, often resulting in mixed genotype populations (intra-host variants) within a single infection. This variation helps influenza to rapidly respond to selection pressures, such as those imposed by the immunological host response and antiviral therapy. We have applied deep sequencing to characterize influenza intra-host variation in a transmission chain consisting of three cases due to oseltamivir-sensitive viruses, and one derived oseltamivir-resistant case. METHODS: Following detection of the A(H1N1)pdm09 infections, we deep-sequenced the complete NA gene from two of the oseltamivir-sensitive virus-infected cases, and all eight gene segments of the viruses causing the remaining two cases. RESULTS: No evidence for the resistance-causing mutation (resulting in NA H275Y substitution) was observed in the oseltamivir-sensitive cases. Furthermore, deep sequencing revealed a subpopulation of oseltamivir-sensitive viruses in the case carrying resistant viruses. We detected higher levels of intra-host variation in the case carrying oseltamivir-resistant viruses than in those infected with oseltamivir-sensitive viruses. CONCLUSIONS: Oseltamivir-resistance was only detected after prophylaxis with oseltamivir, suggesting that the mutation was selected for as a result of antiviral intervention. The persisting oseltamivir-sensitive virus population in the case carrying resistant viruses suggests either that a small proportion survive the treatment, or that the oseltamivir-sensitive virus rapidly re-establishes itself in the virus population after the bottleneck. Moreover, the increased intra-host variation in the oseltamivir-resistant case is consistent with the hypothesis that the population diversity of a RNA virus can increase rapidly following a population bottleneck.

Second-generation sequencing of forensic STRs using the Ion Torrent™ HID STR 10-plex and the Ion PGM™.

Second-generation sequencing (SGS) using Roche/454 and Illumina platforms has proved capable of sequencing the majority of the key forensic genetic STR systems. Given that Roche has announced that the 454 platforms will no longer be supported from 2015, focus should now be shifted to competing SGS platforms, such as the MiSeq (Illumina) and the Ion Personal Genome Machine (Ion PGM™; Thermo Fisher). There are currently several challenges faced with amplicon-based SGS STR typing in forensic genetics, including current lengths of amplicons for CE-typing and lack of uniform data analysis between laboratories. Thermo Fisher has designed a human identification (HID) short tandem repeat (STR) 10-plex panel including amelogenin, CSF1PO, D16S539, D3S1358, D5S818, D7S820, D8S1179, TH01, TPOX and vWA, where the primers have been designed specifically for the purpose of SGS and the data analysis is supported by Ion Torrent™ software. Hence, the combination of the STR 10-plex and the Ion PGM™ represents the first fully integrated SGS STR typing solution from PCR to data analysis. In this study, four experiments were performed to evaluate the alpha-version of the STR 10-plex: (1) typing of control samples; (2) analysis of sensitivity; (3) typing of mixtures; and (4) typing of biological crime case samples. Full profiles and concordant results between replicate SGS runs and CE-typing were observed for all control samples. Full profiles were seen with DNA input down to 50 pg, with the exception of a single locus drop-out in one of the 100 pg dilutions. Mixtures were easily deconvoluted down to 20:1, although alleles from the minor contributor had to be identified manually as some signals were not called by the Ion Torrent™ software. Interestingly, full profiles were obtained for all biological samples from real crime and identification cases, in which only partial profiles were obtained with PCR-CE assays. In conclusion, the Ion Torrent™ HID STR 10-plex panel offers an all-in-one solution from amplification of STRs and amelogenin, and sequencing to data analysis.

Evaluation of the Ion Torrent™ HID SNP 169-plex: A SNP typing assay developed for human identification by second generation sequencing.

The Ion Torrent™ HID SNP assay amplified 136 autosomal SNPs and 33 Y-chromosome markers in one PCR and the markers were subsequently typed using the Ion PGM™ second generation sequencing platform. A total of 51 of the autosomal SNPs were selected from the SNPforID panel that is routinely used in our ISO 17025 accredited laboratory. Concordance between the Ion Torrent™ HID SNP assay and the SNPforID assay was tested by typing 44 Iraqis twice with the Ion Torrent™ HID SNP assay. The same samples were previously typed with the SNPforID assay and the Y-chromosome haplogroups of the individuals were previously identified by typing 45 Y-chromosome SNPs. Full concordance between the assays were obtained except for the SNP genotypes of two SNPs. These SNPs were among the eight SNPs (rs2399332, rs1029047, rs10776839, rs4530059, rs8037429, rs430046, rs1031825 and rs1523537) with inconsistent allele balance among samples. These SNPs should be excluded from the panel. The optimal amount of DNA in the PCR seemed to be ≥0.5ng. Allele drop-outs were rare and only seen in experiments with <0.5ng input DNA and with a coverage of <50reads. No allele drop-in was observed. The great majority of the heterozygote allele balances were between 0.6 and 1.6, which is comparable to the heterozygote balances of STRs typed with PCR-CE. The number of reads with base calls that differed from the genotype call was typically less than five. This allowed detection of 1:100 mixtures with a high degree of certainty in experiments with a high total depth of coverage. In conclusion, the Ion PGM™ is a very promising platform for forensic genetics. However, the secondary sequence analysis software made wrong genotype calls from correctly sequenced alleles. These types of errors must be corrected before the platform can be used in case work. Furthermore, the sequence analysis software should be further developed and include quality settings for each SNP based on validation studies.

Long-term RNA persistence in postmortem contexts.

Ribonucleic acids (RNA) are generally considered fragile molecules that are readily degraded. However, there is growing documentation of long-term (from days to centuries) RNA persistence in a variety of contexts and tissue types, and as such a number of academic disciplines are beginning to exploit degraded RNA. While the reasons for its survival are not fully understood, there are several plausible mechanisms that would safeguard this molecule against degradation. However, after examining the literature available on the postmortem instability and decay mechanisms of RNA, it has become clear that limited experimental studies and no reviews offer an overview of these mechanisms. Hence in this review we outline molecular reasons for RNA surviving long-term postmortem, and provide specific examples of RNA survival in forensic, archival and archaeological contexts. A better understanding of the mechanisms of RNA decay will be crucial for developing expectations on its long-term survival.

Deep sequencing of RNA from ancient maize kernels.

The characterization of biomolecules from ancient samples can shed otherwise unobtainable insights into the past. Despite the fundamental role of transcriptomal change in evolution, the potential of ancient RNA remains unexploited - perhaps due to dogma associated with the fragility of RNA. We hypothesize that seeds offer a plausible refuge for long-term RNA survival, due to the fundamental role of RNA during seed germination. Using RNA-Seq on cDNA synthesized from nucleic acid extracts, we validate this hypothesis through demonstration of partial transcriptomal recovery from two sources of ancient maize kernels. The results suggest that ancient seed transcriptomics may offer a powerful new tool with which to study plant domestication.

Roche genome sequencer FLX based high-throughput sequencing of ancient DNA.

Since the development of so-called "next generation" high-throughput sequencing in 2005, this technology has been applied to a variety of fields. Such applications include disease studies, evolutionary investigations, and ancient DNA. Each application requires a specialized protocol to ensure that the data produced is optimal. Although much of the procedure can be followed directly from the manufacturer's protocols, the key differences lie in the library preparation steps. This chapter presents an optimized protocol for the sequencing of fossil remains and museum specimens, commonly referred to as "ancient DNA," using the Roche GS FLX 454 platform.

A simple method for the parallel deep sequencing of full influenza A genomes.

Given the major threat of influenza A to human and animal health, and its ability to evolve rapidly through mutation and reassortment, tools that enable its timely characterization are necessary to help monitor its evolution and spread. For this purpose, deep sequencing can be a very valuable tool. This study reports a comprehensive method that enables deep sequencing of the complete genomes of influenza A subtypes using the Illumina Genome Analyzer IIx (GAIIx). By using this method, the complete genomes of nine viruses were sequenced in parallel, representing the 2009 pandemic H1N1 virus, H5N1 virus from human and H1N1 virus from swine, on a single lane of a GAIIx flow cell to an average depth of 122-fold. This technique can be applied to cultivated and uncultivated virus.

Application and comparison of large-scale solution-based DNA capture-enrichment methods on ancient DNA.

The development of second-generation sequencing technologies has greatly benefitted the field of ancient DNA (aDNA). Its application can be further exploited by the use of targeted capture-enrichment methods to overcome restrictions posed by low endogenous and contaminating DNA in ancient samples. We tested the performance of Agilent's SureSelect and Mycroarray's MySelect in-solution capture systems on Illumina sequencing libraries built from ancient maize to identify key factors influencing aDNA capture experiments. High levels of clonality as well as the presence of multiple-copy sequences in the capture targets led to biases in the data regardless of the capture method. Neither method consistently outperformed the other in terms of average target enrichment, and no obvious difference was observed either when two tiling designs were compared. In addition to demonstrating the plausibility of capturing aDNA from ancient plant material, our results also enable us to provide useful recommendations for those planning targeted-sequencing on aDNA.

High-throughput sequencing of core STR loci for forensic genetic investigations using the Roche Genome Sequencer FLX platform.

The analysis and profiling of short tandem repeat (STR) loci is routinely used in forensic genetics. Current methods to investigate STR loci, including PCR-based standard fragment analyses and capillary electrophoresis, only provide amplicon lengths that are used to estimate the number of STR repeat units. These methods do not allow for the full resolution of STR base composition that sequencing approaches could provide. Here we present an STR profiling method based on the use of the Roche Genome Sequencer (GS) FLX to simultaneously sequence multiple core STR loci. Using this method in combination with a bioinformatic tool designed specifically to analyze sequence lengths and frequencies, we found that GS FLX STR sequence data are comparable to conventional capillary electrophoresis-based STR typing. Furthermore, we found DNA base substitutions and repeat sequence variations that would not have been identified using conventional STR typing.