The TIGR gene indices: reconstruction and representation of expressed gene sequences

The bs5 allele of the susceptibility gene Bs5 of pepper (Capsicum annuum L.) encoding a natural deletion variant of a CYSTM protein conditions resistance to bacterial spot disease caused by Xanthomonas species

KEY MESSAGE

The bs5 resistance gene against bacterial spot was identified by map-based cloning. The recessive bs5 gene of pepper (Capsicum annuum L.) conditions a non-hypersensitive resistance trait, characterized by a slightly swollen, pale green, photosynthetically active leaf tissue, following Xanthomonas euvesicatoria infection. The isolation of the bs5 gene by map-based cloning revealed that the bs5 protein was shorter by 2 amino acids as compared to the wild type Bs5 protein. The natural 2 amino acid deletion occurred in the cysteine-rich transmembrane domain of the tail-anchored (TA) protein, Ca_CYSTM1. The protein products of the wild type Bs5 and mutant bs5 genes were shown to be located in the cell membrane, indicating an unknown function in this membrane compartment. Successful infection of the Bs5 pepper lines was abolished by the 6 bp deletion in the TM encoding domain of the Ca_CYSTM1 gene in bs5 homozygotes, suggesting, that the resulting resistance might be explained by the lack of entry of the Xanthomonas specific effector molecules into the plant cells.

A transcriptomic view to wounding response in young Scots pine stems

We studied the stress response of five-year-old Scots pine xylem to mechanical wounding using RNA sequencing. In general, we observed a bimodal response in pine xylem after wounding. Transcripts associated with water deficit stress, defence, and cell wall modification were induced at the earliest time point of three hours; at the same time, growth-related processes were down-regulated. A second temporal wave was triggered either at the middle and/or at the late time points (one and four days). Secondary metabolism, such as stilbene and lignan biosynthesis started one day after wounding. Scots pine synthesises the stilbenes pinosylvin and its monomethyl ether both as constitutive and induced defence compounds. Stilbene biosynthesis is induced by wounding, pathogens and UV stress, but is also developmentally regulated when heartwood is formed. Comparison of wounding responses to heartwood formation shows that many induced processes (in addition to stilbene biosynthesis) are similar and relate to defence or desiccation stress, but often specific transcripts are up-regulated in the developmental and wounding induced contexts. Pine resin biosynthesis was not induced in response to wounding, at least not during the first four days.

De novo clustering of long reads by gene from transcriptomics data

Long-read sequencing currently provides sequences of several thousand base pairs. It is therefore possible to obtain complete transcripts, offering an unprecedented vision of the cellular transcriptome. However the literature lacks tools for de novo clustering of such data, in particular for Oxford Nanopore Technologies reads, because of the inherent high error rate compared to short reads. Our goal is to process reads from whole transcriptome sequencing data accurately and without a reference genome in order to reliably group reads coming from the same gene. This de novo approach is therefore particularly suitable for non-model species, but can also serve as a useful pre-processing step to improve read mapping. Our contribution both proposes a new algorithm adapted to clustering of reads by gene and a practical and free access tool that allows to scale the complete processing of eukaryotic transcriptomes. We sequenced a mouse RNA sample using the MinION device. This dataset is used to compare our solution to other algorithms used in the context of biological clustering. We demonstrate that it is the best approach for transcriptomics long reads. When a reference is available to enable mapping, we show that it stands as an alternative method that predicts complementary clusters.

Predictions of novel Schistosoma mansoni - human protein interactions consistent with experimental data

Infection by the human blood fluke, Schistosoma mansoni involves a variety of cross-species protein- protein interactions. The pathogen expresses a diverse arsenal of proteins that facilitate the breach of physical and biochemical barriers present in skin evasion of the immune system, and digestion of human plasma proteins including albumin and hemoglobin, allowing schistosomes to reside in the host for years. However, only a small number of specific interactions between S. mansoni and human proteins have been identified. We present and apply a protocol that generates testable predictions of S. mansoni-human protein interactions. In this study, we have preliminary predictions of novel interactions between schistosome and human proteins relevant to infection and the ability of the parasite to evade the immune system. We applied a computational whole-genome comparative approach to predict potential S. mansoni-human protein interactions based on similarity to known protein complexes. We first predict S. mansoni -human protein interactions based on similarity to known protein complexes. Putative interactions were then scored and assessed using several contextual filters, including the use of annotation automatically derived from literature using a simple natural language processing methodology. Next, in vitro experiments were carried out between schistosome and host proteins to validate several prospective predictions. Our method predicted 7 out of the 10 previously known cross-species interactions involved in pathogenesis between S. mansoni and its human host. Interestingly, two novel putative interactions involving Schistosoma proteins, the cercarial elastase SmCE, and the adult tegument surface protein Sm29, were also predicted and experimentally characterized. Preliminary data suggest that elafin, a host endogenous serine protease inhibitor, may be a novel substrate for SmCE. Additionally, CD59, an inhibitor of the membrane attack complex, could interact with Sm29. Furthermore, the application framework provides an integrated methodology for investigation of host-pathogen interactions and an extensive source of orthogonal data for experimental analysis. We have made the predictions available for community perusal.

The O-methyltransferase PMT2 mediates methylation of pinosylvin in Scots pine

Heartwood extractives are important determinants of the natural durability of pine heartwood. The most important phenolic compounds affecting durability are the stilbenes pinosylvin and its monomethylether, which in addition have important functions as phytoalexins in active defense. A substantial portion of the synthesized pinosylvin is 3-methoxylated but the O-methyltransferase responsible for this modification has not been correctly identified. We studied the expression of the stilbene pathway during heartwood development as well as in response to wounding of xylem and UV-C treatment of needles. We isolated and enzymatically characterized a novel O-methyltransferase, PMT2. The methylated product was verified as pinosylvin monomethylether using ultra performance liquid chromatography-tandem mass spectrometry and high performance liquid chromatography analyses. The PMT2 enzyme was highly specific for stilbenes as substrate, in contrast to caffeoyl-CoA O-methyltransferase (CCoAOMT) and PMT1 that were multifunctional. Expression profile and multifunctional activity of CCoAOMT suggest that it might have additional roles outside lignin biosynthesis. PMT1 is not involved in the stilbene pathway and its biological function remains an open question. We isolated a new specific O-methyltransferase responsible for 3-methoxylation of pinosylvin. Expression of PMT2 closely follows stilbene biosynthesis during developmental and stress induction. We propose that PMT2 is responsible for pinosylvin methylation in Scots pine (Pinus sylvestris), instead of the previously characterized methyltransferase, PMT1.

Aberration of the Citation

Multiple inherent biases related to different citation practices (for e.g., self-citations, negative citations, wrong citations, multi-authorship-biased citations, honorary citations, circumstantial citations, discriminatory citations, selective and arbitrary citations, etc.) make citation-based bibliometrics strongly flawed and defective measures. A paper can be highly cited for a while (for e.g., under circumstantial or transitional knowledge), but years later it may appear that its findings, paradigms, or theories were untrue or invalid anymore. By contrast, a paper may remain shelved or overlooked for years or decades, but new studies or discoveries may actualize its subject at any moment. As citation-based metrics are transformed into "commercial activities," the "citation credit" should be considered on a commercial basis too, in the sense that "citation credit" should be shared out as a "citation dividend" by shareholders (coauthors) averagely or proportionally to their contributions but not fully appropriated by each of them. At equal numbers of citations, the greater number of authors, the lower "citation credit" should be and vice versa. Overlooking the presence of distorted and subjective citation practices makes many people and administrators "obsessed" with the number of citations to such an extent to run after "highly cited" authors and to create specialized citation databases for commercial purposes. Citation-based bibliometrics, however, are unreliable and unscientific measures; citation counts do not mean that a more cited work is of a higher quality or accuracy than a less cited work because citations do not measure the quality or accuracy. Citations do not mean that a highly cited author or journal is more commendable than a less cited author or journal. Citations are not more than countable numbers: no more, no less.

