Genome Analysis of the Broad Host Range Necrotroph Nalanthamala psidii Highlights Genes Associated With Virulence.

Guava wilt disease is caused by the fungus Nalanthamala psidii. The wilt disease results in large-scale destruction of orchards in South Africa, Taiwan, and several Southeast Asian countries. De novo assembly, annotation, and in-depth analysis of the N. psidii genome were carried out to facilitate the identification of characteristics associated with pathogenicity and pathogen evolution. The predicted secretome revealed a range of CAZymes, proteases, lipases and peroxidases associated with plant cell wall degradation, nutrient acquisition, and disease development. Further analysis of the N. psidii carbohydrate-active enzyme profile exposed the broad-spectrum necrotrophic lifestyle of the pathogen, which was corroborated by the identification of putative effectors and secondary metabolites with the potential to induce tissue necrosis and cell surface-dependent immune responses. Putative regulatory proteins including transcription factors and kinases were identified in addition to transporters potentially involved in the secretion of secondary metabolites. Transporters identified included important ABC and MFS transporters involved in the efflux of fungicides. Analysis of the repetitive landscape and the detection of mechanisms linked to reproduction such as het and mating genes rendered insights into the biological complexity and evolutionary potential of N. psidii as guava pathogen. Hence, the assembly and annotation of the N. psidii genome provided a valuable platform to explore the pathogenic potential and necrotrophic lifestyle of the guava wilt pathogen.

Sialotranscriptomics of Rhipicephalus zambeziensis reveals intricate expression profiles of secretory proteins and suggests tight temporal transcriptional regulation during blood-feeding.

BACKGROUND: Ticks secrete a diverse mixture of secretory proteins into the host to evade its immune response and facilitate blood-feeding, making secretory proteins attractive targets for the production of recombinant anti-tick vaccines. The largely neglected tick species, Rhipicephalus zambeziensis, is an efficient vector of Theileria parva in southern Africa but its available sequence information is limited. Next generation sequencing has advanced sequence availability for ticks in recent years and has assisted the characterisation of secretory proteins. This study focused on the de novo assembly and annotation of the salivary gland transcriptome of R. zambeziensis and the temporal expression of secretory protein transcripts in female and male ticks, before the onset of feeding and during early and late feeding. RESULTS: The sialotranscriptome of R. zambeziensis yielded 23,631 transcripts from which 13,584 non-redundant proteins were predicted. Eighty-six percent of these contained a predicted start and stop codon and were estimated to be putatively full-length proteins. A fifth (2569) of the predicted proteins were annotated as putative secretory proteins and explained 52% of the expression in the transcriptome. Expression analyses revealed that 2832 transcripts were differentially expressed among feeding time points and 1209 between the tick sexes. The expression analyses further indicated that 57% of the annotated secretory protein transcripts were differentially expressed. Dynamic expression profiles of secretory protein transcripts were observed during feeding of female ticks. Whereby a number of transcripts were upregulated during early feeding, presumably for feeding site establishment and then during late feeding, 52% of these were downregulated, indicating that transcripts were required at specific feeding stages. This suggested that secretory proteins are under stringent transcriptional regulation that fine-tunes their expression in salivary glands during feeding. No open reading frames were predicted for 7947 transcripts. This class represented 17% of the differentially expressed transcripts, suggesting a potential transcriptional regulatory function of long non-coding RNA in tick blood-feeding. CONCLUSIONS: The assembled sialotranscriptome greatly expands the sequence availability of R. zambeziensis, assists in our understanding of the transcription of secretory proteins during blood-feeding and will be a valuable resource for future vaccine candidate selection.

High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development.

