Peter A. Stockwell

Pseudomonas syringae pv. actinidiae from Recent Outbreaks of Kiwifruit Bacterial Canker Belong to Different Clones That Originated in China

A recently emerged plant disease, bacterial canker of kiwifruit (Actinidia deliciosa and A. chinensis), is caused by Pseudomonas syringae pv. actinidiae (PSA). The disease was first reported in China and Japan in the 1980s. A severe outbreak of PSA began in Italy in 2008 and has spread to other European countries. PSA was found in both New Zealand and Chile in 2010. To study the evolution of the pathogen and analyse the transmission of PSA between countries, genomes of strains from China and Japan (where the genus Actinidia is endemic), Italy, New Zealand and Chile were sequenced. The genomes of PSA strains are very similar. However, all strains from New Zealand share several single nucleotide polymorphisms (SNPs) that distinguish them from all other PSA strains. Similarly, all the PSA strains from the 2008 Italian outbreak form a distinct clonal group and those from Chile form a third group. In addition to the rare SNPs present in the core genomes, there is abundant genetic diversity in a genomic island that is part of the accessory genome. The island from several Chinese strains is almost identical to the island present in the New Zealand strains. The island from a different Chinese strain is identical to the island present in the strains from the recent Italian outbreak. The Chilean strains of PSA carry a third variant of this island. These genomic islands are integrative conjugative elements (ICEs). Sequencing of these ICEs provides evidence of three recent horizontal transmissions of ICE from other strains of Pseudomonas syringae to PSA. The analyses of the core genome SNPs and the ICEs, combined with disease history, all support the hypothesis of an independent Chinese origin for both the Italian and the New Zealand outbreaks and suggest the Chilean strains also originate from China.

Comparison of alignment software for genome-wide bisulphite sequence data.

Recent advances in next generation sequencing (NGS) technology now provide the opportunity to rapidly interrogate the methylation status of the genome. However, there are challenges in handling and interpretation of the methylation sequence data because of its large volume and the consequences of bisulphite modification. We sequenced reduced representation human genomes on the Illumina platform and efficiently mapped and visualized the data with different pipelines and software packages. We examined three pipelines for aligning bisulphite converted sequencing reads and compared their performance. We also comment on pre-processing and quality control of Illumina data. This comparison highlights differences in methods for NGS data processing and provides guidance to advance sequence-based methylation data analysis for molecular biologists.

Transterm: a database to aid the analysis of regulatory sequences in mRNAs