Abscisic-acid-dependent basic leucine zipper (bZIP) transcription factors in plant abiotic stress

One of the major causes of significant crop loss throughout the world is the myriad of environmental stresses including drought, salinity, cold, heavy metal toxicity, and ultraviolet-B (UV-B) rays. Plants as sessile organisms have evolved various effective mechanism which enable them to withstand this plethora of stresses. Most of such regulatory mechanisms usually follow the abscisic-acid (ABA)-dependent pathway. In this review, we have primarily focussed on the basic leucine zipper (bZIP) transcription factors (TFs) activated by the ABA-mediated signalosome. Upon perception of ABA by specialized receptors, the signal is transduced via various groups of Ser/Thr kinases, which phosphorylate the bZIP TFs. Following such post-translational modification of TFs, they are activated so that they bind to specific cis-acting sequences called abscisic-acid-responsive elements (ABREs) or GC-rich coupling elements (CE), thereby influencing the expression of their target downstream genes. Several in silico techniques have been adopted so far to predict the structural features, recognize the regulatory modification sites, undergo phylogenetic analyses, and facilitate genome-wide survey of TF under multiple stresses. Current investigations on the epigenetic regulation that controls greater accessibility of the inducible regions of DNA of the target gene to the bZIP TFs exclusively under stress situations, along with the evolved stress memory responses via genomic imprinting mechanism, have been highlighted. The potentiality of overexpression of bZIP TFs, either in a homologous or in a heterologous background, in generating transgenic plants tolerant to various abiotic stressors have also been addressed by various groups. The present review will provide a coherent documentation on the functional characterization and regulation of bZIP TFs under multiple environmental stresses, with the major goal of generating multiple-stress-tolerant plant cultivars in near future.

Developmental Changes in Scots Pine Transcriptome during Heartwood Formation

Scots pine (Pinus sylvestris L.) wood is desired in woodworking industries due to its favorable timber characteristics and natural durability that is contributed by heartwood extractives. It has been discussed whether the Scots pine heartwood extractives (mainly stilbenes and resin acids) are synthesized in the cells of the transition zone between sapwood and heartwood, or if they are transported from the sapwood. Timing of heartwood formation during the yearly cycle has also not been unambiguously defined. We measured steady-state mRNA levels in Scots pine transition zone and sapwood using RNA sequencing. Year-round expression profiles of selected transcripts were further investigated by quantitative RT-PCR. Differentially accumulating transcripts suggest that, of the Scots pine heartwood extractives, stilbenes are synthesized in situ in the transition zone and gain their carbon-skeletons from Suc and triglycerides. Resin acids, on the other hand, are synthesized early in the spring mainly in the sapwood, meaning that they must be transported to the heartwood transition zone. Heartwood formation is marked by programmed cell death that occurs during the summer months in the transition zone.

Quantitative Expression Analysis in Brassica napus by Northern Blot Analysis and Reverse Transcription-Quantitative PCR in a Complex Experimental Setting

Analysis of gene expression is one of the major ways to better understand plant reactions to changes in environmental conditions. The comparison of many different factors influencing plant growth challenges the gene expression analysis for specific gene-targeted experiments, especially with regard to the choice of suitable reference genes. The aim of this study is to compare expression results obtained by Northern blot, semi-quantitative PCR and RT-qPCR, and to identify a reliable set of reference genes for oilseed rape (Brassica napus L.) suitable for comparing gene expression under complex experimental conditions. We investigated the influence of several factors such as sulfur deficiency, different time points during the day, varying light conditions, and their interaction on gene expression in oilseed rape plants. The expression of selected reference genes was indeed influenced under these conditions in different ways. Therefore, a recently developed algorithm, called GrayNorm, was applied to validate a set of reference genes for normalizing results obtained by Northern blot analysis. After careful comparison of the three methods mentioned above, Northern blot analysis seems to be a reliable and cost-effective alternative for gene expression analysis under a complex growth regime. For using this method in a quantitative way a number of references was validated revealing that for our experiment a set of three references provides an appropriate normalization. Semi-quantitative PCR was prone to many handling errors and difficult to control while RT-qPCR was very sensitive to expression fluctuations of the reference genes.

Heterotrimeric G-proteins in Picea abies and their regulation in response to Heterobasidion annosum s.l. infection

BACKGROUND

Heterotrimeric G-proteins are important signalling switches, present in all eukaryotic kingdoms. In plants they regulate several developmental functions and play an important role in plant-microbe interactions. The current knowledge on plant G-proteins is mostly based on model angiosperms and little is known about the G-protein repertoire and function in other lineages. In this study we investigate the heterotrimeric G-protein subunit repertoire in Pinaceae, including phylogenetic relationships, radiation and sequence diversity levels in relation to other plant linages. We also investigate functional diversification of the G-protein complex in Picea abies by analysing transcriptional regulation of the G-protein subunits in different tissues and in response to pathogen infection.

RESULTS

A full repertoire of G-protein subunits in several conifer species were identified in silico. The full-length P. abies coding regions of one Gα-, one Gβ- and four Gγ-subunits were cloned and sequenced. The phylogenetic analysis of the Gγ-subunits showed that PaGG1 clustered with A-type-like subunits, PaGG3 and PaGG4 clustered with C-type-like subunits, while PaGG2 and its orthologs represented a novel conifer-specific putative Gγ-subunit type. Gene expression analyses by quantitative PCR of P. abies G-protein subunits showed specific up-regulation of the Gα-subunit gene PaGPA1 and the Gγ-subunit gene PaGG1 in response to Heterobasidion annosum sensu lato infection.

CONCLUSIONS

Conifers possess a full repertoire of G-protein subunits. The differential regulation of PaGPA1 and PaGG1 indicates that the heterotrimeric G-protein complex represents a critical linchpin in Heterobasidion annosum s.l. perception and downstream signaling in P. abies.

Molecular marker development from transcript sequences and germplasm evaluation for cultivated peanut (Arachis hypogaea L.)

Molecular markers are important tools for genotyping in genetic studies and molecular breeding. The SSR and SNP are two commonly used marker systems developed from genomic or transcript sequences. The objectives of this study were to: (1) assemble and annotate the publicly available ESTs in Arachis and the in-house short reads, (2) develop and validate SSR and SNP markers, and (3) investigate the genetic diversity and population structure of the peanut breeding lines and the U.S. peanut mini core collection using developed SSR markers. An NCBI EST dataset with 252,951 sequences and an in-house 454 RNAseq dataset with 288,701 sequences were assembled separately after trimming. Transcript sequence comparison and phylogenetic analysis suggested that peanut is closer to cowpea and scarlet bean than to soybean, common bean and Medicago. From these two datasets, 6455 novel SSRs and 11,902 SNPs were identified. Of the discovered SSRs, 380 representing various SSR types were selected for PCR validation. The amplification rate was 89.2 %. Twenty-two (6.5 %) SSRs were polymorphic between at least one pair of four genotypes. Sanger sequencing of PCR products targeting 110 SNPs revealed 13 true SNPs between tetraploid genotypes and 193 homoeologous SNPs within genotypes. Eight out of the 22 polymorphic SSR markers were selected to evaluate the genetic diversity of Florida peanut breeding lines and the U.S. peanut mini core collection. This marker set demonstrated high discrimination power by displaying an average polymorphism information content value of 0.783, a combined probability of identity of 10(-11), and a combined power of exclusion of 0.99991. The structure analysis revealed four sub-populations among the peanut accessions and lines evaluated. The results of this study enriched the peanut genomic resources, provided over 6000 novel SSR markers and the credentials for true peanut SNP marker development, and demonstrated the power of newly developed SSR markers in genotyping peanut germplasm and breeding materials.

