Genomes of 'Candidatus Liberibacter solanacearum' Haplotype A from New Zealand and the United States Suggest Significant Genome Plasticity in the Species.

'Candidatus Liberibacter solanacearum' contains two solanaceous crop-infecting haplotypes, A and B. Two haplotype A draft genomes were assembled and compared with ZC1 (haplotype B), revealing inversion and relocation genomic rearrangements, numerous single-nucleotide polymorphisms, and differences in phage-related regions. Differences in prophage location and sequence were seen both within and between haplotype comparisons. OrthoMCL and BLAST analyses identified 46 putative coding sequences present in haplotype A that were not present in haplotype B. Thirty-eight of these loci were not found in sequences from other Liberibacter spp. Quantitative polymerase chain reaction (qPCR) assays designed to amplify sequences from 15 of these loci were screened against a panel of 'Ca. L. solanacearum'-positive samples to investigate genetic diversity. Seven of the assays demonstrated within-haplotype diversity; five failed to amplify loci in at least one haplotype A sample while three assays produced amplicons from some haplotype B samples. Eight of the loci assays showed consistent A-B differentiation. Differences in genome arrangements, prophage, and qPCR results suggesting locus diversity within the haplotypes provide more evidence for genetic complexity in this emerging bacterial species.

Structure and expression of GSL1 and GSL2 genes encoding gibberellin stimulated-like proteins in diploid and highly heterozygous tetraploid potato reveals their highly conserved and essential status.

BACKGROUND: GSL1 and GSL2, Gibberellin Stimulated-Like proteins (also known as Snakin-1 and Snakin-2), are cysteine-rich peptides from potato (Solanum tuberosum L.) with antimicrobial properties. Similar peptides in other species have been implicated in diverse biological processes and are hypothesised to play a role in several aspects of plant development, plant responses to biotic or abiotic stress through their participation in hormone crosstalk, and redox homeostasis. To help resolve the biological roles of GSL1 and GSL2 peptides we have undertaken an in depth analysis of the structure and expression of these genes in potato. RESULTS: We have characterised the full length genes for both GSL1 (chromosome 4) and GSL2 (chromosome 1) from diploid and tetraploid potato using the reference genome sequence of potato, coupled with further next generation sequencing of four highly heterozygous tetraploid cultivars. The frequency of SNPs in GSL1 and GSL2 were very low with only one SNP every 67 and 53 nucleotides in exon regions of GSL1 and GSL2, respectively. Analysis of comprehensive RNA-seq data substantiated the role of specific promoter motifs in transcriptional control of gene expression. Expression analysis based on the frequency of next generation sequence reads established that GSL2 was expressed at a higher level than GSL1 in 30 out of 32 tissue and treatment libraries. Furthermore, both the GSL1 and GSL2 genes exhibited constitutive expression that was not up regulated in response to biotic or abiotic stresses, hormone treatments or wounding. Potato transformation with antisense knock-down expression cassettes failed to recover viable plants. CONCLUSIONS: The potato GSL1 and GSL2 genes are very highly conserved suggesting they contribute to an important biological function. The known antimicrobial activity of the GSL proteins, coupled with the FPKM analysis from RNA-seq data, implies that both genes contribute to the constitutive defence barriers in potatoes. The lethality of antisense knock-down expression of GSL1 and GSL2, coupled with the rare incidence of SNPs in these genes, suggests an essential role for this gene family. These features are consistent with the GSL protein family playing a role in several aspects of plant development in addition to plant defence against biotic stresses.

In plants, highly expressed genes are the least compact.

In both the monocot rice and the dicot Arabidopsis, highly expressed genes have more and longer introns and a larger primary transcript than genes expressed at a low level: higher expressed genes tend to be less compact than lower expressed genes. In animal genomes, it is the other way round. Although the length differences in plant genes are much smaller than in animals, these findings indicate that plant genes are in this respect different from animal genes. Explanations for the relationship between gene configuration and gene expression in animals might be (or might have been) less important in plants. We speculate that selection, if any, on genome configuration has taken a different turn after the divergence of plants and animals.

High-throughput bioinformatics with the Cyrille2 pipeline system.

