The genome of Tetranychus urticae reveals herbivorous pest adaptations.

The spider mite Tetranychus urticae is a cosmopolitan agricultural pest with an extensive host plant range and an extreme record of pesticide resistance. Here we present the completely sequenced and annotated spider mite genome, representing the first complete chelicerate genome. At 90 megabases T. urticae has the smallest sequenced arthropod genome. Compared with other arthropods, the spider mite genome shows unique changes in the hormonal environment and organization of the Hox complex, and also reveals evolutionary innovation of silk production. We find strong signatures of polyphagy and detoxification in gene families associated with feeding on different hosts and in new gene families acquired by lateral gene transfer. Deep transcriptome analysis of mites feeding on different plants shows how this pest responds to a changing host environment. The T. urticae genome thus offers new insights into arthropod evolution and plant-herbivore interactions, and provides unique opportunities for developing novel plant protection strategies.

Gene silencing pathways found in the green alga Volvox carteri reveal insights into evolution and origins of small RNA systems in plants.

BACKGROUND: Volvox carteri (V. carteri) is a multicellular green alga used as model system for the evolution of multicellularity. So far, the contribution of small RNA pathways to these phenomena is not understood. Thus, we have sequenced V. carteri Argonaute 3 (VcAGO3)-associated small RNAs from different developmental stages. RESULTS: Using this functional approach, we define the Volvox microRNA (miRNA) repertoire and show that miRNAs are not conserved in the closely related unicellular alga Chlamydomonas reinhardtii. Furthermore, we find that miRNAs are differentially expressed during different life stages of V. carteri. In addition to miRNAs, transposon-associated small RNAs or phased siRNA loci, which are common in higher land plants, are highly abundant in Volvox as well. Transposons not only give rise to miRNAs and other small RNAs, they are also targets of small RNAs. CONCLUSION: Our analyses reveal a surprisingly complex small RNA network in Volvox as elaborate as in higher land plants. At least the identified VcAGO3-associated miRNAs are not conserved in C. reinhardtii suggesting fast evolution of small RNA systems. Thus, distinct small RNAs may contribute to multicellularity and also division of labor in reproductive and somatic cells.

The cis-regulatory code of Hox function in Drosophila.

Precise gene expression is a fundamental aspect of organismal function and depends on the combinatorial interplay of transcription factors (TFs) with cis-regulatory DNA elements. While much is known about TF function in general, our understanding of their cell type-specific activities is still poor. To address how widely expressed transcriptional regulators modulate downstream gene activity with high cellular specificity, we have identified binding regions for the Hox TF Deformed (Dfd) in the Drosophila genome. Our analysis of architectural features within Hox cis-regulatory response elements (HREs) shows that HRE structure is essential for cell type-specific gene expression. We also find that Dfd and Ultrabithorax (Ubx), another Hox TF specifying different morphological traits, interact with non-overlapping regions in vivo, despite their similar DNA binding preferences. While Dfd and Ubx HREs exhibit comparable design principles, their motif compositions and motif-pair associations are distinct, explaining the highly selective interaction of these Hox proteins with the regulatory environment. Thus, our results uncover the regulatory code imprinted in Hox enhancers and elucidate the mechanisms underlying functional specificity of TFs in vivo.

Trowel: a fast and accurate error correction module for Illumina sequencing reads.

MOTIVATION: The ability to accurately read the order of nucleotides in DNA and RNA is fundamental for modern biology. Errors in next-generation sequencing can lead to many artifacts, from erroneous genome assemblies to mistaken inferences about RNA editing. Uneven coverage in datasets also contributes to false corrections. RESULT: We introduce Trowel, a massively parallelized and highly efficient error correction module for Illumina read data. Trowel both corrects erroneous base calls and boosts base qualities based on the k-mer spectrum. With high-quality k-mers and relevant base information, Trowel achieves high accuracy for different short read sequencing applications.The latency in the data path has been significantly reduced because of efficient data access and data structures. In performance evaluations, Trowel was highly competitive with other tools regardless of coverage, genome size read length and fragment size. AVAILABILITY AND IMPLEMENTATION: Trowel is written in C++ and is provided under the General Public License v3.0 (GPLv3). It is available at http://trowel-ec.sourceforge.net. CONTACT: euncheon.lim@tue.mpg.de or weigel@tue.mpg.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Transcriptome assemblies for studying sex-biased gene expression in the guppy, Poecilia reticulata.