Using the latest sequencing and optical mapping technologies, we have produced a high-quality de novo assembly of the apple (Malus domestica Borkh.) genome. Repeat sequences, which represented over half of the assembly, provided an unprecedented opportunity to investigate the uncharacterized regions of a tree genome; we identified a new hyper-repetitive retrotransposon sequence that was over-represented in heterochromatic regions and estimated that a major burst of different transposable elements (TEs) occurred 21 million years ago. Notably, the timing of this TE burst coincided with the uplift of the Tian Shan mountains, which is thought to be the center of the location where the apple originated, suggesting that TEs and associated processes may have contributed to the diversification of the apple ancestor and possibly to its divergence from pear. Finally, genome-wide DNA methylation data suggest that epigenetic marks may contribute to agronomically relevant aspects, such as apple fruit development.

Genome-wide association analysis of agronomic traits in wheat under drought-stressed and non-stressed conditions.

This study determined the population structure and genome-wide marker-trait association of agronomic traits of wheat for drought-tolerance breeding. Ninety-three diverse bread wheat genotypes were genotyped using the Diversity Arrays Technology sequencing (DArTseq) protocol. The number of days-to-heading (DTH), number of days-to-maturity (DTM), plant height (PHT), spike length (SPL), number of kernels per spike (KPS), thousand kernel weight (TKW) and grain yield (GYLD), assessed under drought-stressed and non-stressed conditions, were considered for the study. Population structure analysis and genome-wide association mapping were undertaken based on 16,383 silico DArTs loci with < 10% missing data. The population evaluated was grouped into nine distinct genetic structures. Inter-chromosomal linkage disequilibrium showed the existence of linkage decay as physical distance increased. A total of 62 significant (P < 0.001) marker-trait associations (MTAs) were detected explaining more than 20% of the phenotypic variation observed under both drought-stressed and non-stressed conditions. Significant (P < 0.001) MTA event(s) were observed for DTH, PHT, SPL, SPS, and KPS; under both stressed and non-stressed conditions, while additional significant (P < 0.05) associations were observed for TKW, DTM and GYLD under non-stressed condition. The MTAs reported in this population could be useful to initiate marker-assisted selection (MAS) and targeted trait introgression of wheat under drought-stressed and non-stressed conditions, and for fine mapping and cloning of the underlying genes and QTL.

De novo assembly and annotation of the salivary gland transcriptome of Rhipicephalus appendiculatus male and female ticks during blood feeding.

Tick secretory proteins modulate haemostasis, inflammation and immune responses of the host and are attractive recombinant anti-tick vaccine candidates. Yet, many of the proteins have not been characterised due to the limited sequence availability for ticks and other arthropods for homology-based annotation. To address this limitation, we sequenced the salivary glands of the economically important adult male and female Rhipicephalus appendiculatus ticks during feeding. The quality-filtered Illumina sequencing reads were de novo assembled to generate a R. appendiculatus sialotranscriptome of 21,410 transcripts. A non-redundant set of 12,761 R. appendiculatus proteins was predicted from the transcripts, including 2134 putative secretory and 8237 putative housekeeping proteins. Secretory proteins accounted for most of the expression in the salivary gland transcriptome (63%). Of the secretory protein class, the Glycine-rich superfamily contributed 66% and the Lipocalin family 12% of the transcriptome expression. Differential expression analysis identified 1758 female and 2346 male up-regulated transcripts, suggesting varying blood-feeding mechanisms employed between female and male ticks. The sialotranscriptome assembled in this work, greatly improves on the sequence information available for R. appendiculatus and is a valuable resource for potential future vaccine candidate selection.

Phylogenomic analysis reveals deep divergence and recombination in an economically important grapevine virus.