BACKGROUND: Modern omics research involves the application of high-throughput technologies that generate vast volumes of data. These data need to be pre-processed, analyzed and integrated with existing knowledge through the use of diverse sets of software tools, models and databases. The analyses are often interdependent and chained together to form complex workflows or pipelines. Given the volume of the data used and the multitude of computational resources available, specialized pipeline software is required to make high-throughput analysis of large-scale omics datasets feasible. RESULTS: We have developed a generic pipeline system called Cyrille2. The system is modular in design and consists of three functionally distinct parts: 1) a web based, graphical user interface (GUI) that enables a pipeline operator to manage the system; 2) the Scheduler, which forms the functional core of the system and which tracks what data enters the system and determines what jobs must be scheduled for execution, and; 3) the Executor, which searches for scheduled jobs and executes these on a compute cluster. CONCLUSION: The Cyrille2 system is an extensible, modular system, implementing the stated requirements. Cyrille2 enables easy creation and execution of high throughput, flexible bioinformatics pipelines.

Mapping gene regulatory networks from single-cell omics data.

Single-cell techniques are advancing rapidly and are yielding unprecedented insight into cellular heterogeneity. Mapping the gene regulatory networks (GRNs) underlying cell states provides attractive opportunities to mechanistically understand this heterogeneity. In this review, we discuss recently emerging methods to map GRNs from single-cell transcriptomics data, tackling the challenge of increased noise levels and data sparsity compared with bulk data, alongside increasing data volumes. Next, we discuss how new techniques for single-cell epigenomics, such as single-cell ATAC-seq and single-cell DNA methylation profiling, can be used to decipher gene regulatory programmes. We finally look forward to the application of single-cell multi-omics and perturbation techniques that will likely play important roles for GRN inference in the future.

What drives plant stress genes?

Currently, there is a lot of interest in the plant stress response. Using large-scale genomics approaches, more and more genes are being identified that are involved in or even regulate this complex process. The recent boost in expression profile analyses for several plant stress responses has enabled the identification of new promoter elements as important factors in establishing the expression regulatory network controlling plant stress response.

Draft Genome Sequences of Three Pectobacterium Strains Causing Blackleg of Potato: P. carotovorum subsp. brasiliensis ICMP 19477, P. atrosepticum ICMP 1526, and P. carotovorum subsp. carotovorum UGC32.

Blackleg is a disease caused by several species of Pectobacterium that results in losses to potato crops worldwide. Here, we report the draft genomes of three taxonomically and geographically distinct blackleg-causing strains of Pectobacterium: P. carotovorum subsp. brasiliensis ICMP 19477, P. atrosepticum ICMP 1526, and P. carotovorum subsp. carotovorum UGC32. Comparison of these genomes will support the identification of common traits associated with their capacity to cause blackleg.

Allermatch, a webtool for the prediction of potential allergenicity according to current FAO/WHO Codex alimentarius guidelines.

BACKGROUND: Novel proteins entering the food chain, for example by genetic modification of plants, have to be tested for allergenicity. Allermatch http://allermatch.org is a webtool for the efficient and standardized prediction of potential allergenicity of proteins and peptides according to the current recommendations of the FAO/WHO Expert Consultation, as outlined in the Codex alimentarius. DESCRIPTION: A query amino acid sequence is compared with all known allergenic proteins retrieved from the protein databases using a sliding window approach. This identifies stretches of 80 amino acids with more than 35% similarity or small identical stretches of at least six amino acids. The outcome of the analysis is presented in a concise format. The predictive performance of the FAO/WHO criteria is evaluated by screening sets of allergens and non-allergens against the Allermatch databases. Besides correct predictions, both methods are shown to generate false positive and false negative hits and the outcomes should therefore be combined with other methods of allergenicity assessment, as advised by the FAO/WHO. CONCLUSIONS: Allermatch provides an accessible, efficient, and useful webtool for analysis of potential allergenicity of proteins introduced in genetically modified food prior to market release that complies with current FAO/WHO guidelines.

Identification of Mendel's white flower character.