Whole transcriptome analysis with sequencing: methods, challenges and potential solutions

Whole transcriptome analysis plays an essential role in deciphering genome structure and function, identifying genetic networks underlying cellular, physiological, biochemical and biological systems and establishing molecular biomarkers that respond to diseases, pathogens and environmental challenges. Here, we review transcriptome analysis methods and technologies that have been used to conduct whole transcriptome shotgun sequencing or whole transcriptome tag/target sequencing analyses. We focus on how adaptors/linkers are added to both 5' and 3' ends of mRNA molecules for cloning or PCR amplification before sequencing. Challenges and potential solutions are also discussed. In brief, next generation sequencing platforms have accelerated releases of the large amounts of gene expression data. It is now time for the genome research community to assemble whole transcriptomes of all species and collect signature targets for each gene/transcript, and thus use known genes/transcripts to determine known transcriptomes directly in the near future.

Extreme Sensory Complexity Encoded in the 10-Megabase Draft Genome Sequence of the Chromatically Acclimating Cyanobacterium Tolypothrix sp. PCC 7601

Tolypothrix sp. PCC 7601 is a freshwater filamentous cyanobacterium with complex responses to environmental conditions. Here, we present its 9.96-Mbp draft genome sequence, containing 10,065 putative protein-coding sequences, including 305 predicted two-component system proteins and 27 putative phytochrome-class photoreceptors, the most such proteins in any sequenced genome.

Brassica napus L. cultivars show a broad variability in their morphology, physiology and metabolite levels in response to sulfur limitations and to pathogen attack

Under adequate sulfur supply, plants accumulate sulfate in the vacuoles and use sulfur-containing metabolites as storage compounds. Under sulfur-limiting conditions, these pools of stored sulfur-compounds are depleted in order to balance the nitrogen to sulfur ratio for protein synthesis. Stress conditions like sulfur limitation and/or pathogen attack induce changes in the sulfate pool and the levels of sulfur-containing metabolites, which often depend on the ecotypes or cultivars. We are interested in investigating the influence of the genetic background of canola (Brassica napus) cultivars in sulfur-limiting conditions on the resistance against Verticillium longisporum. Therefore, four commercially available B. napus cultivars were analyzed. These high-performing cultivars differ in some characteristics described in their cultivar pass, such as several agronomic traits, differences in the size of the root system, and resistance to certain pathogens, such as Phoma and Verticillium. The objectives of the study were to examine and explore the patterns of morphological, physiological and metabolic diversity in these B. napus cultivars at different sulfur concentrations and in the context of plant defense. Results indicate that the root systems are influenced differently by sulfur deficiency in the cultivars. Total root dry mass and length of root hairs differ not only among the cultivars but also vary in their reaction to sulfur limitation and pathogen attack. As a sensitive indicator of stress, several parameters of photosynthetic activity determined by PAM imaging showed a broad variability among the treatments. These results were supported by thermographic analysis. Levels of sulfur-containing metabolites also showed large variations. The data were interrelated to predict the specific behavior during sulfur limitation and/or pathogen attack. Advice for farming are discussed.

Begin at the beginning: A BAC-end view of the passion fruit (Passiflora) genome

Background

The passion fruit (Passiflora edulis) is a tropical crop of economic importance both for juice production and consumption as fresh fruit. The juice is also used in concentrate blends that are consumed worldwide. However, very little is known about the genome of the species. Therefore, improving our understanding of passion fruit genomics is essential and to some degree a pre-requisite if its genetic resources are to be used more efficiently. In this study, we have constructed a large-insert BAC library and provided the first view on the structure and content of the passion fruit genome, using BAC-end sequence (BES) data as a major resource.

Results

The library consisted of 82,944 clones and its levels of organellar DNA were very low. The library represents six haploid genome equivalents, and the average insert size was 108 kb. To check its utility for gene isolation, successful macroarray screening experiments were carried out with probes complementary to eight Passiflora gene sequences available in public databases. BACs harbouring those genes were used in fluorescent in situ hybridizations and unique signals were detected for four BACs in three chromosomes (n = 9). Then, we explored 10,000 BES and we identified reads likely to contain repetitive mobile elements (19.6% of all BES), simple sequence repeats and putative proteins, and to estimate the GC content (~42%) of the reads. Around 9.6% of all BES were found to have high levels of similarity to plant genes and ontological terms were assigned to more than half of the sequences analysed (940). The vast majority of the top-hits made by our sequences were to Populus trichocarpa (24.8% of the total occurrences), Theobroma cacao (21.6%), Ricinus communis (14.3%), Vitis vinifera (6.5%) and Prunus persica (3.8%).

Conclusions

We generated the first large-insert library for a member of Passifloraceae. This BAC library provides a new resource for genetic and genomic studies, as well as it represents a valuable tool for future whole genome study. Remarkably, a number of BAC-end pair sequences could be mapped to intervals of the sequenced Arabidopsis thaliana, V. vinifera and P. trichocarpa chromosomes, and putative collinear microsyntenic regions were identified.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-816) contains supplementary material, which is available to authorized users.

Deep proteomics of the Xenopus laevis egg using an mRNA-derived reference database

BACKGROUND

Mass spectrometry-based proteomics enables the global identification and quantification of proteins and their posttranslational modifications in complex biological samples. However, proteomic analysis requires a complete and accurate reference set of proteins and is therefore largely restricted to model organisms with sequenced genomes.

RESULTS

Here, we demonstrate the feasibility of deep genome-free proteomics by using a reference proteome derived from heterogeneous mRNA data. We identify more than 11,000 proteins with 99% confidence from the unfertilized Xenopus laevis egg and estimate protein abundance with approximately 2-fold precision. Our reference database outperforms the provisional gene models based on genomic DNA sequencing and references generated by other methods. Surprisingly, we find that many proteins in the egg lack mRNA support and that many of these proteins are found in blood or liver, suggesting that they are taken up from the blood plasma, together with yolk, during oocyte growth and maturation, potentially contributing to early embryogenesis.

CONCLUSION

To facilitate proteomics in nonmodel organisms, we make our platform available as an online resource that converts heterogeneous mRNA data into a protein reference set. Thus, we demonstrate the feasibility and power of genome-free proteomics while shedding new light on embryogenesis in vertebrates.

Gene duplication followed by exon structure divergence substitutes for alternative splicing in zebrafish

In this study we report novel findings regarding the evolutionary relationship between gene duplication and alternative splicing, two processes that increase proteomic diversity. By studying teleost fish, we find that gene duplication followed by exon structure divergence between paralogs, but not gene duplication alone, leads to a significant reduction in alternative splicing, as measured by both the proportion of genes that undergo alternative splicing as well as mean number of transcripts per gene. Additionally, we show that this effect is independent of gene family size and gene function. Furthermore, we provide evidence that the reduction in alternative splicing may be due to the partitioning of ancestral splice forms among the duplicate genes - a form of subfunctionalization. Taken together these results indicate that exon structure evolution subsequent to gene duplication may be a common substitute for alternative splicing.

Putative sugar transporters of the mustard leaf beetle Phaedon cochleariae: their phylogeny and role for nutrient supply in larval defensive glands

Background

Phytophagous insects have emerged successfully on the planet also because of the development of diverse and often astonishing defensive strategies against their enemies. The larvae of the mustard leaf beetle Phaedon cochleariae, for example, secrete deterrents from specialized defensive glands on their back. The secretion process involves ATP-binding cassette transporters. Therefore, sugar as one of the major energy sources to fuel the ATP synthesis for the cellular metabolism and transport processes, has to be present in the defensive glands. However, the role of sugar transporters for the production of defensive secretions was not addressed until now.

Results