The evolutionary history of the exclusively grapevine (Vitis spp.) infecting, grapevine leafroll-associated virus 3 (GLRaV-3) has not been studied extensively, partly due to limited available sequence data. In this study we trace the evolutionary history of GLRaV-3, focussing on isolate GH24, a newly discovered variant. GH24 was discovered through the use of next-generation sequencing (NGS) and the whole genome sequence determined and validated with Sanger sequencing. We assembled an alignment of all 13 available whole genomes of GLRaV-3 isolates and all other publicly available GLRaV-3 sequence data. Using multiple recombination detection methods we identified a clear signal for recombination in one whole genome sequence and further evidence for recombination in two more, including GH24. We inferred phylogenetic trees and networks and estimated the ages of common ancestors of GLRaV-3 clades by means of relaxed clock models calibrated with asynchronous sampling dates. Our results generally confirm previously identified variant groups as well as two new groups (VII and VIII). Higher order groups were defined as supergroups designated A to D. Supergroup A includes variant groups I-V and supergroup B group VI and its related unclassified isolates. Supergroups C and D are less well known, including the newly identified groups VII (including isolate GH24) and VIII respectively. The inferred node ages suggest that the origins of the major groups of GLRaV-3, including isolate GH24, may have occurred prior to worldwide cultivation of grapevines, whilst the current diversity represents closely related isolates that diverged from common ancestors within the last century.

Heteroditopic P,N ligands in gold(I) complexes: synthesis, structure and cytotoxicity.

New heteroditopic, bi- and multidentate imino- and aminophosphine ligands were synthesised and complexed to [AuCl(THT)] (THT=tetrahydrothiophene). X-ray crystallography confirmed Schiff base formation in three products, the successful reduction of the imino-group to the sp(3)-hybridised amine in several instances, and confirmed the formation of mono-gold(I) imino- and aminophosphine complexes for four Au-complexes. Cytotoxicity studies in cancerous and non-cancerous cell lines showed a marked increase in cytotoxicity upon ligand complexation to gold(I). These findings were supported by results from the 60-cell line fingerprint screen of the Developmental Therapeutics Programme of the National Institutes of Health for two promising compounds. The cytotoxicity of some of these ligands and gold(I)complexes is due to the induction of apoptosis. The ligands and gold(I)complexes demonstrated selective toxicity towards specific cell lines, with Jurkat T cells being more sensitive to the cytotoxic effects of these compounds, while the non-cancerous human cell line KMST6 proved more resistant when compared to the cancerous cell lines. Results from the NIH DTP 60 cell-line fingerprint screen support the observed enhancement of cytotoxicity upon gold(I) complexation. One gold(I)complex induced high levels of apoptosis at concentrations of 50 µM in all the cell lines screened in this study, while some of the other compounds selectively induced apoptosis in the cell lines. These results point towards the potential for selective toxicity to cancerous cells through the induction of apoptosis.

High-throughput sequencing reveals small RNAs involved in ASGV infection.

BACKGROUND: Plant small RNAs (sRNAs) associated with virulent virus infections have been reported by previous studies, while the involvement of sRNAs in latent virus infection remains largely uncharacterised. Apple trees show a high degree of resistance and tolerance to viral infections. We analysed two sRNA deep sequencing datasets, prepared from different RNA size fractions, to identify sRNAs involved in Apple stem grooving virus (ASGV) infection. RESULTS: sRNA analysis revealed virus-derived siRNAs (vsiRNAs) originating from two ASGV genetic variants. A vsiRNA profile for one of the ASGV variants was also generated showing an increase in siRNA production towards the 3' end of the virus genome. Virus-derived sRNAs longer than those previously analysed were also observed in the sequencing data. Additionally, tRNA-derived sRNAs were identified and characterised. These sRNAs covered a broad size-range and originated from both ends of the mature tRNAs as well as from their central regions. Several tRNA-derived sRNAs showed differential regulation due to ASGV infection. No changes in microRNA, natural-antisense transcript siRNA, phased-siRNA and repeat-associated siRNA levels were observed. CONCLUSIONS: This study is the first report on the apple sRNA-response to virus infection. The results revealed the vsiRNAs profile of an ASGV variant, as well as the alteration of the tRNA-derived sRNA profile in response to latent virus infection. It also highlights the importance of library preparation in the interpretation of high-throughput sequencing data.