BACKGROUND: The genetic regulation of flower color has been widely studied, notably as a character used by Mendel and his predecessors in the study of inheritance in pea. METHODOLOGY/PRINCIPAL FINDINGS: We used the genome sequence of model legumes, together with their known synteny to the pea genome to identify candidate genes for the A and A2 loci in pea. We then used a combination of genetic mapping, fast neutron mutant analysis, allelic diversity, transcript quantification and transient expression complementation studies to confirm the identity of the candidates. CONCLUSIONS/SIGNIFICANCE: We have identified the pea genes A and A2. A is the factor determining anthocyanin pigmentation in pea that was used by Gregor Mendel 150 years ago in his study of inheritance. The A gene encodes a bHLH transcription factor. The white flowered mutant allele most likely used by Mendel is a simple G to A transition in a splice donor site that leads to a mis-spliced mRNA with a premature stop codon, and we have identified a second rare mutant allele. The A2 gene encodes a WD40 protein that is part of an evolutionarily conserved regulatory complex.

Post alignment clustering procedure for comparative quantitative proteomics LC-MS data.

Comparative LC-MS is a powerful method for detailed quantitative comparison of complex protein mixtures. Dedicated software is required for detection, matching, and alignment of peaks in multiple LC-MS datasets. However, retention time shifts, saturation effects, limitations of experimental accuracy, and possible occurrence of split peaks make it difficult for software to perfectly match all chromatograms. We describe a procedure to assess the above problems and show that dataset quality can be enhanced with the aid of cluster analysis.

DNAVis: interactive visualization of comparative genome annotations.

The software package DNAVis offers a fast, interactive and real-time visualization of DNA sequences and their comparative genome annotations. DNAVis implements advanced methods of information visualization such as linked views, perspective walls and semantic zooming, in addition to the display of heterologous data in dot plot-like matrix views.

Local coexpression domains of two to four genes in the genome of Arabidopsis.

Expression of genes in eukaryotic genomes is known to cluster, but cluster size is generally loosely defined and highly variable. We have here taken a very strict definition of cluster as sets of physically adjacent genes that are highly coexpressed and form so-called local coexpression domains. The Arabidopsis (Arabidopsis thaliana) genome was analyzed for the presence of such local coexpression domains to elucidate its functional characteristics. We used expression data sets that cover different experimental conditions, organs, tissues, and cells from the Massively Parallel Signature Sequencing repository and microarray data (Affymetrix) from a detailed root analysis. With these expression data, we identified 689 and 1,481 local coexpression domains, respectively, consisting of two to four genes with a pairwise Pearson's correlation coefficient larger than 0.7. This number is approximately 1- to 5-fold higher than the numbers expected by chance. A small (5%-10%) yet significant fraction of genes in the Arabidopsis genome is therefore organized into local coexpression domains. These local coexpression domains were distributed over the genome. Genes in such local domains were for the major part not categorized in the same functional category (GOslim). Neither tandemly duplicated genes nor shared promoter sequence nor gene distance explained the occurrence of coexpression of genes in such chromosomal domains. This indicates that other parameters in genes or gene positions are important to establish coexpression in local domains of Arabidopsis chromosomes.

Complete genome sequence of Lactobacillus plantarum WCFS1.

The 3,308,274-bp sequence of the chromosome of Lactobacillus plantarum strain WCFS1, a single colony isolate of strain NCIMB8826 that was originally isolated from human saliva, has been determined, and contains 3,052 predicted protein-encoding genes. Putative biological functions could be assigned to 2,120 (70%) of the predicted proteins. Consistent with the classification of L. plantarum as a facultative heterofermentative lactic acid bacterium, the genome encodes all enzymes required for the glycolysis and phosphoketolase pathways, all of which appear to belong to the class of potentially highly expressed genes in this organism, as was evident from the codon-adaptation index of individual genes. Moreover, L. plantarum encodes a large pyruvate-dissipating potential, leading to various end-products of fermentation. L. plantarum is a species that is encountered in many different environmental niches, and this flexible and adaptive behavior is reflected by the relatively large number of regulatory and transport functions, including 25 complete PTS sugar transport systems. Moreover, the chromosome encodes >200 extracellular proteins, many of which are predicted to be bound to the cell envelope. A large proportion of the genes encoding sugar transport and utilization, as well as genes encoding extracellular functions, appear to be clustered in a 600-kb region near the origin of replication. Many of these genes display deviation of nucleotide composition, consistent with a foreign origin. These findings suggest that these genes, which provide an important part of the interaction of L. plantarum with its environment, form a lifestyle adaptation region in the chromosome.