To identify sugar transporters in P. cochleariae, a transcript catalogue was created by Illumina sequencing of cDNA libraries. A total of 68,667 transcripts were identified and 68 proteins were annotated as either members of the solute carrier 2 (SLC2) family or trehalose transporters. Phylogenetic analyses revealed an extension of the mammalian GLUT6/8 class in insects as well as one group of transporters exhibiting distinctive conserved motifs only present in the insect order Coleoptera. RNA-seq data of samples derived from the defensive glands revealed six transcripts encoding sugar transporters with more than 3,000 counts. Two of them are exclusively expressed in the glandular tissue. Reduction in secretions production was accomplished by silencing two of four selected transporters. RNA-seq experiments of transporter-silenced larvae showed the down-regulation of the silenced transporter but concurrently the up-regulation of other SLC2 transporters suggesting an adaptive system to maintain sugar homeostasis in the defensive glands.

Conclusion

We provide the first comprehensive phylogenetic study of the SLC2 family in a phytophagous beetle species. RNAi and RNA-seq experiments underline the importance of SLC2 transporters in defensive glands to achieve a chemical defense for successful competitive interaction in natural ecosystems.

Establishing genomic tools and resources for Guizotia abyssinica (L.f.) Cass.-the development of a library of expressed sequence tags, microsatellite loci, and the sequencing of its chloroplast genome

We present an EST library, chloroplast genome sequence, and nuclear microsatellite markers that were developed for the semi-domesticated oilseed crop noug (Guizotia abyssinica) from Ethiopia. The EST library consists of 25 711 Sanger reads, assembled into 17 538 contigs and singletons, of which 4781 were functionally annotated using the Arabidopsis Information Resource (TAIR). The age distribution of duplicated genes in the EST library shows evidence of two paleopolyploidizations-a pattern that noug shares with several other species in the Heliantheae tribe (Compositae family). From the EST library, we selected 43 microsatellites and then designed and tested primers for their amplification. The number of microsatellite alleles varied between 2 and 10 (average 4.67), and the average observed and expected heterozygosities were 0.49 and 0.54, respectively. The chloroplast genome was sequenced de novo using Illumina's sequencing technology and completed with traditional Sanger sequencing. No large re-arrangements were found between the noug and sunflower chloroplast genomes, but 1.4% of sites have indels and 1.8% show sequence divergence between the two species. We identified 34 tRNAs, 4 rRNA sequences, and 80 coding sequences, including one region (trnH-psbA) with 15% sequence divergence between noug and sunflower that may be particularly useful for phylogeographic studies in noug and its wild relatives.

Bioinformatics and Gene Network Analyses of the Swine Mammary Gland Transcriptome during Late Gestation

We used the newly-developed Dynamic Impact Approach (DIA) and gene network analysis to study the sow mammary transcriptome at 80, 100, and 110 days of pregnancy. A swine oligoarray with 13,290 inserts was used for transcriptome profiling. An ANOVA with false discovery rate (FDR < 0.15) correction resulted in 1,409 genes with a significant time effect across time comparisons. The DIA uncovered that Fatty acid biosynthesis, Interleukin-4 receptor binding, Galactose metabolism, and mTOR signaling were among the most-impacted pathways. IL-4 receptor binding, ABC transporters, cytokine-cytokine receptor interaction, and Jak-STAT signaling were markedly activated at 110 days compared with 80 and 100 days. Epigenetic and transcription factor regulatory mechanisms appear important in coordinating the final stages of mammary development during pregnancy. Network analysis revealed a crucial role for TP53, ARNT2, E2F4, and PPARG. The bioinformatics analyses revealed a number of pathways and functions that perform an irreplaceable role during late gestation to farrowing.

High-throughput RNA sequencing profiles and transcriptional evidence of aerobic respiratory enzymes in sporulating oocysts and sporozoites of Eimeria tenella

Seven species of Eimeria are responsible for coccidiosis in chickens. Eimeria tenella is one of the most pathogenic parasites since it is associated with high mortality and great economic impact. The life cycle of the parasite includes development in the environment and in the intestinal tract. We conducted RNA sequencing using a next generation sequencer to obtain transcriptome information from the sporulating oocysts, and sporozoites. We collected 2.8 million 75 bp reads of a short-tag sequence, and 25,880 contigs were generated by the Oases assembler. A Blastx search of GenBank databases revealed that 7780 contigs (30.1%) had significant homology with deposited sequence data (E-value <1e-6); among these contigs, 6051 contigs were similar to those of Toxoplasma gondii while only 513 contigs (6.6%) were similar to those of E. tenella. After an orthological analysis conducted with the UniProt database of T. gondii, 6661 contigs were distributed within the categories of cellular components (1528 gene categories), biological processes (861 gene categories), and molecular functions (241 gene categories). The significantly matched contigs contained high numbers of enzymes associated with glycolysis, TCA, and the pentose-phosphate pathway. Most of the enzymes, measured by quantitative reverse transcription-PCR, were up-regulated in sporulating stage. These results suggest that the intracellular carbohydrate amylopectin could be used as an energy source for ATP production including glycolysis and the pentose-phosphate pathway, which generates NADPH and pentoses. Our data also suggest that Eimeria might possess a partial or similar pathway to the TCA cycle essential for aerobic respiration. Furthermore, the newly annotated and non-annotated contigs might contain E. tenella-specific or novel sequences.

Transcriptome-scale homoeolog-specific transcript assemblies of bread wheat

Background

Bread wheat is one of the world’s most important food crops and considerable efforts have been made to develop genomic resources for this species. This includes an on-going project by the International Wheat Genome Sequencing Consortium to assemble its large and complex genome, which is hexaploid and contains three closely related ‘homoeologous’ copies for each chromosome. This multi-national effort avoids the complications polyploidy entails for correct assembly of the genome by sequencing flow-sorted chromosome arms one at a time. Here we report on an alternate approach, a direct homoeolog-specific assembly of the expressed portion of the genome, the transcriptome.

Results

After assessment of the ability of various assemblers to generate homoeolog-specific assemblies, we employed a two-stage assembly process to produce a high-quality assembly of the transcriptome of hexaploid wheat from Roche-454 and Illumina GAII_x paired-end sequence reads. The assembly process made use of a rapid partitioning of expressed sequences into homoeologous clusters, followed by a parallel high-fidelity assembly of each cluster on a 1150-processor compute cloud. We assessed assembly quality through comparison to known wheat gene sequences and found that in ca. 98.5% of cases the assembly was sufficiently accurate for homoeologous triplets to be cleanly separated into either two or three separate contigs. Comparison to publicly available transcript collections suggests that the assembly covers ~75-80% of the complete transcriptome.

Conclusions

This work therefore describes the first homoeolog-specific sequence assembly of the wheat transcriptome and provides a reference transcriptome for future wheat research. Furthermore, our assembly methodology is transferable to other polyploid organisms.

In silico identification of known osmotic stress responsive genes from Arabidopsis in soybean and Medicago

Plants experience various environmental stresses, but tolerance to these adverse conditions is a very complex phenomenon. The present research aimed to evaluate a set of genes involved in osmotic response, comparing soybean and medicago with the well-described Arabidopsis thaliana model plant. Based on 103 Arabidopsis proteins from 27 categories of osmotic stress response, comparative analyses against Genosoja and Medicago truncatula databases allowed the identification of 1,088 soybean and 1,210 Medicago sequences. The analysis showed a high number of sequences and high diversity, comprising genes from all categories in both organisms. Genes with unknown function were among the most representative, followed by transcription factors, ion transport proteins, water channel, plant defense, protein degradation, cellular structure, organization & biogenesis and senescence. An analysis of sequences with unknown function allowed the annotation of 174 soybean and 217 Medicago sequences, most of them concerning transcription factors. However, for about 30% of the sequences no function could be attributed using in silico procedures. The establishment of a gene set involved in osmotic stress responses in soybean and barrel medic will help to better understand the survival mechanisms for this type of stress condition in legumes.