BACKGROUND: Sexually dimorphic phenotypes are generally associated with differential gene expression between the sexes. The study of molecular evolution and genomic location of these differentially expressed, or sex-biased, genes is important for understanding inter-sexual divergence under sex-specific selection pressures. Teleost fish provide a unique opportunity to examine this divergence in the presence of variable sex-determination mechanisms of recent origin. The guppy, Poecilia reticulata, displays sexual dimorphism in size, ornaments, and behavior, traits shaped by natural and sexual selection in the wild. RESULTS: To gain insight into molecular mechanisms underlying the guppy's sexual dimorphism, we assembled a reference transcriptome combining genome-independent as well as genome-guided assemblies and analyzed sex-biased gene expression between different tissues of adult male and female guppies. We found tissue-associated sex-biased expression of genes related to pigmentation, signal transduction, and spermatogenesis in males; and growth, cell-division, extra-cellular matrix organization, nutrient transport, and folliculogenesis in females. While most sex-biased genes were randomly distributed across linkage groups, we observed accumulation of ovary-biased genes on the sex linkage group, LG12. Both testis-biased and ovary-biased genes showed a significantly higher rate of non-synonymous to synonymous substitutions (dN/dS) compared to unbiased genes. However, in somatic tissues only female-biased genes, including those co-expressed in multiple tissues, showed elevated ratios of non-synonymous substitutions. CONCLUSIONS: Our work identifies a set of annotated gene products that are candidate factors affecting sexual dimorphism in guppies. The differential genomic distribution of gonad-biased genes provides evidence for sex-specific selection pressures acting on the nascent sex chromosomes of the guppy. The elevated rates of evolution of testis-biased and female-biased genes indicate differing evolution under distinct selection pressures on the reproductive versus non-reproductive tissues.

A spatial and temporal map of C. elegans gene expression.

The C. elegans genome has been completely sequenced, and the developmental anatomy of this model organism is described at single-cell resolution. Here we utilize strategies that exploit this precisely defined architecture to link gene expression to cell type. We obtained RNAs from specific cells and from each developmental stage using tissue-specific promoters to mark cells for isolation by FACS or for mRNA extraction by the mRNA-tagging method. We then generated gene expression profiles of more than 30 different cells and developmental stages using tiling arrays. Machine-learning-based analysis detected transcripts corresponding to established gene models and revealed novel transcriptionally active regions (TARs) in noncoding domains that comprise at least 10% of the total C. elegans genome. Our results show that about 75% of transcripts with detectable expression are differentially expressed among developmental stages and across cell types. Examination of known tissue- and cell-specific transcripts validates these data sets and suggests that newly identified TARs may exercise cell-specific functions. Additionally, we used self-organizing maps to define groups of coregulated transcripts and applied regulatory element analysis to identify known transcription factor- and miRNA-binding sites, as well as novel motifs that likely function to control subsets of these genes. By using cell-specific, whole-genome profiling strategies, we have detected a large number of novel transcripts and produced high-resolution gene expression maps that provide a basis for establishing the roles of individual genes in cellular differentiation.

Receptor-type guanylate cyclase is required for carbon dioxide sensation by Caenorhabditis elegans.

CO(2) is both a critical regulator of animal physiology and an important sensory cue for many animals for host detection, food location, and mate finding. The free-living soil nematode Caenorhabditis elegans shows CO(2) avoidance behavior, which requires a pair of ciliated sensory neurons, the BAG neurons. Using in vivo calcium imaging, we show that CO(2) specifically activates the BAG neurons and that the CO(2)-sensing function of BAG neurons requires TAX-2/TAX-4 cyclic nucleotide-gated ion channels and the receptor-type guanylate cyclase GCY-9. Our results delineate a molecular pathway for CO(2) sensing and suggest that activation of a receptor-type guanylate cyclase is an evolutionarily conserved mechanism by which animals detect environmental CO(2).

Reference-guided assembly of four diverse Arabidopsis thaliana genomes.

