Plasticity of maternal environment-dependent expression-QTLs of tomato seeds.

Seeds are essential for plant reproduction, survival, and dispersal. Germination ability and successful establishment of young seedlings strongly depend on seed quality and on environmental factors such as nutrient availability. In tomato (Solanum lycopersicum) and many other species, seed quality and seedling establishment characteristics are determined by genetic variation, as well as the maternal environment in which the seeds develop and mature. The genetic contribution to variation in seed and seedling quality traits and environmental responsiveness can be estimated at transcriptome level in the dry seed by mapping genomic loci that affect gene expression (expression QTLs) in contrasting maternal environments. In this study, we applied RNA-sequencing to construct a linkage map and measure gene expression of seeds of a tomato recombinant inbred line (RIL) population derived from a cross between S. lycopersicum (cv. Moneymaker) and S. pimpinellifolium (G1.1554). The seeds matured on plants cultivated under different nutritional environments, i.e., on high phosphorus or low nitrogen. The obtained single-nucleotide polymorphisms (SNPs) were subsequently used to construct a genetic map. We show how the genetic landscape of plasticity in gene regulation in dry seeds is affected by the maternal nutrient environment. The combined information on natural genetic variation mediating (variation in) responsiveness to the environment may contribute to knowledge-based breeding programs aiming to develop crop cultivars that are resilient to stressful environments.

Genetic mechanisms underlying increased microalgal thermotolerance, maximal growth rate, and yield on light following adaptive laboratory evolution.

BACKGROUND: Adaptive laboratory evolution (ALE) is a powerful method for strain optimization towards abiotic stress factors and for identifying adaptation mechanisms. In this study, the green microalga Picochlorum sp. BPE23 was cultured under supra-optimal temperature to force genetic adaptation. The robustness and adaptive capacity of Picochlorum strains turned them into an emerging model for evolutionary studies on abiotic stressors such as temperature, salinity, and light. RESULTS: Mutant strains showed an expanded maximal growth temperature of 44.6 °C, whereas the maximal growth temperature of the wild-type strain was 42 °C. Moreover, at the optimal growth temperature of 38 °C, the biomass yield on light was 22.3% higher, and the maximal growth rate was 70.5% higher than the wild type. Genome sequencing and transcriptome analysis were performed to elucidate the mechanisms behind the improved phenotype. A de novo assembled phased reference genome allowed the identification of 21 genic mutations involved in various processes. Moreover, approximately half of the genome contigs were found to be duplicated or even triplicated in all mutants, suggesting a causal role in adaptation. CONCLUSIONS: The developed tools and mutant strains provide a strong framework from whereupon Picochlorum sp. BPE23 can be further developed. Moreover, the extensive strain characterization provides evidence of how microalgae evolve to supra-optimal temperature and to photobioreactor growth conditions. With this study, microalgal evolutionary mechanisms were identified by combining ALE with genome sequencing.

Dormancy cycling: translation-related transcripts are the main difference between dormant and non-dormant seeds in the field.

Primary seed dormancy is a mechanism that orchestrates the timing of seed germination in order to prevent out-of-season germination. Secondary dormancy can be induced in imbibed seeds when they encounter prolonged unfavourable conditions. Secondary dormancy is not induced during dry storage, and therefore the mechanisms underlying this process have remained largely unexplored. Here, a 2-year seed burial experiment in which dormancy cycling was studied at the physiological and transcriptional level is presented. For these analyses six different Arabidopsis thaliana genotypes were used: Landsberg erecta (Ler) and the dormancy associated DELAY OF GERMINATION (DOG) near-isogenic lines 1, 2, 3, 6 and 22 (NILDOG1, 2, 3, 6 and 22). The germination potential of seeds exhumed from the field showed that these seeds go through dormancy cycling and that the dynamics of this cycling is genotype dependent. RNA-seq analysis revealed large transcriptional changes during dormancy cycling, especially at the time points preceding shifts in dormancy status. Dormancy cycling is driven by soil temperature and the endosperm is important in the perception of the environment. Genes that are upregulated in the low- to non-dormant stages are enriched for genes involved in translation, indicating that the non-dormant seeds are prepared for rapid seed germination.