A second generation framework for the analysis of microsatellites in expressed sequence tags and the development of EST-SSR markers for a conifer, Cryptomeria japonica

Background

Microsatellites or simple sequence repeats (SSRs) in expressed sequence tags (ESTs) are useful resources for genome analysis because of their abundance, functionality and polymorphism. The advent of commercial second generation sequencing machines has lead to new strategies for developing EST-SSR markers, necessitating the development of bioinformatic framework that can keep pace with the increasing quality and quantity of sequence data produced. We describe an open scheme for analyzing ESTs and developing EST-SSR markers from reads collected by Sanger sequencing and pyrosequencing of sugi (Cryptomeria japonica).

Results

We collected 141,097 sequence reads by Sanger sequencing and 1,333,444 by pyrosequencing. After trimming contaminant and low quality sequences, 118,319 Sanger and 1,201,150 pyrosequencing reads were passed to the MIRA assembler, generating 81,284 contigs that were analysed for SSRs. 4,059 SSRs were found in 3,694 (4.54%) contigs, giving an SSR frequency lower than that in seven other plant species with gene indices (5.4–21.9%). The average GC content of the SSR-containing contigs was 41.55%, compared to 40.23% for all contigs. Tri-SSRs were the most common SSRs; the most common motif was AT, which was found in 655 (46.3%) di-SSRs, followed by the AAG motif, found in 342 (25.9%) tri-SSRs. Most (72.8%) tri-SSRs were in coding regions, but 55.6% of the di-SSRs were in non-coding regions; the AT motif was most abundant in 3′ untranslated regions. Gene ontology (GO) annotations showed that six GO terms were significantly overrepresented within SSR-containing contigs. Forty–four EST-SSR markers were developed from 192 primer pairs using two pipelines: read2Marker and the newly-developed CMiB, which combines several open tools. Markers resulting from both pipelines showed no differences in PCR success rate and polymorphisms, but PCR success and polymorphism were significantly affected by the expected PCR product size and number of SSR repeats, respectively. EST-SSR markers exhibited less polymorphism than genomic SSRs.

Conclusions

We have created a new open pipeline for developing EST-SSR markers and applied it in a comprehensive analysis of EST-SSRs and EST-SSR markers in C. japonica. The results will be useful in genomic analyses of conifers and other non-model species.

How can microarrays unlock asthma?

Asthma is a complex disease regulated by the interplay of a large number of underlying mechanisms which contribute to the overall pathology. Despite various breakthroughs identifying genes related to asthma, our understanding of the importance of the genetic background remains limited. Although current therapies for asthma are relatively effective, subpopulations of asthmatics do not respond to these regimens. By unlocking the role of these underlying mechanisms, a source of novel and more effective treatments may be identified. In the new age of high-throughput technologies, gene-expression microarrays provide a quick and effective method of identifying novel genes and pathways, which would be impossible to discover using an individual gene screening approach. In this review we follow the history of expression microarray technologies and describe their contributions to advancing our current knowledge and understanding of asthma pathology.

Changes in gene expression in a porcine preadipocyte cell line during differentiation

Adipocyte differentiation plays an important role in the formation of fat tissues in pigs and affects meat quality and productivity. Clarification of the nature of the pig genes that participate in adipocyte differentiation will provide a clue to the regulation of fat content and thickness in pig carcases by dietary control; it will also help to find target genes for exploring potentially useful polymorphisms for molecular breeding aimed at fat traits. We constructed a DNA oligomer microarray based on pig transcripts, and we used the array to investigate time-dependent changes in gene expression in the PSPA porcine preadipocyte cell line during differentiation into adipocytes. We selected genes with markedly altered expression (at least fivefold difference in comparison with expression in undifferentiated cells) and classified them into five groups according to gene expression pattern. In the early stage after stimulation of adipocyte differentiation, we observed up-regulation of many genes encoding proteins involved in regulating cell proliferation and transcription. Among the probes corresponding to transcripts that showed marked changes in expression, 27 were located within previously reported QTL regions for traits related to adipose tissues. These results will be valuable resources for finding the genes responsible for fat-related traits that have been identified in previous studies using various pig resource families.

RNA-seq discovery, functional characterization, and comparison of sesquiterpene synthases from Solanum lycopersicum and Solanum habrochaites trichomes

Solanum lycopersicum and Solanum habrochaites (f. typicum) accession PI127826 emit a variety of sesquiterpenes. To identify terpene synthases involved in the production of these volatile sesquiterpenes, we used massive parallel pyrosequencing (RNA-seq) to obtain the transcriptome of the stem trichomes from these plants. This approach resulted initially in the discovery of six sesquiterpene synthase cDNAs from S. lycopersicum and five from S. habrochaites. Searches of other databases and the S. lycopersicum genome resulted in the discovery of two additional sesquiterpene synthases expressed in trichomes. The sesquiterpene synthases from S. lycopersicum and S. habrochaites have high levels of protein identity. Several of them appeared to encode for non-functional proteins. Functional recombinant proteins produced germacrenes, β-caryophyllene/α-humulene, viridiflorene and valencene from (E,E)-farnesyl diphosphate. However, the activities of these enzymes do not completely explain the differences in sesquiterpene production between the two tomato plants. RT-qPCR confirmed high levels of expression of most of the S. lycopersicum sesquiterpene synthases in stem trichomes. In addition, one sesquiterpene synthase was induced by jasmonic acid, while another appeared to be slightly repressed by the treatment. Our data provide a foundation to study the evolution of terpene synthases in cultivated and wild tomato.

EST dataset of pituitary and identification of somatolactin and novel genes in Chinese sturgeon, Acipenser sinensis

Chinese sturgeon (Acipenser sinensis) is a rare and endangered species and also an important resource for the sturgeon aquaculture industry, however, a few genes have been identified in this species. We report here construction of a pituitary cDNA library from a 24 years old female Chinese sturgeon just after its spawning, and obtained 2,025 ESTs from the library. 885 unique sequences were identified, which were categorized into 12 functional groups. More than half of the unique sequences (57%) do not match with annotated sequences in the public databases. Three of these novel genes were further identified. Notably, a full-length of cDNA (1,143 bp) encoding somatolactin of 232 amino acids was identified. Phylogenetic analysis showed 97% amino acid identity with White sturgeon somatolactin. RT-PCR analysis indicated that the somatolactin mRNA was only detected in pituitary. Pituitary-specific expression of the somatolactin suggested that the protein may play important physiological functions in pituitary-endocrine system of the Chinese sturgeon.

Analysis of expressed sequence tags from the marine microalga Pseudochattonella farcimen (Dictyochophyceae)

Pseudochattonella farcimen (Eikrem, Edvardsen, et Throndsen) is a unicellular alga belonging to the Dictyochophyceae (Heterokonta). It forms recurring blooms in Scandinavian coastal waters, and has been associated to fish mortality. Here we report the sequencing and analysis of 10,368 expressed sequence tags (ESTs) corresponding to 8,149 unique gene models from this species. Compared to EST libraries from other heterokonts, P. farcimen contains a high number of genes with functions related to cell communication and signaling. We found several genes encoding proteins related to fatty acid metabolism, including eight fatty acid desaturases and two phospholipase A2 genes. Three desaturases are highly similar to Δ4-desaturases from haptophytes. P. farcimen also possesses three putative polyketide synthases (PKSs), belonging to two different families. Some of these genes may have been acquired via horizontal gene transfer by a common ancestor of brown algae and dictyochophytes, together with genes involved in mannitol metabolism, which are also present in P. farcimen. Our findings may explain the unusual fatty acid profile previously observed in P. farcimen, and are discussed from an evolutionary perspective and in relation to the ichthyotoxicity of this alga.

