The Oxford Nanopore Technologies (ONT) MinION is a new sequencing technology that potentially offers read lengths of tens of kilobases (kb) limited only by the length of DNA molecules presented to it. The device has a low capital cost, is by far the most portable DNA sequencer available, and can produce data in real-time. It has numerous prospective applications including improving genome sequence assemblies and resolution of repeat-rich regions. Before such a technology is widely adopted, it is important to assess its performance and limitations in respect of throughput and accuracy. In this study we assessed the performance of the MinION by re-sequencing three bacterial genomes, with very different nucleotide compositions ranging from 28.6% to 70.7%; the high G + C strain was underrepresented in the sequencing reads. We estimate the error rate of the MinION (after base calling) to be 38.2%. Mean and median read lengths were 2 kb and 1 kb respectively, while the longest single read was 98 kb. The whole length of a 5 kb rRNA operon was covered by a single read. As the first nanopore-based single molecule sequencer available to researchers, the MinION is an exciting prospect; however, the current error rate limits its ability to compete with existing sequencing technologies, though we do show that MinION sequence reads can enhance contiguity of de novo assembly when used in conjunction with Illumina MiSeq data.

Horizontal gene transfer (HGT) can radically alter the genomes of microorganisms, providing the capacity to adapt to new lifestyles, environments, and hosts. However, the extent of HGT between eukaryotes is unclear. Using whole-genome, gene-by-gene phylogenetic analysis we demonstrate an extensive pattern of cross-kingdom HGT between fungi and oomycetes. Comparative genomics, including the de novo genome sequence of Hyphochytrium catenoides , a free-living sister of the oomycetes, shows that these transfers largely converge within the radiation of oomycetes that colonize plant tissues. The repertoire of HGTs includes a large number of putatively secreted proteins; for example, 7.6% of the secreted proteome of the sudden oak death parasite Phytophthora ramorum has been acquired from fungi by HGT. Transfers include gene products with the capacity to break down plant cell walls and acquire sugars, nucleic acids, nitrogen, and phosphate sources from the environment. Predicted HGTs also include proteins implicated in resisting plant defense mechanisms and effector proteins for attacking plant cells. These data are consistent with the hypothesis that some oomycetes became successful plant parasites by multiple acquisitions of genes from fungi.

The rice blast fungus Magnaporthe oryzae is one of the most significant pathogens affecting global food security. To cause rice blast disease the fungus elaborates a specialised infection structure called an appressorium. Here, we report genome wide transcriptional profile analysis of appressorium development using next generation sequencing (NGS). We performed both RNA-Seq and High-Throughput SuperSAGE analysis to compare the utility of these procedures for identifying differential gene expression in M. oryzae. We then analysed global patterns of gene expression during appressorium development. We show evidence for large-scale gene expression changes, highlighting the role of autophagy, lipid metabolism and melanin biosynthesis in appressorium differentiation. We reveal the role of the Pmk1 MAP kinase as a key global regulator of appressorium-associated gene expression. We also provide evidence for differential expression of transporter-encoding gene families and specific high level expression of genes involved in quinate uptake and utilization, consistent with pathogen-mediated perturbation of host metabolism during plant infection. When considered together, these data provide a comprehensive high-resolution analysis of gene expression changes associated with cellular differentiation that will provide a key resource for understanding the biology of rice blast disease.

Charcot-Marie-Tooth disease is characterized by length-dependent axonal degeneration with distal sensory loss and weakness, deep-tendon-reflex abnormalities, and skeletal deformities. It is caused by mutations in more than 40 genes. We investigated a four-generation family with 23 members affected by the axonal form (type 2), for which the common causes had been excluded by Sanger sequencing. Exome sequencing of three affected individuals separated by eight meioses identified a single shared novel heterozygous variant, c.917A>G, in DYNC1H1, which encodes the cytoplasmic dynein heavy chain 1 (here, novel refers to a variant that has not been seen in dbSNP131or the August 2010 release of the 1000 Genomes project). Testing of six additional affected family members showed cosegregation and a maximum LOD score of 3.6. The shared DYNC1H1 gene variant is a missense substitution, p.His306Arg, at a highly conserved residue within the homodimerization domain. Three mouse models with different mutations within this domain have previously been reported with age-related progressive loss of muscle bulk and locomotor ability. Cytoplasmic dynein is a large multisubunit motor protein complex and has a key role in retrograde axonal transport in neurons. Our results highlight the importance of dynein and retrograde axonal transport in neuronal function in humans.