Extending the sRNAome of apple by next-generation sequencing.

The global importance of apple as a fruit crop necessitates investigations into molecular aspects of the processes that influence fruit quality and yield, including plant development, fruit ripening and disease resistance. In order to study and understand biological processes it is essential to recognise the range of molecules, which influence these processes. Small non-coding RNAs are regulatory agents involved in diverse plant activities, ranging from development to stress response. The occurrence of these molecules in apple leaves was studied by means of next-generation sequencing. 85 novel microRNA (miRNA) gene loci were predicted and characterized along with known miRNA loci. Both cis- and trans-natural antisense transcript pairs were identified. Although the trans-overlapping regions were enriched in small RNA (sRNA) production, cis-overlaps did not seem to agree. More than 150 phased regions were also identified, and for a small subset of these, potential miRNAs that could initiate phasing, were revealed. Repeat-associated siRNAs, which are generated from repetitive genomic regions such as transposons, were also analysed. For this group almost all available repeat sequences, associated with the apple genome and present in Repbase, were found to produce siRNAs. Results from this study extend our current knowledge on apple sRNAs and their precursors significantly. A rich molecular resource has been created and is available to the research community to serve as a baseline for future studies.

A glance at quality score: implication for de novo transcriptome reconstruction of Illumina reads.

Downstream analyses of short-reads from next-generation sequencing platforms are often preceded by a pre-processing step that removes uncalled and wrongly called bases. Standard approaches rely on their associated base quality scores to retain the read or a portion of it when the score is above a predefined threshold. It is difficult to differentiate sequencing error from biological variation without a reference using quality scores. The effects of quality score based trimming have not been systematically studied in de novo transcriptome assembly. Using RNA-Seq data produced from Illumina, we teased out the effects of quality score based filtering or trimming on de novo transcriptome reconstruction. We showed that assemblies produced from reads subjected to different quality score thresholds contain truncated and missing transfrags when compared to those from untrimmed reads. Our data supports the fact that de novo assembling of untrimmed data is challenging for de Bruijn graph assemblers. However, our results indicates that comparing the assemblies from untrimmed and trimmed read subsets can suggest appropriate filtering parameters and enable selection of the optimum de novo transcriptome assembly in non-model organisms.

Diversity arrays technology (DArT) markers in apple for genetic linkage maps.

Diversity Arrays Technology (DArT) provides a high-throughput whole-genome genotyping platform for the detection and scoring of hundreds of polymorphic loci without any need for prior sequence information. The work presented here details the development and performance of a DArT genotyping array for apple. This is the first paper on DArT in horticultural trees. Genetic mapping of DArT markers in two mapping populations and their integration with other marker types showed that DArT is a powerful high-throughput method for obtaining accurate and reproducible marker data, despite the low cost per data point. This method appears to be suitable for aligning the genetic maps of different segregating populations. The standard complexity reduction method, based on the methylation-sensitive PstI restriction enzyme, resulted in a high frequency of markers, although there was 52-54% redundancy due to the repeated sampling of highly similar sequences. Sequencing of the marker clones showed that they are significantly enriched for low-copy, genic regions. The genome coverage using the standard method was 55-76%. For improved genome coverage, an alternative complexity reduction method was examined, which resulted in less redundancy and additional segregating markers. The DArT markers proved to be of high quality and were very suitable for genetic mapping at low cost for the apple, providing moderate genome coverage. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-011-9579-5) contains supplementary material, which is available to authorized users.

Solution structure of RING finger-like domain of retinoblastoma-binding protein-6 (RBBP6) suggests it functions as a U-box.