Global analysis of proline-rich tandem repeat proteins reveals broad phylogenetic diversity in plant secretomes

Cell walls, constructed by precisely choreographed changes in the plant secretome, play critical roles in plant cell physiology and development. Along with structural polysaccharides, secreted proline-rich Tandem Repeat Proteins (TRPs) are important for cell wall function, yet the evolutionary diversity of these structural TRPs remains virtually unexplored. Using a systems-level computational approach to analyze taxonomically diverse plant sequence data, we identified 31 distinct Pro-rich TRP classes targeted for secretion. This analysis expands upon the known phylogenetic diversity of extensins, the most widely studied class of wall structural proteins, and demonstrates that extensins evolved before plant vascularization. Our results also show that most Pro-rich TRP classes have unexpectedly restricted evolutionary distributions, revealing considerable differences in plant secretome signatures that define unexplored diversity.

Molecular adaptation in flowering and symbiotic recognition pathways: insights from patterns of polymorphism in the legume Medicago truncatula

Background

We studied patterns of molecular adaptation in the wild Mediterranean legume Medicago truncatula. We focused on two phenotypic traits that are not functionally linked: flowering time and perception of symbiotic microbes. Phenology is an important fitness component, especially for annual plants, and many instances of molecular adaptation have been reported for genes involved in flowering pathways. While perception of symbiotic microbes is also integral to adaptation in many plant species, very few reports of molecular adaptation exist for symbiotic genes. Here we used data from 57 individuals and 53 gene fragments to quantify the overall strength of both positive and purifying selection in M. truncatula and asked if footprints of positive selection can be detected at key genes of rhizobia recognition pathways.

Results

We examined nucleotide variation among 57 accessions from natural populations in 53 gene fragments: 5 genes involved in nitrogen-fixing bacteria recognition, 11 genes involved in flowering, and 37 genes used as control loci. We detected 1757 polymorphic sites yielding an average nucleotide diversity (pi) of 0.003 per site. Non-synonymous variation is under sizable purifying selection with 90% of amino-acid changing mutations being strongly selected against. Accessions were structured in two groups consistent with geographical origins. Each of these two groups harboured an excess of rare alleles, relative to expectations of a constant-sized population, suggesting recent population expansion. Using coalescent simulations and an approximate Bayesian computation framework we detected several instances of genes departing from selective neutrality within each group and showed that the polymorphism of two nodulation and four flowering genes has probably been shaped by recent positive selection.

Conclusion

We quantify the intensity of purifying selection in the M. truncatula genome and show that putative footprints of natural selection can be detected at different time scales in both flowering and symbiotic pathways.

Bioinformatic analysis of ESTs collected by Sanger and pyrosequencing methods for a keystone forest tree species: oak

BACKGROUND

The Fagaceae family comprises about 1,000 woody species worldwide. About half belong to the Quercus family. These oaks are often a source of raw material for biomass wood and fiber. Pedunculate and sessile oaks, are among the most important deciduous forest tree species in Europe. Despite their ecological and economical importance, very few genomic resources have yet been generated for these species. Here, we describe the development of an EST catalogue that will support ecosystem genomics studies, where geneticists, ecophysiologists, molecular biologists and ecologists join their efforts for understanding, monitoring and predicting functional genetic diversity.

RESULTS

We generated 145,827 sequence reads from 20 cDNA libraries using the Sanger method. Unexploitable chromatograms and quality checking lead us to eliminate 19,941 sequences. Finally a total of 125,925 ESTs were retained from 111,361 cDNA clones. Pyrosequencing was also conducted for 14 libraries, generating 1,948,579 reads, from which 370,566 sequences (19.0%) were eliminated, resulting in 1,578,192 sequences. Following clustering and assembly using TGICL pipeline, 1,704,117 EST sequences collapsed into 69,154 tentative contigs and 153,517 singletons, providing 222,671 non-redundant sequences (including alternative transcripts). We also assembled the sequences using MIRA and PartiGene software and compared the three unigene sets. Gene ontology annotation was then assigned to 29,303 unigene elements. Blast search against the SWISS-PROT database revealed putative homologs for 32,810 (14.7%) unigene elements, but more extensive search with Pfam, Refseq_protein, Refseq_RNA and eight gene indices revealed homology for 67.4% of them. The EST catalogue was examined for putative homologs of candidate genes involved in bud phenology, cuticle formation, phenylpropanoids biosynthesis and cell wall formation. Our results suggest a good coverage of genes involved in these traits. Comparative orthologous sequences (COS) with other plant gene models were identified and allow to unravel the oak paleo-history. Simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) were searched, resulting in 52,834 SSRs and 36,411 SNPs. All of these are available through the Oak Contig Browser http://genotoul-contigbrowser.toulouse.inra.fr:9092/Quercus_robur/index.html.

CONCLUSIONS

This genomic resource provides a unique tool to discover genes of interest, study the oak transcriptome, and develop new markers to investigate functional diversity in natural populations.

Genomic tools development for Aquilegia: construction of a BAC-based physical map

BACKGROUND

The genus Aquilegia, consisting of approximately 70 taxa, is a member of the basal eudicot lineage, Ranuculales, which is evolutionarily intermediate between monocots and core eudicots, and represents a relatively unstudied clade in the angiosperm phylogenetic tree that bridges the gap between these two major plant groups. Aquilegia species are closely related and their distribution covers highly diverse habitats. These provide rich resources to better understand the genetic basis of adaptation to different pollinators and habitats that in turn leads to rapid speciation. To gain insights into the genome structure and facilitate gene identification, comparative genomics and whole-genome shotgun sequencing assembly, BAC-based genomics resources are of crucial importance.

RESULTS

BAC-based genomic resources, including two BAC libraries, a physical map with anchored markers and BAC end sequences, were established from A. formosa. The physical map was composed of a total of 50,155 BAC clones in 832 contigs and 3939 singletons, covering 21X genome equivalents. These contigs spanned a physical length of 689.8 Mb (~2.3X of the genome) suggesting the complex heterozygosity of the genome. A set of 197 markers was developed from ESTs induced by drought-stress, or involved in anthocyanin biosynthesis or floral development, and was integrated into the physical map. Among these were 87 genetically mapped markers that anchored 54 contigs, spanning 76.4 Mb (25.5%) across the genome. Analysis of a selection of 12,086 BAC end sequences (BESs) from the minimal tiling path (MTP) allowed a preview of the Aquilegia genome organization, including identification of transposable elements, simple sequence repeats and gene content. Common repetitive elements previously reported in both monocots and core eudicots were identified in Aquilegia suggesting the value of this genome in connecting the two major plant clades. Comparison with sequenced plant genomes indicated a higher similarity to grapevine (Vitis vinifera) than to rice and Arabidopsis in the transcriptomes.

CONCLUSIONS

The A. formosa BAC-based genomic resources provide valuable tools to study Aquilegia genome. Further integration of other existing genomics resources, such as ESTs, into the physical map should enable better understanding of the molecular mechanisms underlying adaptive radiation and elaboration of floral morphology.

Generation and analysis of ESTs from strawberry (Fragaria xananassa) fruits and evaluation of their utility in genetic and molecular studies

BACKGROUND

Cultivated strawberry is a hybrid octoploid species (Fragaria xananassa Duchesne ex. Rozier) whose fruit is highly appreciated due to its organoleptic properties and health benefits. Despite recent studies on the control of its growth and ripening processes, information about the role played by different hormones on these processes remains elusive. Further advancement of this knowledge is hampered by the limited sequence information on genes from this species, despite the abundant information available on genes from the wild diploid relative Fragaria vesca. However, the diploid species, or one ancestor, only partially contributes to the genome of the cultivated octoploid. We have produced a collection of expressed sequence tags (ESTs) from different cDNA libraries prepared from different fruit parts and developmental stages. The collection has been analysed and the sequence information used to explore the involvement of different hormones in fruit developmental processes, and for the comparison of transcripts in the receptacle of ripe fruits of diploid and octoploid species. The study is particularly important since the commercial fruit is indeed an enlarged flower receptacle with the true fruits, the achenes, on the surface and connected through a network of vascular vessels to the central pith.