Detection of QTLs for genotype × environment interactions in tomato seeds and seedlings.

Seed quality and seedling establishment are the most important factors affecting successful crop development. They depend on the genetic background and are acquired during seed maturation and therefor, affected by the maternal environment under which the seeds develop. There is little knowledge about the genetic and environmental factors that affect seed quality and seedling establishment. The aim of this study is to identify the loci and possible molecular mechanisms involved in acquisition of seed quality and how these are controlled by adverse maternal conditions. For this, we used a tomato recombinant inbred line (RIL) population consisting of 100 lines which were grown under two different nutritional environmental conditions, high phosphate and low nitrate. Most of the seed germination traits such as maximum germination percentage (Gmax ), germination rate (t50 ) and uniformity (U8416 ) showed ample variation between genotypes and under different germination conditions. This phenotypic variation leads to identification of quantitative trait loci (QTLs) which were dependent on genetic factors, but also on the interaction with the maternal environment (QTL × E). Further studies of these QTLs may ultimately help to predict the effect of different maternal environmental conditions on seed quality and seedling establishment which will be very useful to improve the production of high-performance seeds.

A multi-parent recombinant inbred line population of C. elegans allows identification of novel QTLs for complex life history traits.

BACKGROUND: The nematode Caenorhabditis elegans has been extensively used to explore the relationships between complex traits, genotypes, and environments. Complex traits can vary across different genotypes of a species, and the genetic regulators of trait variation can be mapped on the genome using quantitative trait locus (QTL) analysis of recombinant inbred lines (RILs) derived from genetically and phenotypically divergent parents. Most RILs have been derived from crossing two parents from globally distant locations. However, the genetic diversity between local C. elegans populations can be as diverse as between global populations and could thus provide means of identifying genetic variation associated with complex traits relevant on a broader scale. RESULTS: To investigate the effect of local genetic variation on heritable traits, we developed a new RIL population derived from 4 parental wild isolates collected from 2 closely located sites in France: Orsay and Santeuil. We crossed these 4 genetically diverse parental isolates to generate a population of 200 multi-parental RILs and used RNA-seq to obtain sequence polymorphisms identifying almost 9000 SNPs variable between the 4 genotypes with an average spacing of 11 kb, doubling the mapping resolution relative to currently available RIL panels for many loci. The SNPs were used to construct a genetic map to facilitate QTL analysis. We measured life history traits such as lifespan, stress resistance, developmental speed, and population growth in different environments, and found substantial variation for most traits. We detected multiple QTLs for most traits, including novel QTLs not found in previous QTL analysis, including those for lifespan and pathogen responses. This shows that recombining genetic variation across C. elegans populations that are in geographical close proximity provides ample variation for QTL mapping. CONCLUSION: Taken together, we show that using more parents than the classical two parental genotypes to construct a RIL population facilitates the detection of QTLs and that the use of wild isolates facilitates the detection of QTLs. The use of multi-parent RIL populations can further enhance our understanding of local adaptation and life history trade-offs.

Induction of desiccation tolerance in desiccation sensitive Citrus limon seeds.

Many economically important perennial species bear recalcitrant seeds, including tea, coffee, cocoa, mango, citrus, rubber, oil palm and coconut. Orthodox seeds can be dried almost completely without losing viability, but so-called recalcitrant seeds have a very limited storage life and die upon drying below a higher critical moisture content than orthodox seeds. As a result, the development of long-term storage methods for recalcitrant seeds is compromised. Lowering this critical moisture content would be very valuable since dry seed storage is the safest, most convenient and cheapest method for conserving plant genetic resources. Therefore, we have attempted to induce desiccation tolerance (DT) in the desiccation sensitive seeds of Citrus limon. We show that DT can be induced by paclobutrazol (an inhibitor of gibberellin biosynthesis) and we studied its associated transcriptome to delineate the molecular mechanisms underlying this induction of DT. Paclobutrazol not only interfered with gibberellin related gene expression but also caused extensive changes in expression of genes involved in the biosynthesis and signaling of other hormones. Paclobutrazol induced a transcriptomic switch encompassing suppression of biotic- and induction of abiotic responses. We hypothesize that this is the main driver of the induction of DT by paclobutrazol in C. limon seeds.