Robust seasonal dynamics in microbial community composition have previously been observed in the English Channel L4 marine observatory. These could be explained either by seasonal changes in the taxa present at the L4 site, or by the continuous modulation of abundance of taxa within a persistent microbial community. To test these competing hypotheses, deep sequencing of 16S rRNA from one randomly selected time point to a depth of 10 729 927 reads was compared with an existing taxonomic survey data covering 6 years. When compared against the 6-year survey of 72 shallow sequenced time points, the deep sequenced time point maintained 95.4% of the combined shallow OTUs. Additionally, on average, 99.75%±0.06 (mean±s.d.) of the operational taxonomic units found in each shallow sequenced sample were also found in the single deep sequenced sample. This suggests that the vast majority of taxa identified in this ecosystem are always present, but just in different proportions that are predictable. Thus observed changes in community composition are actually variations in the relative abundance of taxa, not, as was previously believed, demonstrating extinction and recolonization of taxa in the ecosystem through time.

The Cystic Fibrosis (CF) lung harbors a complex, polymicrobial ecosystem, in which Pseudomonas aeruginosa is capable of sustaining chronic infections, which are highly resistant to multiple antibiotics. Here, we investigate the phenotypic and genotypic diversity of 44 morphologically identical P. aeruginosa isolates taken from a single CF patient sputum sample. Comprehensive phenotypic analysis of isolates revealed large variances and trade-offs in growth, virulence factors and quorum sensing (QS) signals. Whole genome analysis of 22 isolates revealed high levels of intra-isolate diversity ranging from 5 to 64 SNPs and that recombination and not spontaneous mutation was the dominant driver of diversity in this population. Furthermore, phenotypic differences between isolates were not linked to mutations in known genes but were statistically associated with distinct recombination events. We also assessed antibiotic susceptibility of all isolates. Resistance to antibiotics significantly increased when multiple isolates were mixed together. Our results highlight the significant role of recombination in generating phenotypic and genetic diversification during in vivo chronic CF infection. We also discuss (i) how these findings could influence how patient-to-patient transmission studies are performed using whole genome sequencing, and (ii) the need to refine antibiotic susceptibility testing in sputum samples taken from patients with CF.

During selenate respiration by Thauera selenatis, the reduction of selenate results in the formation of intracellular selenium (Se) deposits that are ultimately secreted as Se nanospheres of approximately 150 nm in diameter. We report that the Se nanospheres are associated with a protein of approximately 95 kDa. Subsequent experiments to investigate the expression and secretion profile of this protein have demonstrated that it is up-regulated and secreted in response to increasing selenite concentrations. The protein was purified from Se nanospheres, and peptide fragments from a tryptic digest were used to identify the gene in the draft T. selenatis genome. A matched open reading frame was located, encoding a protein with a calculated mass of 94.5 kDa. N-terminal sequence analysis of the mature protein revealed no cleavable signal peptide, suggesting that the protein is exported directly from the cytoplasm. The protein has been called Se factor A (SefA), and homologues of known function have not been reported previously. The sefA gene was cloned and expressed in Escherichia coli, and the recombinant His-tagged SefA purified. In vivo experiments demonstrate that SefA forms larger (approximately 300 nm) Se nanospheres in E. coli when treated with selenite, and these are retained within the cell. In vitro assays demonstrate that the formation of Se nanospheres upon the reduction of selenite by glutathione are stabilized by the presence of SefA. The role of SefA in selenium nanosphere assembly has potential for exploitation in bionanomaterial fabrication.nanoparticles | biomineralization | anaerobic respiration

Microarray and next-generation sequencing techniques which allow whole genome analysis of chromatin structure and sequence-specific protein binding are revolutionizing our view of chromosome architecture and function. However, many current methods in this field rely on biochemical purification of highly specific fractions of DNA prepared from chromatin digested with either micrococcal nuclease or DNaseI and are restricted in the parameters they can measure. Here, we show that a broad size-range of genomic DNA species, produced by partial micrococcal nuclease digestion of chromatin, can be sequenced using paired-end mode next-generation technology. The paired sequence reads, rather than DNA molecules, can then be size-selected and mapped as particle classes to the target genome. Using budding yeast as a model, we show that this approach reveals position and structural information for a spectrum of nuclease resistant complexes ranging from transcription factor-bound DNA elements up to mono- and poly-nucleosomes. We illustrate the utility of this approach in visualizing the MNase digestion landscape of protein-coding gene transcriptional start sites, and demonstrate a comparative analysis which probes the function of the chromatin-remodelling transcription factor Cbf1p.