RESULTS

We have sequenced over 4,500 ESTs from Fragaria xananassa, thus doubling the number of ESTs available in the GenBank of this species. We then assembled this information together with that available from F. xananassa resulting a total of 7,096 unigenes. The identification of SSRs and SNPs in many of the ESTs allowed their conversion into functional molecular markers. The availability of libraries prepared from green growing fruits has allowed the cloning of cDNAs encoding for genes of auxin, ethylene and brassinosteroid signalling processes, followed by expression studies in selected fruit parts and developmental stages. In addition, the sequence information generated in the project, jointly with previous information on sequences from both F. xananassa and F. vesca, has allowed designing an oligo-based microarray that has been used to compare the transcriptome of the ripe receptacle of the diploid and octoploid species. Comparison of the transcriptomes, grouping the genes by biological processes, points to differences being quantitative rather than qualitative.

CONCLUSIONS

The present study generates essential knowledge and molecular tools that will be useful in improving investigations at the molecular level in cultivated strawberry (F. xananassa). This knowledge is likely to provide useful resources in the ongoing breeding programs. The sequence information has already allowed the development of molecular markers that have been applied to germplasm characterization and could be eventually used in QTL analysis. Massive transcription analysis can be of utility to target specific genes to be further studied, by their involvement in the different plant developmental processes.

Ancient genome duplications during the evolution of kiwifruit (Actinidia) and related Ericales

BACKGROUND AND AIMS

To assess the number and phylogenetic distribution of large-scale genome duplications in the ancestry of Actinidia, publicly available expressed sequenced tags (ESTs) for members of the Actinidiaceae and related Ericales, including tea (Camellia sinensis), were analysed.

METHODS

Synonymous divergences (K(s)) were calculated for all duplications within gene families and examined for evidence of large-scale duplication events. Phylogenetic comparisons for a selection of orthologues among several related species in Ericales and two outgroups permitted placement of duplication events in relation to lineage divergences. Gene ontology (GO) categories were analysed for each whole-genome duplication (WGD) and the whole transcriptome.

KEY RESULTS

Evidence for three ancient WGDs in Actinidia was found. Analyses of paleologue GO categories indicated a different pattern of retained genes for each genome duplication, but a pattern consistent with the dosage-balance hypothesis among all retained paleologues.

CONCLUSIONS

This study provides evidence for one independent WGD in the ancestry of Actinidia (Ad-alpha), a WGD shared by Actinidia and Camellia (Ad-beta), and the well-established At-gamma WGD that occurred prior to the divergence of all taxa examined. More ESTs in other taxa are needed to elucidate which groups in Ericales share the Ad-beta or Ad-alpha duplications and their impact on diversification.

RNA-Seq Atlas of Glycine max: a guide to the soybean transcriptome

BACKGROUND

Next generation sequencing is transforming our understanding of transcriptomes. It can determine the expression level of transcripts with a dynamic range of over six orders of magnitude from multiple tissues, developmental stages or conditions. Patterns of gene expression provide insight into functions of genes with unknown annotation.

RESULTS

The RNA Seq-Atlas presented here provides a record of high-resolution gene expression in a set of fourteen diverse tissues. Hierarchical clustering of transcriptional profiles for these tissues suggests three clades with similar profiles: aerial, underground and seed tissues. We also investigate the relationship between gene structure and gene expression and find a correlation between gene length and expression. Additionally, we find dramatic tissue-specific gene expression of both the most highly-expressed genes and the genes specific to legumes in seed development and nodule tissues. Analysis of the gene expression profiles of over 2,000 genes with preferential gene expression in seed suggests there are more than 177 genes with functional roles that are involved in the economically important seed filling process. Finally, the Seq-atlas also provides a means of evaluating existing gene model annotations for the Glycine max genome.

CONCLUSIONS

This RNA-Seq atlas extends the analyses of previous gene expression atlases performed using Affymetrix GeneChip technology and provides an example of new methods to accommodate the increase in transcriptome data obtained from next generation sequencing. Data contained within this RNA-Seq atlas of Glycine max can be explored at http://www.soybase.org/soyseq.

The ANISEED database: digital representation, formalization, and elucidation of a chordate developmental program

Developmental biology aims to understand how the dynamics of embryonic shapes and organ functions are encoded in linear DNA molecules. Thanks to recent progress in genomics and imaging technologies, systemic approaches are now used in parallel with small-scale studies to establish links between genomic information and phenotypes, often described at the subcellular level. Current model organism databases, however, do not integrate heterogeneous data sets at different scales into a global view of the developmental program. Here, we present a novel, generic digital system, NISEED, and its implementation, ANISEED, to ascidians, which are invertebrate chordates suitable for developmental systems biology approaches. ANISEED hosts an unprecedented combination of anatomical and molecular data on ascidian development. This includes the first detailed anatomical ontologies for these embryos, and quantitative geometrical descriptions of developing cells obtained from reconstructed three-dimensional (3D) embryos up to the gastrula stages. Fully annotated gene model sets are linked to 30,000 high-resolution spatial gene expression patterns in wild-type and experimentally manipulated conditions and to 528 experimentally validated cis-regulatory regions imported from specialized databases or extracted from 160 literature articles. This highly structured data set can be explored via a Developmental Browser, a Genome Browser, and a 3D Virtual Embryo module. We show how integration of heterogeneous data in ANISEED can provide a system-level understanding of the developmental program through the automatic inference of gene regulatory interactions, the identification of inducing signals, and the discovery and explanation of novel asymmetric divisions.

Comparative microarray analysis of Rhipicephalus (Boophilus) microplus expression profiles of larvae pre-attachment and feeding adult female stages on Bos indicus and Bos taurus cattle

Background

Rhipicephalus (Boophilus) microplus is an obligate blood feeder which is host specific to cattle. Existing knowledge pertaining to the host or host breed effects on tick transcript expression profiles during the tick - host interaction is poor.

Results

Global analysis of gene expression changes in whole R. microplus ticks during larval, pre-attachment and early adult stages feeding on Bos indicus and Bos taurus cattle were compared using gene expression microarray analysis. Among the 13,601 R. microplus transcripts from BmiGI Version 2 we identified 297 high and 17 low expressed transcripts that were significantly differentially expressed between R. microplus feeding on tick resistant cattle [Bos indicus (Brahman)] compared to R. microplus feeding on tick susceptible cattle [Bos taurus (Holstein-Friesian)] (p ≤ 0.001). These include genes encoding enzymes involved in primary metabolism, and genes related to stress, defence, cell wall modification, cellular signaling, receptor, and cuticle formation. Microarrays were validated by qRT-PCR analysis of selected transcripts using three housekeeping genes as normalization controls.

Conclusion

The analysis of all tick stages under survey suggested a coordinated regulation of defence proteins, proteases and protease inhibitors to achieve successful attachment and survival of R. microplus on different host breeds, particularly Bos indicus cattle. R. microplus ticks demonstrate different transcript expression patterns when they encounter tick resistant and susceptible breeds of cattle. In this study we provide the first transcriptome evidence demonstrating the influence of tick resistant and susceptible cattle breeds on transcript expression patterns and the molecular physiology of ticks during host attachment and feeding.

The microarray data used in this analysis have been submitted to NCBI GEO database under accession number GSE20605 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE20605.

Using the DFCI gene index databases for biological discovery