AraQTL - workbench and archive for systems genetics in Arabidopsis thaliana.

Genetical genomics studies uncover genome-wide genetic interactions between genes and their transcriptional regulators. High-throughput measurement of gene expression in recombinant inbred line populations has enabled investigation of the genetic architecture of variation in gene expression. This has the potential to enrich our understanding of the molecular mechanisms affected by and underlying natural variation. Moreover, it contributes to the systems biology of natural variation, as a substantial number of experiments have resulted in a valuable amount of interconnectable phenotypic, molecular and genotypic data. A number of genetical genomics studies have been published for Arabidopsis thaliana, uncovering many expression quantitative trait loci (eQTLs). However, these complex data are not easily accessible to the plant research community, leaving most of the valuable genetic interactions unexplored as cross-analysis of these studies is a major effort. We address this problem with AraQTL (http://www.bioinformatics.nl/Ara QTL/), an easily accessible workbench and database for comparative analysis and meta-analysis of all published Arabidopsis eQTL datasets. AraQTL provides a workbench for comparing, re-using and extending upon the results of these experiments. For example, one can easily screen a physical region for specific local eQTLs that could harbour candidate genes for phenotypic QTLs, or detect gene-by-environment interactions by comparing eQTLs under different conditions.

Tulipa gesneriana and Lilium longiflorum PEBP Genes and Their Putative Roles in Flowering Time Control.

Floral induction in Tulipa gesneriana and Lilium longiflorum is triggered by contrasting temperature conditions, high and low temperature, respectively. In Arabidopsis, the floral integrator FLOWERING LOCUS T (FT), a member of the PEBP (phosphatidyl ethanolamine-binding protein) gene family, is a key player in flowering time control. In this study, one PEBP gene was identified and characterized in lily (LlFT) and three PEBP genes were isolated from tulip (TgFT1, TgFT2 and TgFT3). Overexpression of these genes in Arabidopsis thaliana resulted in an early flowering phenotype for LlFT and TgFT2, but a late flowering phenotype for TgFT1 and TgFT3. Overexpression of LlFT in L. longiflorum also resulted in an early flowering phenotype, confirming its proposed role as a flowering time-controlling gene. The tulip PEBP genes TgFT2 and TgFT3 have a similar expression pattern in tulip, but show opposite effects on the timing of flowering in Arabidopsis. Therefore, the difference between these two proteins was further investigated by interchanging amino acids thought to be important for the FT function. This resulted in the conversion of phenotypes in Arabidopsis upon overexpressing the substituted TgFT2 and TgFT3 genes, revealing the importance of these interchanged amino acid residues. Based on all obtained results, we hypothesize that LlFT is involved in creating meristem competence to flowering-related cues in lily, and TgFT2 is considered to act as a florigen involved in the floral induction in tulip. The function of TgFT3 remains unclear, but, based on our observations and phylogenetic analysis, we propose a bulb-specific function for this gene.

Molecular Regulation of Temperature-Dependent Floral Induction in Tulipa gesneriana.

The vegetative-to-reproductive phase change in tulip (Tulipa gesneriana) is promoted by increasing temperatures during spring. The warm winters of recent years interfere with this process and are calling for new adapted cultivars. A better understanding of the underlying molecular mechanisms would be of help, but unlike the model plant Arabidopsis (Arabidopsis thaliana), very little is known about floral induction in tulip. To shed light on the gene regulatory network controlling flowering in tulip, RNA sequencing was performed on meristem-enriched tissue collected under two contrasting temperature conditions, low and high. The start of reproductive development correlated with rounding of the shoot apical meristem and induction of TGSQA expression, a tulip gene with a high similarity to Arabidopsis APETALA1 Gene Ontology enrichment analysis of differentially expressed genes showed the overrepresentation of genes potentially involved in floral induction, bulb maturation, and dormancy establishment. Expression analysis revealed that TERMINAL FLOWER1 (TgTFL1) and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1-like1 (TgSOC1-like1) might be repressors, whereas TgSOC1-like2 likely is an activator, of flowering. Subsequently, the flowering time-associated expression of eight potential flowering time genes was confirmed in three tulip cultivars grown in the field. Additionally, heterologous functional analyses in Arabidopsis resulted in flowering time phenotypes in line with TgTFL1 being a floral repressor and TgSOC1-like2 being a floral activator in tulip. Taken together, we have shown that long before morphological changes occur in the shoot apical meristem, the expression of floral repressors in tulip is suppressed by increased ambient temperatures, leading either directly or indirectly to the activation of potential flowering activators shortly before the commencement of the phase change.

Elucidating and mining the Tulipa and Lilium transcriptomes.

Genome sequencing remains a challenge for species with large and complex genomes containing extensive repetitive sequences, of which the bulbous and monocotyledonous plants tulip and lily are examples. In such a case, sequencing of only the active part of the genome, represented by the transcriptome, is a good alternative to obtain information about gene content. In this study we aimed to generate a high quality transcriptome of tulip and lily and to make this data available as an open-access resource via a user-friendly web-based interface. The Illumina HiSeq 2000 platform was applied and the transcribed RNA was sequenced from a collection of different lily and tulip tissues, respectively. In order to obtain good transcriptome coverage and to facilitate effective data mining, assembly was done using different filtering parameters for clearing out contamination and noise of the RNAseq datasets. This analysis revealed limitations of commonly applied methods and parameter settings used in de novo transcriptome assembly. The final created transcriptomes are publicly available via a user friendly Transcriptome browser ( http://www.bioinformatics.nl/bulbs/db/species/index ). The usefulness of this resource has been exemplified by a search for all potential transcription factors in lily and tulip, with special focus on the TCP transcription factor family. This analysis and other quality parameters point out the quality of the transcriptomes, which can serve as a basis for further genomics studies in lily, tulip, and bulbous plants in general.

Effects of Parental Temperature and Nitrate on Seed Performance are Reflected by Partly Overlapping Genetic and Metabolic Pathways.

Seed performance is affected by the seed maturation environment, and previously we have shown that temperature, nitrate and light intensity were the most influential environmental factors affecting seed performance. Seeds developed in these environments were selected to assess the underlying metabolic pathways, using a combination of transcriptomics and metabolomics. These analyses revealed that the effects of the parental temperature and nitrate environments were reflected by partly overlapping genetic and metabolic networks, as indicated by similar changes in the expression levels of metabolites and transcripts. Nitrogen metabolism-related metabolites (asparagine, γ-aminobutyric acid and allantoin) were significantly decreased in both low temperature (15 °C) and low nitrate (N0) maturation environments. Correspondingly, nitrogen metabolism genes (ALLANTOINASE, NITRATE REDUCTASE 1, NITRITE REDUCTASE 1 and NITRILASE 4) were differentially regulated in the low temperature and nitrate maturation environments, as compared with control conditions. High light intensity during seed maturation increased galactinol content, and displayed a high correlation with seed longevity. Low light had a genotype-specific effect on cell surface-encoding genes in the DELAY OF GERMINATION 6-near isogenic line (NILDOG6). Overall, the integration of phenotypes, metabolites and transcripts led to new insights into the regulation of seed performance.

Understanding and identifying amino acid repeats.

Amino acid repeats (AARs) are abundant in protein sequences. They have particular roles in protein function and evolution. Simple repeat patterns generated by DNA slippage tend to introduce length variations and point mutations in repeat regions. Loss of normal and gain of abnormal function owing to their variable length are potential risks leading to diseases. Repeats with complex patterns mostly refer to the functional domain repeats, such as the well-known leucine-rich repeat and WD repeat, which are frequently involved in protein­protein interaction. They are mainly derived from internal gene duplication events and stabilized by 'gate-keeper' residues, which play crucial roles in preventing inter-domain aggregation. AARs are widely distributed in different proteomes across a variety of taxonomic ranges, and especially abundant in eukaryotic proteins. However, their specific evolutionary and functional scenarios are still poorly understood. Identifying AARs in protein sequences is the first step for the further investigation of their biological function and evolutionary mechanism. In principle, this is an NP-hard problem, as most of the repeat fragments are shaped by a series of sophisticated evolutionary events and become latent periodical patterns. It is not possible to define a uniform criterion for detecting and verifying various repeat patterns. Instead, different algorithms based on different strategies have been developed to cope with different repeat patterns. In this review, we attempt to describe the amino acid repeat-detection algorithms currently available and compare their strategies based on an in-depth analysis of the biological significance of protein repeats.

A gene co-expression network predicts functional genes controlling the re-establishment of desiccation tolerance in germinated Arabidopsis thaliana seeds.

MAIN CONCLUSION: During re-establishment of desiccation tolerance (DT), early events promote initial protection and growth arrest, while late events promote stress adaptation and contribute to survival in the dry state. Mature seeds of Arabidopsis thaliana are desiccation tolerant, but they lose desiccation tolerance (DT) while progressing to germination. Yet, there is a small developmental window during which DT can be rescued by treatment with abscisic acid (ABA). To gain temporal resolution and identify relevant genes in this process, data from a time series of microarrays were used to build a gene co-expression network. The network has two regions, namely early response (ER) and late response (LR). Genes in the ER region are related to biological processes, such as dormancy, acquisition of DT and drought, amplification of signals, growth arrest and induction of protection mechanisms (such as LEA proteins). Genes in the LR region lead to inhibition of photosynthesis and primary metabolism, promote adaptation to stress conditions and contribute to seed longevity. Phenotyping of 12 hubs in relation to re-establishment of DT with T-DNA insertion lines indicated a significant increase in the ability to re-establish DT compared with the wild-type in the lines cbsx4, at3g53040 and at4g25580, suggesting the operation of redundant and compensatory mechanisms. Moreover, we show that re-establishment of DT by polyethylene glycol and ABA occurs through partially overlapping mechanisms. Our data confirm that co-expression network analysis is a valid approach to examine data from time series of transcriptome analysis, as it provides promising insights into biologically relevant relations that help to generate new information about the roles of certain genes for DT.

Discovery of T cell epitopes implementing HLA-peptidomics into a reverse immunology approach.

T cell recognition of minor histocompatibility Ags (MiHA) plays an important role in the graft-versus-tumor effect of allogeneic stem cell transplantation. Selective infusion of T cells reactive for hematopoiesis-restricted MiHA presented in the context of HLA class I or II molecules may help to separate the graft-versus-tumor effects from graft-versus-host disease effects after allogeneic stem cell transplantation. Over the years, increasing numbers of MiHA have been identified by forward immunology approaches, and the relevance of these MiHA has been illustrated by correlation with clinical outcome. As the tissue distribution of MiHA affects the clinical outcome of T cell responses against these Ags, it would be beneficial to identify additional predefined MiHA that are exclusively expressed on hematopoietic cells. Therefore, several reverse immunology approaches have been explored for the prediction of MiHA. Thus far, these approaches frequently resulted in the identification of T cells directed against epitopes that are not naturally processed and presented. In this study we established a method for the identification of biologically relevant MiHA, implementing mass spectrometry-based HLA-peptidomics into a reverse immunology approach. For this purpose, HLA class I binding peptides were eluted from transformed B cells, analyzed by mass spectrometry, and matched with a database dedicated to identifying polymorphic peptides. This process resulted in a set of 40 MiHA candidates that were evaluated in multiple selection steps. The identification of LB-NISCH-1A demonstrated the technical feasibility of our approach. On the basis of these results, we present an approach that can be of value for the efficient identification of MiHA or other T cell epitopes.

The human leukocyte antigen-presented ligandome of B lymphocytes.

Peptides presented by human leukocyte antigen (HLA) molecules on the cell surface play a crucial role in adaptive immunology, mediating the communication between T cells and antigen presenting cells. Knowledge of these peptides is of pivotal importance in fundamental studies of T cell action and in cellular immunotherapy and transplantation. In this paper we present the in-depth identification and relative quantification of 14,500 peptide ligands constituting the HLA ligandome of B cells. This large number of identified ligands provides general insight into the presented peptide repertoire and antigen presentation. Our uniquely large set of HLA ligands allowed us to characterize in detail the peptides constituting the ligandome in terms of relative abundance, peptide length distribution, physicochemical properties, binding affinity to the HLA molecule, and presence of post-translational modifications. The presented B-lymphocyte ligandome is shown to be a rich source of information by the presence of minor histocompatibility antigens, virus-derived epitopes, and post-translationally modified HLA ligands, and it can be a good starting point for solving a wealth of specific immunological questions. These HLA ligands can form the basis for reversed immunology approaches to identify T cell epitopes based not on in silico predictions but on the bona fide eluted HLA ligandome.

QualitySNPng: a user-friendly SNP detection and visualization tool.

QualitySNPng is a new software tool for the detection and interactive visualization of single-nucleotide polymorphisms (SNPs). It uses a haplotype-based strategy to identify reliable SNPs; it is optimized for the analysis of current RNA-seq data; but it can also be used on genomic DNA sequences derived from next-generation sequencing experiments. QualitySNPng does not require a sequenced reference genome and delivers reliable SNPs for di- as well as polyploid species. The tool features a user-friendly interface, multiple filtering options to handle typical sequencing errors, support for SAM and ACE files and interactive visualization. QualitySNPng produces high-quality SNP information that can be used directly in genotyping by sequencing approaches for application in QTL and genome-wide association mapping as well as to populate SNP arrays. The software can be used as a stand-alone application with a graphical user interface or as part of a pipeline system like Galaxy. Versions for Windows, Mac OS X and Linux, as well as the source code, are available from http://www.bioinformatics.nl/QualitySNPng.

Beyond genomic variation--comparison and functional annotation of three Brassica rapa genomes: a turnip, a rapid cycling and a Chinese cabbage.

BACKGROUND: Brassica rapa is an economically important crop species. During its long breeding history, a large number of morphotypes have been generated, including leafy vegetables such as Chinese cabbage and pakchoi, turnip tuber crops and oil crops. RESULTS: To investigate the genetic variation underlying this morphological variation, we re-sequenced, assembled and annotated the genomes of two B. rapa subspecies, turnip crops (turnip) and a rapid cycling. We then analysed the two resulting genomes together with the Chinese cabbage Chiifu reference genome to obtain an impression of the B. rapa pan-genome. The number of genes with protein-coding changes between the three genotypes was lower than that among different accessions of Arabidopsis thaliana, which can be explained by the smaller effective population size of B. rapa due to its domestication. Based on orthology to a number of non-brassica species, we estimated the date of divergence among the three B. rapa morphotypes at approximately 250,000 YA, far predating Brassica domestication (5,000-10,000 YA). CONCLUSIONS: By analysing genes unique to turnip we found evidence for copy number differences in peroxidases, pointing to a role for the phenylpropanoid biosynthesis pathway in the generation of morphological variation. The estimated date of divergence among three B. rapa morphotypes implies that prior to domestication there was already considerably divergence among B. rapa genotypes. Our study thus provides two new B. rapa reference genomes, delivers a set of computer tools to analyse the resulting pan-genome and uses these to shed light on genetic drivers behind the rich morphological variation found in B. rapa.

Interaction between parental environment and genotype affects plant and seed performance in Arabidopsis.

Seed performance after dispersal is highly dependent on parental environmental cues, especially during seed formation and maturation. Here we examine which environmental factors are the most dominant in this respect and whether their effects are dependent on the genotypes under investigation. We studied the influence of light intensity, photoperiod, temperature, nitrate, and phosphate during seed development on five plant attributes and thirteen seed attributes, using 12 Arabidopsis genotypes that have been reported to be affected in seed traits. As expected, the various environments during seed development resulted in changed plant and/or seed performances. Comparative analysis clearly indicated that, overall, temperature plays the most dominant role in both plant and seed performance, whereas light has a prominent impact on plant traits. In comparison to temperature and light, nitrate mildly affected some of the plant and seed traits while phosphate had even less influence on those traits. Moreover, clear genotype-by-environment interactions were identified. This was shown by the fact that individual genotypes responded differentially to the environmental conditions. Low temperature significantly increased seed dormancy and decreased seed longevity of NILDOG1 and cyp707a1-1, whereas low light intensity increased seed dormancy and decreased seed longevity of NILDOG3 and NILDOG6. This also indicates that different genetic and molecular pathways are involved in the plant and seed responses. By identifying environmental conditions that affect the dormancy vs longevity correlation in the same way as previously identified naturally occurring loci, we have identified selective forces that probably shaped evolution for these important seed traits.

Promoter propagation in prokaryotes.

Transcriptional activation or 'rewiring' of silent genes is an important, yet poorly understood, phenomenon in prokaryotic genomes. Anecdotal evidence coming from experimental evolution studies in bacterial systems has shown the promptness of adaptation upon appropriate selective pressure. In many cases, a partial or complete promoter is mobilized to silent genes from elsewhere in the genome. We term hereafter such recruited regulatory sequences as Putative Mobile Promoters (PMPs) and we hypothesize they have a large impact on rapid adaptation of novel or cryptic functions. Querying all publicly available prokaryotic genomes (1362) uncovered >4000 families of highly conserved PMPs (50 to 100 long with ≥80% nt identity) in 1043 genomes from 424 different genera. The genomes with the largest number of PMP families are Anabaena variabilis (28 families), Geobacter uraniireducens (27 families) and Cyanothece PCC7424 (25 families). Family size varied from 2 to 93 homologous promoters (in Desulfurivibrio alkaliphilus). Some PMPs are present in particular species, but some are conserved across distant genera. The identified PMPs represent a conservative dataset of very recent or conserved events of mobilization of non-coding DNA and thus they constitute evidence of an extensive reservoir of recyclable regulatory sequences for rapid transcriptional rewiring.

ProRepeat: an integrated repository for studying amino acid tandem repeats in proteins.

ProRepeat (http://prorepeat.bioinformatics.nl/) is an integrated curated repository and analysis platform for in-depth research on the biological characteristics of amino acid tandem repeats. ProRepeat collects repeats from all proteins included in the UniProt knowledgebase, together with 85 completely sequenced eukaryotic proteomes contained within the RefSeq collection. It contains non-redundant perfect tandem repeats, approximate tandem repeats and simple, low-complexity sequences, covering the majority of the amino acid tandem repeat patterns found in proteins. The ProRepeat web interface allows querying the repeat database using repeat characteristics like repeat unit and length, number of repetitions of the repeat unit and position of the repeat in the protein. Users can also search for repeats by the characteristics of repeat containing proteins, such as entry ID, protein description, sequence length, gene name and taxon. ProRepeat offers powerful analysis tools for finding biological interesting properties of repeats, such as the strong position bias of leucine repeats in the N-terminus of eukaryotic protein sequences, the differences of repeat abundance among proteomes, the functional classification of repeat containing proteins and GC content constrains of repeats' corresponding codons.