We present whole-genome assemblies of four divergent Arabidopsis thaliana strains that complement the 125-Mb reference genome sequence released a decade ago. Using a newly developed reference-guided approach, we assembled large contigs from 9 to 42 Gb of Illumina short-read data from the Landsberg erecta (Ler-1), C24, Bur-0, and Kro-0 strains, which have been sequenced as part of the 1,001 Genomes Project for this species. Using alignments against the reference sequence, we first reduced the complexity of the de novo assembly and later integrated reads without similarity to the reference sequence. As an example, half of the noncentromeric C24 genome was covered by scaffolds that are longer than 260 kb, with a maximum of 2.2 Mb. Moreover, over 96% of the reference genome was covered by the reference-guided assembly, compared with only 87% with a complete de novo assembly. Comparisons with 2 Mb of dideoxy sequence reveal that the per-base error rate of the reference-guided assemblies was below 1 in 10,000. Our assemblies provide a detailed, genomewide picture of large-scale differences between A. thaliana individuals, most of which are difficult to access with alignment-consensus methods only. We demonstrate their practical relevance in studying the expression differences of polymorphic genes and show how the analysis of sRNA sequencing data can lead to erroneous conclusions if aligned against the reference genome alone. Genome assemblies, raw reads, and further information are accessible through http://1001genomes.org/projects/assemblies.html.

A gene expression map of Arabidopsis thaliana development.

Regulatory regions of plant genes tend to be more compact than those of animal genes, but the complement of transcription factors encoded in plant genomes is as large or larger than that found in those of animals. Plants therefore provide an opportunity to study how transcriptional programs control multicellular development. We analyzed global gene expression during development of the reference plant Arabidopsis thaliana in samples covering many stages, from embryogenesis to senescence, and diverse organs. Here, we provide a first analysis of this data set, which is part of the AtGenExpress expression atlas. We observed that the expression levels of transcription factor genes and signal transduction components are similar to those of metabolic genes. Examining the expression patterns of large gene families, we found that they are often more similar than would be expected by chance, indicating that many gene families have been co-opted for specific developmental processes.

Global effects of the small RNA biogenesis machinery on the Arabidopsis thaliana transcriptome.

In Arabidopsis thaliana, four different dicer-like (DCL) proteins have distinct but partially overlapping functions in the biogenesis of microRNAs (miRNAs) and siRNAs from longer, noncoding precursor RNAs. To analyze the impact of different components of the small RNA biogenesis machinery on the transcriptome, we subjected dcl and other mutants impaired in small RNA biogenesis to whole-genome tiling array analysis. We compared both protein-coding genes and noncoding transcripts, including most pri-miRNAs, in two tissues and several stress conditions. Our analysis revealed a surprising number of common targets in dcl1 and dcl2 dcl3 dcl4 triple mutants. Furthermore, our results suggest that the DCL1 is not only involved in miRNA action but also contributes to silencing of a subset of transposons, apparently through an effect on DNA methylation.

AGRONOMICS1: a new resource for Arabidopsis transcriptome profiling.

Transcriptome profiling has become a routine tool in biology. For Arabidopsis (Arabidopsis thaliana), the Affymetrix ATH1 expression array is most commonly used, but it lacks about one-third of all annotated genes present in the reference strain. An alternative are tiling arrays, but previous designs have not allowed the simultaneous analysis of both strands on a single array. We introduce AGRONOMICS1, a new Affymetrix Arabidopsis microarray that contains the complete paths of both genome strands, with on average one 25mer probe per 35-bp genome sequence window. In addition, the new AGRONOMICS1 array contains all perfect match probes from the original ATH1 array, allowing for seamless integration of the very large existing ATH1 knowledge base. The AGRONOMICS1 array can be used for diverse functional genomics applications such as reliable expression profiling of more than 30,000 genes, detection of alternative splicing, and chromatin immunoprecipitation coupled to microarrays (ChIP-chip). Here, we describe the design of the array and compare its performance with that of the ATH1 array. We find results from both microarrays to be of similar quality, but AGRONOMICS1 arrays yield robust expression information for many more genes, as expected. Analysis of the ATH1 probes on AGRONOMICS1 arrays produces results that closely mirror those of ATH1 arrays. Finally, the AGRONOMICS1 array is shown to be useful for ChIP-chip experiments. We show that heterochromatic H3K9me2 is strongly confined to the gene body of target genes in euchromatic chromosome regions, suggesting that spreading of heterochromatin is limited outside of pericentromeric regions.

Stress-induced changes in the Arabidopsis thaliana transcriptome analyzed using whole-genome tiling arrays.

The responses of plants to abiotic stresses are accompanied by massive changes in transcriptome composition. To provide a comprehensive view of stress-induced changes in the Arabidopsis thaliana transcriptome, we have used whole-genome tiling arrays to analyze the effects of salt, osmotic, cold and heat stress as well as application of the hormone abscisic acid (ABA), an important mediator of stress responses. Among annotated genes in the reference strain Columbia we have found many stress-responsive genes, including several transcription factor genes as well as pseudogenes and transposons that have been missed in previous analyses with standard expression arrays. In addition, we report hundreds of newly identified, stress-induced transcribed regions. These often overlap with known, annotated genes. The results are accessible through the Arabidopsis thaliana Tiling Array Express (At-TAX) homepage, which provides convenient tools for displaying expression values of annotated genes, as well as visualization of unannotated transcribed regions along each chromosome.