Retinoblastoma-binding protein-6 (RBBP6) plays a facilitating role, through its RING finger-like domain, in the ubiquitination of p53 by Hdm2 that is suggestive of E4-like activity. Although the presence of eight conserved cysteine residues makes it highly probable that the RING finger-like domain coordinates two zinc ions, analysis of the primary sequence suggests an alternative classification as a member of the U-box family, the members of which do not bind zinc ions. We show here that despite binding two zinc ions, the domain adopts a homodimeric structure highly similar to those of a number of U-boxes. Zinc ions could be replaced by cadmium ions without significantly disrupting the structure or the stability of the domain, although the rate of substitution was an order of magnitude slower than any previous measurement, suggesting that the structure is particularly stable, a conclusion supported by the high thermal stability of the domain. A hallmark of U-box-containing proteins is their association with chaperones, with which they cooperate in eliminating irretrievably unfolded proteins by tagging them for degradation by the proteasome. Using a yeast two-hybrid screen, we show that RBBP6 interacts with chaperones Hsp70 and Hsp40 through its N-terminal ubiquitin-like domain. Taken together with the structural similarities to U-box-containing proteins, our data suggest that RBBP6 plays a role in chaperone-mediated ubiquitination and possibly in protein quality control.

The genome of woodland strawberry (Fragaria vesca).

The woodland strawberry, Fragaria vesca (2n = 2x = 14), is a versatile experimental plant system. This diminutive herbaceous perennial has a small genome (240 Mb), is amenable to genetic transformation and shares substantial sequence identity with the cultivated strawberry (Fragaria × ananassa) and other economically important rosaceous plants. Here we report the draft F. vesca genome, which was sequenced to ×39 coverage using second-generation technology, assembled de novo and then anchored to the genetic linkage map into seven pseudochromosomes. This diploid strawberry sequence lacks the large genome duplications seen in other rosids. Gene prediction modeling identified 34,809 genes, with most being supported by transcriptome mapping. Genes critical to valuable horticultural traits including flavor, nutritional value and flowering time were identified. Macrosyntenic relationships between Fragaria and Prunus predict a hypothetical ancestral Rosaceae genome that had nine chromosomes. New phylogenetic analysis of 154 protein-coding genes suggests that assignment of Populus to Malvidae, rather than Fabidae, is warranted.

Deep sequencing analysis of viruses infecting grapevines: Virome of a vineyard.

Double stranded RNA, isolated from 44 pooled randomly selected vines from a diseased South African vineyard, has been used in a deep sequencing analysis to build a census of the viral population. The dsRNA was sequenced in an unbiased manner using the sequencing-by-synthesis technology offered by the Illumina Genome Analyzer II and yielded 837 megabases of metagenomic sequence data. Four known viral pathogens were identified. It was found that Grapevine leafroll-associated virus 3 (GLRaV-3) is the most prevalent species, constituting 59% of the total reads, followed by Grapevine rupestris stem pitting-associated virus and Grapevine virus A. Grapevine virus E, a virus not previously reported in South African vineyards, was identified in the census. Viruses not previously identified in grapevine were also detected. The second most prevalent virus detected was a member of the Chrysoviridae family similar to Penicillium chrysogenum virus. Sequences aligning to two other mycoviruses were also detected.

Genome-wide SNP identification by high-throughput sequencing and selective mapping allows sequence assembly positioning using a framework genetic linkage map.

BACKGROUND: Determining the position and order of contigs and scaffolds from a genome assembly within an organism's genome remains a technical challenge in a majority of sequencing projects. In order to exploit contemporary technologies for DNA sequencing, we developed a strategy for whole genome single nucleotide polymorphism sequencing allowing the positioning of sequence contigs onto a linkage map using the bin mapping method. RESULTS: The strategy was tested on a draft genome of the fungal pathogen Venturia inaequalis, the causal agent of apple scab, and further validated using sequence contigs derived from the diploid plant genome Fragaria vesca. Using our novel method we were able to anchor 70% and 92% of sequences assemblies for V. inaequalis and F. vesca, respectively, to genetic linkage maps. CONCLUSIONS: We demonstrated the utility of this approach by accurately determining the bin map positions of the majority of the large sequence contigs from each genome sequence and validated our method by mapping single sequence repeat markers derived from sequence contigs on a full mapping population.