Messenger RNAs, in addition to coding for proteins, may contain regulatory elements that affect how the protein is translated. These include protein and microRNA-binding sites. Transterm (http://mRNA.otago.ac.nz/Transterm.html) is a database of regions and elements that affect translation with two major unique components. The first is integrated results of analysis of general features that affect translation (initiation, elongation, termination) for species or strains in Genbank, processed through a standard pipeline. The second is curated descriptions of experimentally determined regulatory elements that function as translational control elements in mRNAs. Transterm focuses on protein binding sites, particularly those in 3′-untranslated regions (3′-UTR). For this release the interface has been extensively updated based on user feedback. The data is now accessible by strain rather than species, for example there are 10 Escherichia coli strains (genomes) analysed separately. In addition to providing a repository of data, the database also provides tools for users to query their own mRNA sequences. Users can search sequences for Transterm or user defined regulatory elements, including protein or miRNA targets. Transterm also provides a central core of links to related resources for complementary analyses.

The translational termination signal database

The Translational Termination Database (TransTerm) consists of the immediate context sequences around the natural termination codons from 45 organisms, and summary tables. The influence of termination codon context on their effectivness as stop signals has been widely documented. The SPECIES--TRI.DAT table shows trinucleotide stop codon usage in each organism and for comparison the occurrence of these sequences in the noncoding region. The SPECIES--TETRA.DAT table contains is a similar table of tetranucleotide stop signal usage. The database is available from EMBL.

DMAP: differential methylation analysis package for RRBS and WGBS data

MOTIVATION The rapid development of high-throughput sequencing technologies has enabled epigeneticists to quantify DNA methylation on a massive scale. Progressive increase in sequencing capacity present challenges in terms of processing analysis and the interpretation of the large amount of data; investigating differential methylation between genome-scale data from multiple samples highlights this challenge. RESULTS We have developed a differential methylation analysis package (DMAP) to generate coverage-filtered reference methylomes and to identify differentially methylated regions across multiple samples from reduced representation bisulphite sequencing and whole genome bisulphite sequencing experiments. We introduce a novel fragment-based approach for investigating DNA methylation patterns for reduced representation bisulphite sequencing data. Further, DMAP provides the identity of gene and CpG features and distances to the differentially methylated regions in a format that is easily analyzed with limited bioinformatics knowledge. AVAILABILITY AND IMPLEMENTATION The software has been implemented in C and has been written to ensure portability between different platforms. The source code and documentation is freely available (DMAP: as compressed TAR archive folder) from http://biochem.otago.ac.nz/research/databases-software/. Two test datasets are also available for download from the Web site. Test dataset 1 contains reads from chromosome 1 of a patient and a control, which is used for comparative analysis in the current article. Test dataset 2 contains reads from a part of chromosome 21 of three disease and three control samples for testing the operation of DMAP, especially for the analysis of variance. Example commands for the analyses are included.

Self-programmable, self-assembling two-dimensional genetic matter.

Putative two-dimensional coding systems can beconstructed from aqueous solutions of purine andpyrimidine nucleic acid bases evaporated at moderatetemperatures on the surfaces of inorganic solids. Theresultant structures are monolayers which are formedspontaneously by molecular self-assembly and they havebeen observed with molecular resolution by scanningtunnelling microscopy (STM). When formed fromsolutions of a single base, the monolayers of adenineand uracil have crystalline characteristics and theSTM images can be interpreted in terms of thegeometrical placement of planar arranged moleculesthat interact laterally by intermolecular hydrogenbonding. When formed from solutions containing amixture of adenine and uracil, the monolayers haveaperiodic structures. Small crystalline domainswithin these monolayers can be interpreted in terms ofthe single phase configurations of the molecules andthe remaining aperiodic structures can presumably beinterpreted, geometrically, in terms of the 21theoretically possible adenine-adenine, uracil-uraciland adenine-uracil hydrogen bonding interactions. Wepropose that combinatorial arrangements of planararranged purine and pyrimidine bases could provide thenecessary complexity to act as a primitive geneticmechanism and may have relevance to the origin of life.

Is gene deletion in eukaryotes sequence-dependent A study of nine deletion junctions and nineteen other deletion breakpoints in intron 7 of the human dystrophin gene

Although large deletions comprise 65% of the mutations that underlie most cases of Duchenne and Becker muscular dystrophies, the DNA sequence characteristics of the deletions and the molecular processes leading to their formation are largely unknown. Intron 7 of the human dystrophin gene is unusually large (110 kb) and a substantial number of deletions have been identified with endpoints within this intron. The distribution of 28 deletion endpoints was mapped to local sequence elements by PCR. The break points were distributed among unique sequence, LINE-1, Alu, MIR, MER and microsatellite sequences with frequencies expected from the frequency of those sequences in the intron. Thus, deletions in this intron are not associated primarily with any one of those sequences in the intron. Nine deletion junctions were amplified and sequenced. Eight were deletions between DNA sequences with minimal homology (0-4 bp) and are therefore unlikely to be products of homologous recombination. In the ninth case, a complex rearrangement was found to be consistent with unequal recombinational exchange between two Alu sequences coupled with a duplication. We have hypothesized that a paucity of matrix attachment regions in this very large intron expanded by the insertion of many mobile elements might provoke a chromatin structure that stimulates deletions (McNaughton et al., 1997, Genomics 40, 294-304). The data presented here are consistent with that idea and demonstrate that the deletion sequences are not usually produced by homologous DNA misalignments.

Technical Considerations for Reduced Representation Bisulfite Sequencing with Multiplexed Libraries

Reduced representation bisulfite sequencing (RRBS), which couples bisulfite conversion and next generation sequencing, is an innovative method that specifically enriches genomic regions with a high density of potential methylation sites and enables investigation of DNA methylation at single-nucleotide resolution. Recent advances in the Illumina DNA sample preparation protocol and sequencing technology have vastly improved sequencing throughput capacity. Although the new Illumina technology is now widely used, the unique challenges associated with multiplexed RRBS libraries on this platform have not been previously described. We have made modifications to the RRBS library preparation protocol to sequence multiplexed libraries on a single flow cell lane of the Illumina HiSeq 2000. Furthermore, our analysis incorporates a bioinformatics pipeline specifically designed to process bisulfite-converted sequencing reads and evaluate the output and quality of the sequencing data generated from the multiplexed libraries. We obtained an average of 42 million paired-end reads per sample for each flow-cell lane, with a high unique mapping efficiency to the reference human genome. Here we provide a roadmap of modifications, strategies, and trouble shooting approaches we implemented to optimize sequencing of multiplexed libraries on an a RRBS background.

Mapping the zebrafish brain methylome using reduced representation bisulfite sequencing.

Reduced representation bisulfite sequencing (RRBS) has been used to profile DNA methylation patterns in mammalian genomes such as human, mouse and rat. The methylome of the zebrafish, an important animal model, has not yet been characterized at base-pair resolution using RRBS. Therefore, we evaluated the technique of RRBS in this model organism by generating four single-nucleotide resolution DNA methylomes of adult zebrafish brain. We performed several simulations to show the distribution of fragments and enrichment of CpGs in different in silico reduced representation genomes of zebrafish. Four RRBS brain libraries generated 98 million sequenced reads and had higher frequencies of multiple mapping than equivalent human RRBS libraries. The zebrafish methylome indicates there is higher global DNA methylation in the zebrafish genome compared with its equivalent human methylome. This observation was confirmed by RRBS of zebrafish liver. High coverage CpG dinucleotides are enriched in CpG island shores more than in the CpG island core. We found that 45% of the mapped CpGs reside in gene bodies, and 7% in gene promoters. This analysis provides a roadmap for generating reproducible base-pair level methylomes for zebrafish using RRBS and our results provide the first evidence that RRBS is a suitable technique for global methylation analysis in zebrafish.

A homologue of retroviral pseudoproteases in the parapoxvirus, orf virus

The nucleotide sequence of a near-terminal region of orf virus DNA was determined. Examination of the sequence revealed an open reading frame encoding a peptide with significant amino acid homology to the pseudoprotease domains recently identified in a number of retroviruses including mouse mammary tumor virus, simian Mason-Pfizer virus, maedi-visna virus, and equine infectious anaemia virus. The orf virus pseudoprotease shares up to 28% amino acid homology with retroviral pseudoproteases and appears to be a discrete transcriptional unit rather than a subunit of a larger polypeptide as is the case in retroviruses. The sharing of amino acid composition across such wide taxonomic boundaries suggests that this polypeptide has a functional significance in both retroviruses and poxviruses.