Algorithm-Based Liquid Formulation Development Including a DoE Concept Predicts Long-Term Viral Vector Stability.

Specifically tailored amino acid-based formulations were previously shown to have a high potential to avoid stress-mediated degradation of complex molecules such as monoclonal antibodies and viral vectors. By using adenovirus 5 (Ad5) as a model, we studied whether such formulations may also efficiently protect viral vectors in thermal stress experiments and during long-term liquid storage. Algorithm-based amino acid preselection using an excipient database and subsequent application of design of experiments (DoE) in combination with a 37°C challenging model enabled the prediction of long-term storage stability of Ad5. By statistical analysis of the Ad5 infectivity, amino acids with significant influence on Ad5 stability were detected after 2 and 3 weeks of liquid storage at 37°C. Ad5 formulations comprising positively selected amino acids did not reveal any loss of infectivity after 24 months in liquid storage at 5°C. By contrast, a 2 log reduction after 3 months and complete loss of infectivity after 18 months was observed with a standard viral vector formulation. By an optimization round, we designed a simple and well-balanced formulation avoiding MgCl2, previously considered essential in Ad5 formulations. This work demonstrates the efficacy of an algorithm-based development approach in the formulation development for viral vectors.

ESTs and EST-linked polymorphisms for genetic mapping and phylogenetic reconstruction in the guppy, Poecilia reticulata.

BACKGROUND: The guppy, Poecilia reticulata, is a well-known model organism for studying inheritance and variation of male ornamental traits as well as adaptation to different river habitats. However, genomic resources for studying this important model were not previously widely available. RESULTS: With the aim of generating molecular markers for genetic mapping of the guppy, cDNA libraries were constructed from embryos and different adult organs to generate expressed sequence tags (ESTs). About 18,000 ESTs were annotated according to BLASTN and BLASTX results and the sequence information from the 3' UTRs was exploited to generate PCR primers for re-sequencing of genomic DNA from different wild type strains. By comparison of EST-linked genomic sequences from at least four different ecotypes, about 1,700 polymorphisms were identified, representing about 400 distinct genes. Two interconnected MySQL databases were built to organize the ESTs and markers, respectively. A robust phylogeny of the guppy was reconstructed, based on 10 different nuclear genes. CONCLUSION: Our EST and marker databases provide useful tools for genetic mapping and phylogenetic studies of the guppy.

Opsin gene duplication and diversification in the guppy, a model for sexual selection.

Identification of genes that control variation in adaptive characters is a prerequisite for understanding the processes that drive sexual and natural selection. Male coloration and female colour perception play important roles in mate choice in the guppy (Poecilia reticulata), a model organism for studies of natural and sexual selection. We examined a potential source for the known variation in colour perception, by analysing genomic and complementary DNA sequences of genes that code for visual pigment proteins. We find high sequence variability, both within and between populations, and expanded copy number for long-wave sensitive (LWS) opsin genes. Alleles with non-synonymous changes that suggest dissimilar spectral tuning properties occur in the same population and even in the same individual, and the high frequency of non-synonymous substitutions argues for diversifying selection acting on these proteins. Therefore, variability in tuning amino acids is partitioned within individuals and populations of the guppy, in contrast to variability for LWS at higher taxonomic levels in cichlids, a second model system for differentiation owing to sexual selection. Since opsin variability parallels the extreme male colour polymorphism within guppy populations, we suggest that mate choice has been a major factor driving the coevolution of opsins and male ornaments in this species.

Protective Role of Complement C3 Against Cytokine-Mediated ß-Cell Apoptosis.

Type 1 diabetes is a chronic autoimmune disease characterized by pancreatic islet inflammation and ß-cell destruction by proinflammatory cytokines and other mediators. Based on RNA sequencing and protein-protein interaction analyses of human islets exposed to proinflammatory cytokines, we identified complement C3 as a hub for some of the effects of cytokines. The proinflammatory cytokines interleukin-1ß plus interferon-γ increase C3 expression in rodent and human pancreatic ß-cells, and C3 is detected by histology in and around the islets of diabetic patients. Surprisingly, C3 silencing exacerbates apoptosis under both basal condition and following exposure to cytokines, and it increases chemokine expression upon cytokine treatment. C3 exerts its prosurvival effects via AKT activation and c-Jun N-terminal kinase inhibition. Exogenously added C3 also protects against cytokine-induced ß-cell death and partially rescues the deleterious effects of inhibition of endogenous C3. These data suggest that locally produced C3 is an important prosurvival mechanism in pancreatic ß-cells under a proinflammatory assault.