The DFCI Gene Index Web pages provide access to analyses of ESTs and gene sequences for nearly 114 species, as well as a number of resources derived from these. Each species-specific database is presented using a common format with a home page. A variety of methods exist that allow users to search each species-specific database. Methods implemented currently include nucleotide or protein sequence queries using WU-BLAST, text-based searches using various sequence identifiers, searches by gene, tissue and library name, and searches using functional classes through Gene Ontology assignments. This protocol provides guidance for using the Gene Index Databases to extract information.

Dating and functional characterization of duplicated genes in the apple (Malus domestica Borkh.) by analyzing EST data

BACKGROUND

Gene duplication is central to genome evolution. In plants, genes can be duplicated through small-scale events and large-scale duplications often involving polyploidy. The apple belongs to the subtribe Pyrinae (Rosaceae), a diverse lineage that originated via allopolyploidization. Both small-scale duplications and polyploidy may have been important mechanisms shaping the genome of this species.

RESULTS

This study evaluates the gene duplication and polyploidy history of the apple by characterizing duplicated genes in this species using EST data. Overall, 68% of the apple genes were clustered into families with a mean copy-number of 4.6. Analysis of the age distribution of gene duplications supported a continuous mode of small-scale duplications, plus two episodes of large-scale duplicates of vastly different ages. The youngest was consistent with the polyploid origin of the Pyrinae 37-48 MYBP, whereas the older may be related to gamma-triplication; an ancient hexapolyploidization previously characterized in the four sequenced eurosid genomes and basal to the eurosid-asterid divergence. Duplicated genes were studied for functional diversification with an emphasis on young paralogs; those originated during or after the formation of the Pyrinae lineage. Unequal assignment of single-copy genes and gene families to Gene Ontology categories suggested functional bias in the pattern of gene retention of paralogs. Young paralogs related to signal transduction, metabolism, and energy pathways have been preferentially retained. Non-random retention of duplicated genes seems to have mediated the expansion of gene families, some of which may have substantially increased their members after the origin of the Pyrinae. The joint analysis of over-duplicated functional categories and phylogenies, allowed evaluation of the role of both polyploidy and small-scale duplications during this process. Finally, gene expression analysis indicated that 82% of duplicated genes, including 80% of young paralogs, showed uncorrelated expression profiles, suggesting extensive subfunctionalization and a role of gene duplication in the acquisition of novel patterns of gene expression.

CONCLUSIONS

This study reports a genome-wide analysis of the mode of gene duplication in the apple, and provides evidence for its role in genome functional diversification by characterising three major processes: selective retention of paralogs, amplification of gene families, and changes in gene expression.

Characterization of duplicate gene evolution in the recent natural allopolyploid Tragopogon miscellus by next-generation sequencing and Sequenom iPLEX MassARRAY genotyping

Tragopogon miscellus (Asteraceae) is an evolutionary model for the study of natural allopolyploidy, but until now has been under-resourced as a genetic model. Using 454 and Illumina expressed sequence tag sequencing of the parental diploid species of T. miscellus, we identified 7782 single nucleotide polymorphisms that differ between the two progenitor genomes present in this allotetraploid. Validation of a sample of 98 of these SNPs in genomic DNA using Sequenom MassARRAY iPlex genotyping confirmed 92 SNP markers at the genomic level that were diagnostic for the two parental genomes. In a transcriptome profile of 2989 SNPs in a single T. miscellus leaf, using Illumina sequencing, 69% of SNPs showed approximately equal expression of both homeologs (duplicate homologous genes derived from different parents), 22% showed apparent differential expression and 8.5% showed apparent silencing of one homeolog in T. miscellus. The majority of cases of homeolog silencing involved the T. dubius SNP homeolog (164/254; 65%) rather than the T. pratensis homeolog (90/254). Sequenom analysis of genomic DNA showed that in a sample of 27 of the homeologs showing apparent silencing, 23 (85%) were because of genomic homeolog loss. These methods could be applied to any organism, allowing efficient and cost-effective generation of genetic markers.

SNP identification, verification, and utility for population genetics in a non-model genus

Background

By targeting SNPs contained in both coding and non-coding areas of the genome, we are able to identify genetic differences and characterize genome-wide patterns of variation among individuals, populations and species. We investigated the utility of 454 sequencing and MassARRAY genotyping for population genetics in natural populations of the teleost, Fundulus heteroclitus as well as closely related Fundulus species (F. grandis, F. majalis and F. similis).

Results

We used 454 pyrosequencing and MassARRAY genotyping technology to identify and type 458 genome-wide SNPs and determine genetic differentiation within and between populations and species of Fundulus. Specifically, pyrosequencing identified 96 putative SNPs across coding and non-coding regions of the F. heteroclitus genome: 88.8% were verified as true SNPs with MassARRAY. Additionally, putative SNPs identified in F. heteroclitus EST sequences were verified in most (86.5%) F. heteroclitus individuals; fewer were genotyped in F. grandis (74.4%), F. majalis (72.9%), and F. similis (60.7%) individuals. SNPs were polymorphic and showed latitudinal clinal variation separating northern and southern populations and established isolation by distance in F. heteroclitus populations. In F. grandis, SNPs were less polymorphic but still established isolation by distance. Markers differentiated species and populations.

Conclusions

In total, these approaches were used to quickly determine differences within the Fundulus genome and provide markers for population genetic studies.

wDBTF: an integrated database resource for studying wheat transcription factor families

BACKGROUND

Transcription factors (TFs) regulate gene expression by interacting with promoters of their target genes and are classified into families based on their DNA-binding domains. Genes coding for TFs have been identified in the sequences of model plant genomes. The rice (Oryza sativa spp. japonica) genome contains 2,384 TF gene models, which represent the mRNA transcript of a locus, classed into 63 families.

RESULTS

We have created an extensive list of wheat (Triticum aestivum L) TF sequences based on sequence homology with rice TFs identified and classified in the Database of Rice Transcription Factors (DRTF). We have identified 7,112 wheat sequences (contigs and singletons) from a dataset of 1,033,960 expressed sequence tag and mRNA (ET) sequences available. This number is about three times the number of TFs in rice so proportionally is very similar if allowance is made for the hexaploidy of wheat. Of these sequences 3,820 encode gene products with a DNA-binding domain and thus were confirmed as potential regulators. These 3,820 sequences were classified into 40 families and 84 subfamilies and some members defined orphan families. The results were compiled in the Database of Wheat Transcription Factor (wDBTF), an inventory available on the web http://wwwappli.nantes.inra.fr:8180/wDBFT/. For each accession, a link to its library source and its Affymetrix identification number is provided. The positions of Pfam (protein family database) motifs were given when known.

CONCLUSIONS

wDBTF collates 3,820 wheat TF sequences validated by the presence of a DNA-binding domain out of 7,112 potential TF sequences identified from publicly available gene expression data. We also incorporated in silico expression data on these TFs into the database. Thus this database provides a major resource for systematic studies of TF families and their expression in wheat as illustrated here in a study of DOF family members expressed during seed development.

The sulfate transporter family in wheat: tissue-specific gene expression in relation to nutrition

Sulfate uptake and distribution in plants are managed by the differential expression of a family of transporters, developmentally, spatially, and in response to sulfur nutrition. Elucidation of the signaling pathways involved requires a knowledge of the component parts and their interactions. Here, the expression patterns of the full complement of sulfate transporters in wheat, as influenced by development and sulfur nutrition, are described. The 10 wheat sulfate transporters characterized here are compared to the gene families for both rice and Brachypodium, for whom full genome information is available. Expression is reported in young seedlings with a focus on roles in uptake from nutrient solution and differential expression in relation to sulfate deprivation. In addition, patterns of expression in all organs at the grain filling stage are reported and indicate differential responses to nutritional signals of the individual transporters in specific tissues and an overall coordination of uptake, storage, and remobilization to deliver sulfur to the developing grain.