In vitro evaluation of dichloro-bis(pyrazole)palladium(II) and dichloro-bis(pyrazole)platinum(II) complexes as anticancer agents.

INTRODUCTION: Cisplatin (cis-diamminedichloroplatinum) was first identified for its anti-bacterial activity, and was later also shown to be an efficient anticancer agent. However, the therapeutic use of this anticancer drug is somewhat limited by its toxic side effects, which include nephrotoxicity, nausea, and vomiting. Furthermore the development of drug-resistant tumours is commonly observed following therapy with cisplatin. Hence there is a need for improved platinum derived drugs to overcome these limitations. AIMS: Apoptosis contributes significantly to the cytotoxic effects of anticancer agents such as cisplatin; therefore in this study the potential anticancer properties of a series of pyrazole palladium(II) and platinum(II) complexes, [(3,5-R(2)pz)(2)PdCl(2)] [R = H (1), R = Me (2)] and [(3,5-R(2)pz)(2)PtCl(2)] [R = H (3), R = Me (4)], were evaluated by assessment of their pro-apoptotic activity. METHODS: The induction of apoptosis was measured in CHO cells by the detection of phosphatidylserine (PS) exposure using the annexin V and APOPercentage assays; DNA fragmentation using the Terminal deoxynucleotide transferase dUTP Nick End Labelling (TUNEL) assay; and the detection of activated caspase-3. RESULTS: The platinum complexes were shown to be considerably more active than the palladium complexes, with complex 3 demonstrating the highest level of cytotoxic and pro-apoptotic activity. The LD(50) values for complex 3 and cisplatin were 20 and 70 microM, respectively, demonstrating that the cytotoxic activity for complex 3 was three times higher than for cisplatin. Various human cancer cell lines, including CaSki, HeLa, as well as the p53 mutant Jurkat T cell line were also shown to be susceptible to complex 3. CONCLUSIONS: Collectively, this in vitro study provides insights into action of palladium and platinum complexes and demonstrates the potential use of these compounds, and in particular complex 3, in the development of new anticancer agents.

DWNN, a novel ubiquitin-like domain, implicates RBBP6 in mRNA processing and ubiquitin-like pathways.

BACKGROUND: RBBP6 is a 250 kDa splicing-associated protein that has been identified as an E3 ligase due to the presence of a RING finger domain. In humans and mice it interacts with both p53 and Rb, and plays a role in the induction of apoptosis and regulation of the cell cycle. RBBP6 has recently been shown to be highly up-regulated in oesophageal cancer, and to be a promising target for immunotherapy against the disease. RESULTS: We show here using heteronuclear NMR that the N-terminal 81 amino acids of RBBP6 constitute a novel ubiquitin-like domain, which we have called the DWNN domain. The domain lacks conserved equivalents of K48 and K63, although the equivalents of K6 and K29 are highly, although not absolutely, conserved. The di-glycine motif that is characteristic of proteins involved in ubiquitination is found in the human and mouse form of the domain, although it is not present in all organisms. It forms part of a three-domain form of RBBP6 containing the DWNN domain, a zinc knuckle and a RING finger domain, which is found in all eukaryotic genomes so far examined, in the majority of cases at single copy number. The domain is also independently expressed in vertebrates as a single domain protein. CONCLUSION: DWNN is a novel ubiquitin-like domain found only at the N-terminus of the RBBP6 family of splicing-associated proteins. The ubiquitin-like structure of the domain greatly increases the likelihood that RBBP6 functions through some form of ubiquitin-like modification. Furthermore, the fact that the DWNN domain is independently expressed in higher vertebrates leads us to propose that the domain may itself function as a novel ubiquitin-like modifier of other proteins.