Genome BLAST distance phylogenies inferred from whole plastid and whole mitochondrion genome sequences.

BACKGROUND: Phylogenetic methods which do not rely on multiple sequence alignments are important tools in inferring trees directly from completely sequenced genomes. Here, we extend the recently described Genome BLAST Distance Phylogeny (GBDP) strategy to compute phylogenetic trees from all completely sequenced plastid genomes currently available and from a selection of mitochondrial genomes representing the major eukaryotic lineages. BLASTN, TBLASTX, or combinations of both are used to locate high-scoring segment pairs (HSPs) between two sequences from which pairwise similarities and distances are computed in different ways resulting in a total of 96 GBDP variants. The suitability of these distance formulae for phylogeny reconstruction is directly estimated by computing a recently described measure of "treelikeness", the so-called delta value, from the respective distance matrices. Additionally, we compare the trees inferred from these matrices using UPGMA, NJ, BIONJ, FastME, or STC, respectively, with the NCBI taxonomy tree of the taxa under study. RESULTS: Our results indicate that, at this taxonomic level, plastid genomes are much more valuable for inferring phylogenies than are mitochondrial genomes, and that distances based on breakpoints are of little use. Distances based on the proportion of "matched" HSP length to average genome length were best for tree estimation. Additionally we found that using TBLASTX instead of BLASTN and, particularly, combining TBLASTX and BLASTN leads to a small but significant increase in accuracy. Other factors do not significantly affect the phylogenetic outcome. The BIONJ algorithm results in phylogenies most in accordance with the current NCBI taxonomy, with NJ and FastME performing insignificantly worse, and STC performing as well if applied to high quality distance matrices. delta values are found to be a reliable predictor of phylogenetic accuracy. CONCLUSION: Using the most treelike distance matrices, as judged by their delta values, distance methods are able to recover all major plant lineages, and are more in accordance with Apicomplexa organelles being derived from "green" plastids than from plastids of the "red" type. GBDP-like methods can be used to reliably infer phylogenies from different kinds of genomic data. A framework is established to further develop and improve such methods. delta values are a topology-independent tool of general use for the development and assessment of distance methods for phylogenetic inference.

Convergent targeting of a common host protein-network by pathogen effectors from three kingdoms of life.

While conceptual principles governing plant immunity are becoming clear, its systems-level organization and the evolutionary dynamic of the host-pathogen interface are still obscure. We generated a systematic protein-protein interaction network of virulence effectors from the ascomycete pathogen Golovinomyces orontii and Arabidopsis thaliana host proteins. We combined this data set with corresponding data for the eubacterial pathogen Pseudomonas syringae and the oomycete pathogen Hyaloperonospora arabidopsidis. The resulting network identifies host proteins onto which intraspecies and interspecies pathogen effectors converge. Phenotyping of 124 Arabidopsis effector-interactor mutants revealed a correlation between intraspecies and interspecies convergence and several altered immune response phenotypes. Several effectors and the most heavily targeted host protein colocalized in subnuclear foci. Products of adaptively selected Arabidopsis genes are enriched for interactions with effector targets. Our data suggest the existence of a molecular host-pathogen interface that is conserved across Arabidopsis accessions, while evolutionary adaptation occurs in the immediate network neighborhood of effector targets.

The Capsella rubella genome and the genomic consequences of rapid mating system evolution.

The shift from outcrossing to selfing is common in flowering plants, but the genomic consequences and the speed at which they emerge remain poorly understood. An excellent model for understanding the evolution of self fertilization is provided by Capsella rubella, which became self compatible <200,000 years ago. We report a C. rubella reference genome sequence and compare RNA expression and polymorphism patterns between C. rubella and its outcrossing progenitor Capsella grandiflora. We found a clear shift in the expression of genes associated with flowering phenotypes, similar to that seen in Arabidopsis, in which self fertilization evolved about 1 million years ago. Comparisons of the two Capsella species showed evidence of rapid genome-wide relaxation of purifying selection in C. rubella without a concomitant change in transposable element abundance. Overall we document that the transition to selfing may be typified by parallel shifts in gene expression, along with a measurable reduction of purifying selection.