Identification of tissue-specific, abiotic stress-responsive gene expression patterns in wine grape (Vitis vinifera L.) based on curation and mining of large-scale EST data sets

Transcriptomic Analysis Elaborates the Resistance Mechanism of Grapevine Rootstocks against Salt Stress

Grapes are subject to a wide range of climatic conditions during their life cycle, but the use of rootstocks can effectively ameliorate the effects of abiotic stress. However, the tolerance mechanism of different grape rootstock varieties varies under various stresses, and systematic research on this aspect is limited. On the basis of previous research, transcriptome sequencing was performed on three tolerant grape rootstock varieties (3309C, 520A, 1103P) and three intolerant grape rootstock varieties (5BB, 101-14, Beta). In total, 56,478,468 clean reads were obtained. One hundred and ten genes only existed in all combinations during P1 with a downregulated trend, and 178 genes existed only in P1 of tolerant grape rootstock varieties. Salt treatment firstly affected the photosynthesis of leaves, and tolerant varieties weakened or even eliminated this effect through their own mechanisms in the later stage. Tolerant varieties mobilized a large number of MFs during the P2 stage, such as hydrolase activity, carboxypeptidase activity, and dioxygenase activity. Carbon metabolism was significantly enriched in P1, while circadian rhythm and flavonoid biosynthesis were only enriched in tolerant varieties. In the intolerant varieties, photosynthesis-related pathways were always the most significantly enriched. There were large differences in the gene expression of the main signal pathways related to salt stress in different varieties. Salt stress affected the expression of genes related to plant abiotic stress, biotic stress, transcription factors, hormones, and secondary metabolism. Tolerant varieties mobilized more bHLH, WRKY, and MYB transcription factors to respond to salt stress than intolerant varieties. In the tolerant rootstocks, SOS was co-expressed. Among these, SOS1 and SOS2 were upregulated, and the SOS3 and SOS5 components were downregulated. The genes of heat shock proteins and the phenylalanine pathway were upregulated in the tolerant varieties. These findings outline a tolerance mechanism model for rootstocks for coping with osmotic stress, providing important information for improving the resistance of grapes under global climate change.

Tracking the time-dependent and tissue-specific processes of arsenic accumulation and stress responses in rice (Oryza sativa L.)

The present study analysed time (0.5 h to 24 h) and tissue [roots, old leaves (OL) and young leaves (YL)] dependent nature of arsenic (As) accumulation and ensuing responses in two contrasting varieties of rice (Oryza sativa L.); Pooja (tolerant) and CO-50 (moderately sensitive). Arsenic accumulation was 5.4-, 4.7- and 7.3-fold higher at 24 h in roots, OL and YL, respectively of var. CO-50 than that in var. Pooja. Arsenic accumulation in YL depicted a delayed accumulation; at 2 h onwards in var. Pooja (0.23 µg g^-1 dw) while at 1 h onwards in var. CO50 (0.26 µg g^-1 dw). The responses of oxidative stress parameters, antioxidant enzymes, metabolites and ions were also found to be tissue- and time-dependent and depicted differential pattern in the two varieties. Among hormone, salicylic acid and abscisic acid showed variable response in var. Pooja and var. CO-50. Metabolite analysis depicted an involvement of various metabolites in As stress responses of two varieties. In conclusion, an early sensing of the As stress, proper coordination of hormones, biochemical responses, ionic and metabolic profiles allowed var. Pooja to resist As stress and reduce As accumulation more effectively as compared to that of var. CO-50.

Identification and validation of new reference genes for accurate quantitative reverse transcriptase-PCR normalization in the Antarctic plant Colobanthus quitensis under abiotic stress conditions

The Antarctic ecotype of Colobanthus quitensis is a vascular plant highly adapted to the harsh environmental conditions of Maritime Antarctica which is now facing with the rapid local warming experienced in the Antarctic Peninsula during the last decades. Thus, the identification of the molecular mechanisms leading to the adaptation to this warming trend is a new target for modern cell physiology. The selection of suitable reference genes for quantification of key stress-responsive genes through quantitative Reverse Transcriptase-Polymerase Chain Reaction (qRT-PCR) is important to ensure accurate and reliable results. In this study, we evaluated the expression stability of eleven candidate genes in C. quitensis under different abiotic stress conditions using geNorm and RefFinder tools. The statistical analysis showed that the appropriate reference genes varied depending on the experimental conditions, even if EF1α and PP2Acs ranked as the most stable reference genes when all stress conditions were considered. To further validate the stability of the selected reference genes, the expression patterns of C. quitensis catalase gene (CqCAT) was analyzed. The reference genes validated in this study will be useful for improving the accuracy of qRT-PCR analysis for gene expression studies of the Antarctic ecotype of C. quitensis and could be extended to other ecotypes adapted to low temperatures.

Salt-Responsive Genes are Differentially Regulated at the Chromatin Levels Between Seedlings and Roots in Rice

The elucidation of epigenetic responses of salt-responsive genes facilitates understanding of the underlying mechanisms that confer salt tolerance in rice. However, it is still largely unknown how epigenetic mechanisms are associated with the expression of salt-responsive genes in rice and other crops. In this study, we reported tissue-specific gene expression and tissue-specific changes in chromatin modifications or signatures between seedlings and roots in response to salt treatment. Our study indicated that among six of individual mark examined (H3K4me3, H3K27me3, H4K12ac, H3K9ac, H3K27ac and H3K36me3), a positive association between salt-related changes in histone marks and the expression of differentially expressed genes (DEGs) was observed only for H3K9ac and H4K12ac in seedlings and H3K36me3 in roots. In contrast, chromatin states (CSs) with combinations of six histone modification marks played crucial roles in the differential expression of salt-responsive genes between seedlings and roots. Most importantly, CS7 containing the bivalent marks H3K4me3 and H3K27me3, with a mutual exclusion of functions with each other, displayed distinct functions in the expression of DEGs in both tissues. Specifically, H3K27me3 in CS7 mainly suppressed the expression of DEGs in roots, while H3K4me3 affected the expression of down- and up-regulated genes, possibly by antagonizing the repressive role of H3K27me3 in seedlings. Our findings indicate distinct impacts of the CSs on the differential expression of salt-responsive genes between seedlings and roots in rice, which provides an important background for understanding chromatin-based epigenetic mechanisms that might confer salt tolerance in plants.

Expression of a Grape VqSTS36-Increased Resistance to Powdery Mildew and Osmotic Stress in Arabidopsis but Enhanced Susceptibility to Botrytis cinerea in Arabidopsis and Tomato

Stilbene synthase genes make a contribution to improving the tolerances of biotic and abiotic stress in plants. However, the mechanisms mediated by these STS genes remain unclear. To provide insight into the role of STS genes defense against biotic and abiotic stress, we overexpressed VqSTS36 in Arabidopsis thaliana and tomato (Micro-Tom) via Agrobacterium-mediated transformation. VqSTS36-transformed Arabidopsis lines displayed an increased resistance to powdery mildew, but both VqSTS36-transformed Arabidopsis and tomato lines showed the increased susceptibility to Botrytis cinerea. Besides, transgenic Arabidopsis lines were found to confer tolerance to salt and drought stress in seed and seedlings. When transgenic plants were treated with a different stress, qPCR assays of defense-related genes in transgenic Arabidopsis and tomato suggested that VqSTS36 played a specific role in different phytohormone-related pathways, including salicylic acid, jasmonic acid, and abscisic acid signaling pathways. All of these results provided a better understanding of the mechanism behind the role of VqSTS36 in biotic and abiotic stress.

Systems biology approach in plant abiotic stresses

Plant abiotic stresses are the major constraint on plant growth and development, causing enormous crop losses across the world. Plants have unique features to defend themselves against these challenging adverse stress conditions. They modulate their phenotypes upon changes in physiological, biochemical, molecular and genetic information, thus making them tolerant against abiotic stresses. It is of paramount importance to determine the stress-tolerant traits of a diverse range of genotypes of plant species and integrate those traits for crop improvement. Stress-tolerant traits can be identified by conducting genome-wide analysis of stress-tolerant genotypes through the highly advanced structural and functional genomics approach. Specifically, whole-genome sequencing, development of molecular markers, genome-wide association studies and comparative analysis of interaction networks between tolerant and susceptible crop varieties grown under stress conditions can greatly facilitate discovery of novel agronomic traits that protect plants against abiotic stresses.

The LATERAL ORGAN BOUNDARIES Domain gene family in grapevine: genome-wide characterization and expression analyses during developmental processes and stress responses

LATERAL ORGAN BOUNDARIES (LOB) DOMAIN (LBD) constitute a family of plant-specific transcription factors with key roles in the regulation of plant organ development, pollen development, plant regeneration, pathogen response, and anthocyanin and nitrogen metabolisms. However, the role of LBDs in fruit ripening and in grapevine (Vitis vinifera L.) development and stress responses is poorly documented. By performing a model curation of LBDs in the latest genome annotation 50 genes were identified. Phylogenetic analysis showed that LBD genes can be grouped into two classes mapping on 16 out of the 19 V. vinifera chromosomes. New gene subclasses were identified that have not been characterized in other species. Segmental and tandem duplications contributed significantly to the expansion and evolution of the LBD gene family in grapevine as noticed for other species. The analysis of cis-regulatory elements and transcription factor binding sites in the VviLBD promoter regions suggests the involvement of several hormones in the regulation of LBDs expression. Expression profiling suggest the involvement of LBD transcription factors in grapevine development, berry ripening and stress responses. Altogether this study provides valuable information and robust candidate genes for future functional analysis aiming to clarify mechanisms responsible for the onset of fruit ripening and fruit defense strategies.

RNA-Sequencing Reveals Biological Networks during Table Grapevine ('Fujiminori') Fruit Development

Grapevine berry development is a complex and genetically controlled process, with many morphological, biochemical and physiological changes occurring during the maturation process. Research carried out on grapevine berry development has been mainly concerned with wine grape, while barely focusing on table grape. 'Fujiminori' is an important table grapevine cultivar, which is cultivated in most provinces of China. In order to uncover the dynamic networks involved in anthocyanin biosynthesis, cell wall development, lipid metabolism and starch-sugar metabolism in 'Fujiminori' fruit, we employed RNA-sequencing (RNA-seq) and analyzed the whole transcriptome of grape berry during development at the expanding period (40 days after full bloom, 40DAF), véraison period (65DAF), and mature period (90DAF). The sequencing depth in each sample was greater than 12×, and the expression level of nearly half of the expressed genes were greater than 1. Moreover, greater than 64% of the clean reads were aligned to the Vitis vinifera reference genome, and 5,620, 3,381, and 5,196 differentially expressed genes (DEGs) were identified between different fruit stages, respectively. Results of the analysis of DEGs showed that the most significant changes in various processes occurred from the expanding stage to the véraison stage. The expression patterns of F3'H and F3'5'H were crucial in determining red or blue color of the fruit skin. The dynamic networks of cell wall development, lipid metabolism and starch-sugar metabolism were also constructed. A total of 4,934 SSR loci were also identified from 4,337 grapevine genes, which may be helpful for the development of phylogenetic analysis in grapevine and other fruit trees. Our work provides the foundation for developmental research of grapevine fruit as well as other non-climacteric fruits.

Current Knowledge and Computational Techniques for Grapevine Meta-Omics Analysis

Growing grapevine (Vitis vinifera) is a key contribution to the economy of many countries. Tools provided by genomics and bioinformatics did help researchers in obtaining biological knowledge about the different cultivars. Several genetic markers for common diseases were identified. Recently, the impact of microbiome has been proved to be of fundamental importance both in humans and in plants for its ability to confer protection or induce diseases. In this review we report current knowledge about grapevine microbiome, together with a description of the available computational methodologies for meta-omics analysis.

Structural and Functional Analysis of the GRAS Gene Family in Grapevine Indicates a Role of GRAS Proteins in the Control of Development and Stress Responses

GRAS transcription factors are involved in many processes of plant growth and development (e.g., axillary shoot meristem formation, root radial patterning, nodule morphogenesis, arbuscular development) as well as in plant disease resistance and abiotic stress responses. However, little information is available concerning this gene family in grapevine (Vitis vinifera L.), an economically important woody crop. We performed a model curation of GRAS genes identified in the latest genome annotation leading to the identification of 52 genes. Gene models were improved and three new genes were identified that could be grapevine- or woody-plant specific. Phylogenetic analysis showed that GRAS genes could be classified into 13 groups that mapped on the 19 V. vinifera chromosomes. Five new subfamilies, previously not characterized in other species, were identified. Multiple sequence alignment showed typical GRAS domain in the proteins and new motifs were also described. As observed in other species, both segmental and tandem duplications contributed significantly to the expansion and evolution of the GRAS gene family in grapevine. Expression patterns across a variety of tissues and upon abiotic and biotic conditions revealed possible divergent functions of GRAS genes in grapevine development and stress responses. By comparing the information available for tomato and grapevine GRAS genes, we identified candidate genes that might constitute conserved transcriptional regulators of both climacteric and non-climacteric fruit ripening. Altogether this study provides valuable information and robust candidate genes for future functional analysis aiming at improving the quality of fleshy fruits.

Comprehensive transcript profiling of two grapevine rootstock genotypes contrasting in drought susceptibility links the phenylpropanoid pathway to enhanced tolerance

In light of ongoing climate changes in wine-growing regions, the selection of drought-tolerant rootstocks is becoming a crucial factor for developing a sustainable viticulture. In this study, M4, a new rootstock genotype that shows tolerance to drought, was compared from a genomic and transcriptomic point of view with the less drought-tolerant genotype 101.14. The root and leaf transcriptome of both 101.14 and the M4 rootstock genotype was analysed, following exposure to progressive drought conditions. Multifactorial analyses indicated that stress treatment represents the main factor driving differential gene expression in roots, whereas in leaves the genotype is the prominent factor. Upon stress, M4 roots and leaves showed a higher induction of resveratrol and flavonoid biosynthetic genes, respectively. The higher expression of VvSTS genes in M4, confirmed by the accumulation of higher levels of resveratrol in M4 roots compared with 101.14, was coupled to an up-regulation of several VvWRKY transcription factors. Interestingly, VvSTS promoter analyses performed on both the resequenced genomes highlighted a significantly higher number of W-BOX elements in the tolerant genotype. It is proposed that the elevated synthesis of resveratrol in M4 roots upon water stress could enhance the plant's ability to cope with the oxidative stress usually associated with water deficit.

A knowledge base for Vitis vinifera functional analysis

Background

Vitis vinifera (Grapevine) is the most important fruit species in the modern world. Wine and table grapes sales contribute significantly to the economy of major wine producing countries. The most relevant goals in wine production concern quality and safety. In order to significantly improve the achievement of these objectives and to gain biological knowledge about cultivars, a genomic approach is the most reliable strategy. The recent grapevine genome sequencing offers the opportunity to study the potential roles of genes and microRNAs in fruit maturation and other physiological and pathological processes. Although several systems allowing the analysis of plant genomes have been reported, none of them has been designed specifically for the functional analysis of grapevine genomes of cultivars under environmental stress in connection with microRNA data.

Description

Here we introduce a novel knowledge base, called BIOWINE, designed for the functional analysis of Vitis vinifera genomes of cultivars present in Sicily. The system allows the analysis of RNA-seq experiments of two different cultivars, namely Nero d'Avola and Nerello Mascalese. Samples were taken under different climatic conditions of phenological phases, diseases, and geographic locations. The BIOWINE web interface is equipped with data analysis modules for grapevine genomes. In particular users may analyze the current genome assembly together with the RNA-seq data through a customized version of GBrowse. The web interface allows users to perform gene set enrichment by exploiting third-party databases.

Conclusions

BIOWINE is a knowledge base implementing a set of bioinformatics tools for the analysis of grapevine genomes. The system aims to increase our understanding of the grapevine varieties and species of Sicilian products focusing on adaptability to different climatic conditions, phenological phases, diseases, and geographic locations.

The elucidation of stress memory inheritance in Brassica rapa plants

Plants are able to maintain the memory of stress exposure throughout their ontogenesis and faithfully propagate it into the next generation. Recent evidence argues for the epigenetic nature of this phenomenon. Small RNAs (smRNAs) are one of the vital epigenetic factors because they can both affect gene expression at the place of their generation and maintain non-cell-autonomous gene regulation. Here, we have made an attempt to decipher the contribution of smRNAs to the heat-shock-induced transgenerational inheritance in Brassica rapa plants using sequencing technology. To do this, we have generated comprehensive profiles of a transcriptome and a small RNAome (smRNAome) from somatic and reproductive tissues of stressed plants and their untreated progeny. We have demonstrated that the highest tissue-specific alterations in the transcriptome and smRNAome profile are detected in tissues that were not directly exposed to stress, namely, in the endosperm and pollen. Importantly, we have revealed that the progeny of stressed plants exhibit the highest fluctuations at the smRNAome level but not at the transcriptome level. Additionally, we have uncovered the existence of heat-inducible and transgenerationally transmitted tRNA-derived small RNA fragments in plants. Finally, we suggest that miR168 and braAGO1 are involved in the stress-induced transgenerational inheritance in plants.

Shoot chloride exclusion and salt tolerance in grapevine is associated with differential ion transporter expression in roots

BACKGROUND

Salt tolerance in grapevine is associated with chloride (Cl-) exclusion from shoots; the rate-limiting step being the passage of Cl- between the root symplast and xylem apoplast. Despite an understanding of the physiological mechanism of Cl- exclusion in grapevine, the molecular identity of membrane proteins that control this process have remained elusive. To elucidate candidate genes likely to control Cl- exclusion, we compared the root transcriptomes of three Vitis spp. with contrasting shoot Cl- exclusion capacities using a custom microarray.

RESULTS

When challenged with 50 mM Cl-, transcriptional changes of genotypes 140 Ruggeri (shoot Cl- excluding rootstock), K51-40 (shoot Cl- including rootstock) and Cabernet Sauvignon (intermediate shoot Cl- excluder) differed. The magnitude of salt-induced transcriptional changes in roots correlated with the amount of Cl- accumulated in shoots. Abiotic-stress responsive transcripts (e.g. heat shock proteins) were induced in 140 Ruggeri, respiratory transcripts were repressed in Cabernet Sauvignon, and the expression of hypersensitive response and ROS scavenging transcripts was altered in K51-40. Despite these differences, no obvious Cl- transporters were identified. However, under control conditions where differences in shoot Cl- exclusion between rootstocks were still significant, genes encoding putative ion channels SLAH3, ALMT1 and putative kinases SnRK2.6 and CPKs were differentially expressed between rootstocks, as were members of the NRT1 (NAXT1 and NRT1.4), and CLC families.

CONCLUSIONS

These results suggest that transcriptional events contributing to the Cl- exclusion mechanism in grapevine are not stress-inducible, but constitutively different between contrasting varieties. We have identified individual genes from large families known to have members with roles in anion transport in other plants, as likely candidates for controlling anion homeostasis and Cl- exclusion in Vitis species. We propose these genes as priority candidates for functional characterisation to determine their role in chloride transport in grapevine and other plants.

Microarray analysis of differentially expressed genes engaged in fruit development between table and wine grape

Microarray analysis of genes can provide individual gene-expression profiles and new insights for elucidating biological mechanisms responsible for fruit development. To obtain an overall view on expression profiles of metabolism-related genes involved in fruit development of table and wine grapes, a microarray system comprising 15,403 ESTs was used to compare the expressed genes. The expression patterns from the microarray analysis were validated with quantitative real-time polymerase chain reaction analysis of 18 selected genes of interest. During the entire fruit development stage, 2,493 genes exhibited at least 2.0-fold differences in expression levels with 1,244 genes being up-regulated and 1,249 being down-regulated. Following gene ontology analysis, only 929 differentially expressed (including 403 up-regulated and 526 down-regulated) genes were annotated in table and wine grapes. These differentially expressed genes were found to be mainly involved in carbohydrate metabolism, biosynthesis of secondary metabolites as well as energy, lipid and amino acid metabolism via KEGG. Our results provide new insights into the molecular mechanisms and expression profiles of genes in the fruit development stage of table and wine grapes.

Nitrogen starvation, salt and heat stress in coffee (Coffea arabica L.): identification and validation of new genes for qPCR normalization

Abiotic stresses are among the most important factors that affect food production. One important step to face these environmental challenges is the transcriptional modulation. Quantitative real-time PCR is a rapid, sensitive, and reliable method for the detection of mRNAs and it has become a powerful tool to mitigate plant stress tolerance; however, suitable reference genes are required for data normalization. Reference genes for coffee plants during nitrogen starvation, salinity and heat stress have not yet been reported. We evaluated the expression stability of ten candidate reference genes using geNorm PLUS, NormFinder, and BestKeeper softwares, in plants submitted to nitrogen starvation, salt and heat stress. EF1, EF1α, GAPDH, MDH, and UBQ10 were ranked as the most stable genes in all stresses and software analyses, while RPL39 and RPII were classified as the less reliable references. For reference gene validation, the transcriptional pattern of a Coffea non-symbiotic hemoglobin (CaHb1) was analyzed using the two new recommended and the most unstable gene references for normalization. The most unstable gene may lead to incorrect interpretation of CaHb1 transcriptional analysis. Here, we recommend two new reference genes in Coffea for use in data normalization in abiotic stresses: MDH and EF1.

Salt-dependent increase in triterpenoids is reversible upon transfer to fresh water in mangrove plants Kandelia candel and Bruguiera gymnorrhiza

This study examined the salinity dependence of triterpenoid content and triterpenoid synthase gene expression in mangrove plants, Kandelia candel and Bruguiera gymnorrhiza (Rhizophoraceae) after long-term exposure to salinity and subsequent re-adaptation. Seedlings of the two mangrove species grown in varying salt concentrations for 4 months were divided into two treatment groups and grown for another 4 months, one group continued under the respective saline condition and the other in fresh water for re-adaptation. The total content of triterpenoids increased with increasing salinity in roots and leaves of K. candel, but only in roots in B. gymnorrhiza. This increase was reversed to a variable extent, depending on the species and organ, after transfer to fresh water. In contrast, the total content of phytosterols showed no correlation with salinity throughout the experiment. The increase in total triterpenoids was accompanied by an up-regulation of several triterpenoid synthase genes: KcMS, a multifunctional triterpenoid synthase, in roots and leaves of K. candel and BgLUS, a lupeol synthase, and BgbAS, a β-amyrin synthase, in roots of B. gymnorrhiza. The expression of root KcCAS, a cycloartenol synthase, which is involved in phytosterol biosynthesis, was not modulated by the salinity conditions but decreased with increasing salinity in leaves, followed by the restoration to the initial level after transfer to fresh water. The concentrations of individual triterpenoids, but not of phytosterols, in the roots positively correlated with the salinity. These results reinforced the importance of triterpenoids in the adaptation of mangroves to withstand salt and/or water stress.

Meta-analytical biomarker search of EST expression data reveals three differentially expressed candidates

Background

Researches have been conducted for the identification of differentially expressed genes (DEGs) by generating and mining of cDNA expressed sequence tags (ESTs) for more than a decade. Although the availability of public databases make possible the comprehensive mining of DEGs among the ESTs from multiple tissue types, existing studies usually employed statistics suitable only for two categories. Multi-class test has been developed to enable the finding of tissue specific genes, but subsequent search for cancer genes involves separate two-category test only on the ESTs of the tissue of interest. This constricts the amount of data used. On the other hand, simple pooling of cancer and normal genes from multiple tissue types runs the risk of Simpson's paradox. Here we presented a different approach which searched for multi-cancer DEG candidates by analyzing all pertinent ESTs in all categories and narrowing down the cancer biomarker candidates via integrative analysis with microarray data and selection of secretory and membrane protein genes as well as incorporation of network analysis. Finally, the differential expression patterns of three selected cancer biomarker candidates were confirmed by real-time qPCR analysis.

Results

Seven hundred and twenty three primary DEG candidates (p-value < 0.05 and lower bound of confidence interval of odds ratio ≧ 1.65) were selected from a curated EST database with the application of Cochran-Mantel-Haenszel statistic (CMH). GeneGO analysis results indicated this set as neoplasm enriched. Cross-examination with microarray data further narrowed the list down to 235 genes, among which 96 had membrane or secretory annotations. After examined the candidates in protein interaction network, public tissue expression databases, and literatures, we selected three genes for further evaluation by real-time qPCR with eight major normal and cancer tissues. The higher-than-normal tissue expression of COL3A1, DLG3, and RNF43 in some of the cancer tissues is in agreement with our in silico predictions.

Conclusions

Searching digitized transcriptome using CMH enabled us to identify multi-cancer differentially expressed gene candidates. Our methodology demonstrated simultaneously analysis for cancer biomarkers of multiple tissue types with the EST data. With the revived interest in digitizing the transcriptomes by NGS, cancer biomarkers could be more precisely detected from the ESTs. The three candidates identified in this study, COL3A1, DLG3, and RNF43, are valuable targets for further evaluation with a larger sample size of normal and cancer tissue or serum samples.

Small RNA and transcriptome deep sequencing proffers insight into floral gene regulation in Rosa cultivars

Background

Roses (Rosa sp.), which belong to the family Rosaceae, are the most economically important ornamental plants—making up 30% of the floriculture market. However, given high demand for roses, rose breeding programs are limited in molecular resources which can greatly enhance and speed breeding efforts. A better understanding of important genes that contribute to important floral development and desired phenotypes will lead to improved rose cultivars. For this study, we analyzed rose miRNAs and the rose flower transcriptome in order to generate a database to expound upon current knowledge regarding regulation of important floral characteristics. A rose genetic database will enable comprehensive analysis of gene expression and regulation via miRNA among different Rosa cultivars.

Results

We produced more than 0.5 million reads from expressed sequences, totalling more than 110 million bp. From these, we generated 35,657, 31,434, 34,725, and 39,722 flower unigenes from Rosa hybrid: ‘Vital’, ‘Maroussia’, and ‘Sympathy’ and Rosa rugosa Thunb. , respectively. The unigenes were assigned functional annotations, domains, metabolic pathways, Gene Ontology (GO) terms, Plant Ontology (PO) terms, and MIPS Functional Catalogue (FunCat) terms. Rose flower transcripts were compared with genes from whole genome sequences of Rosaceae members (apple, strawberry, and peach) and grape. We also produced approximately 40 million small RNA reads from flower tissue for Rosa, representing 267 unique miRNA tags. Among identified miRNAs, 25 of them were novel and 242 of them were conserved miRNAs. Statistical analyses of miRNA profiles revealed both shared and species-specific miRNAs, which presumably effect flower development and phenotypes.

Conclusions

In this study, we constructed a Rose miRNA and transcriptome database, and we analyzed the miRNAs and transcriptome generated from the flower tissues of four Rosa cultivars. The database provides a comprehensive genetic resource which can be used to better understand rose flower development and to identify candidate genes for important phenotypes.

ST proteins, a new family of plant tandem repeat proteins with a DUF2775 domain mainly found in Fabaceae and Asteraceae

BACKGROUND

Many proteins with tandem repeats in their sequence have been described and classified according to the length of the repeats: I) Repeats of short oligopeptides (from 2 to 20 amino acids), including structural cell wall proteins and arabinogalactan proteins. II) Repeats that range in length from 20 to 40 residues, including proteins with a well-established three-dimensional structure often involved in mediating protein-protein interactions. (III) Longer repeats in the order of 100 amino acids that constitute structurally and functionally independent units. Here we analyse ShooT specific (ST) proteins, a family of proteins with tandem repeats of unknown function that were first found in Leguminosae, and their possible similarities to other proteins with tandem repeats.

RESULTS

ST protein sequences were only found in dicotyledonous plants, limited to several plant families, mainly the Fabaceae and the Asteraceae. ST mRNAs accumulate mainly in the roots and under biotic interactions. Most ST proteins have one or several Domain(s) of Unknown Function 2775 (DUF2775). All deduced ST proteins have a signal peptide, indicating that these proteins enter the secretory pathway, and the mature proteins have tandem repeat oligopeptides that share a hexapeptide (E/D)FEPRP followed by 4 partially conserved amino acids, which could determine a putative N-glycosylation signal, and a fully conserved tyrosine. In a phylogenetic tree, the sequences clade according to taxonomic group. A possible involvement in symbiosis and abiotic stress as well as in plant cell elongation is suggested, although different STs could play different roles in plant development.

CONCLUSIONS

We describe a new family of proteins called ST whose presence is limited to the plant kingdom, specifically to a few families of dicotyledonous plants. They present 20 to 40 amino acid tandem repeat sequences with different characteristics (signal peptide, DUF2775 domain, conservative repeat regions) from the described group of 20 to 40 amino acid tandem repeat proteins and also from known cell wall proteins with repeat sequences. Several putative roles in plant physiology can be inferred from the characteristics found.

Yellow lupin (Lupinus luteus L.) transcriptome sequencing: molecular marker development and comparative studies

Background

Yellow lupin (Lupinus luteus L.) is a minor legume crop characterized by its high seed protein content. Although grown in several temperate countries, its orphan condition has limited the generation of genomic tools to aid breeding efforts to improve yield and nutritional quality. In this study, we report the construction of 454-expresed sequence tag (EST) libraries, carried out comparative studies between L. luteus and model legume species, developed a comprehensive set of EST-simple sequence repeat (SSR) markers, and validated their utility on diversity studies and transferability to related species.

Results

Two runs of 454 pyrosequencing yielded 205 Mb and 530 Mb of sequence data for L1 (young leaves, buds and flowers) and L2 (immature seeds) EST- libraries. A combined assembly (L1L2) yielded 71,655 contigs with an average contig length of 632 nucleotides. L1L2 contigs were clustered into 55,309 isotigs. 38,200 isotigs translated into proteins and 8,741 of them were full length. Around 57% of L. luteus sequences had significant similarity with at least one sequence of Medicago, Lotus, Arabidopsis, or Glycine, and 40.17% showed positive matches with all of these species. L. luteus isotigs were also screened for the presence of SSR sequences. A total of 2,572 isotigs contained at least one EST-SSR, with a frequency of one SSR per 17.75 kbp. Empirical evaluation of the EST-SSR candidate markers resulted in 222 polymorphic EST-SSRs. Two hundred and fifty four (65.7%) and 113 (30%) SSR primer pairs were able to amplify fragments from L. hispanicus and L. mutabilis DNA, respectively. Fifty polymorphic EST-SSRs were used to genotype a sample of 64 L. luteus accessions. Neighbor-joining distance analysis detected the existence of several clusters among L. luteus accessions, strongly suggesting the existence of population subdivisions. However, no clear clustering patterns followed the accession’s origin.

Conclusion

L. luteus deep transcriptome sequencing will facilitate the further development of genomic tools and lupin germplasm. Massive sequencing of cDNA libraries will continue to produce raw materials for gene discovery, identification of polymorphisms (SNPs, EST-SSRs, INDELs, etc.) for marker development, anchoring sequences for genome comparisons and putative gene candidates for QTL detection.

Differences in gene expression between strabismic and normal human extraocular muscles

PURPOSE

Strabismic extraocular muscles (EOMs) differ from normal EOMs in structural and functional properties, but the gene expression profile of these two types of EOM has not been examined. Differences in gene expression may inform about causes and effects of the strabismic condition in humans.

METHODS

EOM samples were obtained during corrective surgery from patients with horizontal strabismus and from deceased organ donors with normal EOMs. Microarrays and quantitative PCR identified significantly up- and down-regulated genes in EOM samples. Analysis was performed on probe sets with more than 3-fold differential expression between normal and strabismic samples, with an adjusted P value of ≤ 0.05.

RESULTS

Microarray analysis showed that 604 genes in these samples had significantly different expression. Expression predominantly was upregulated in genes involved in extracellular matrix structure, and down-regulated in genes related to contractility. Expression of genes associated with signaling, calcium handling, mitochondria function and biogenesis, and energy homeostasis also was significantly different between normal and strabismic EOM. Skeletal muscle PCR array identified 22 (25%) of 87 muscle-specific genes that were significantly down-regulated in strabismic EOMs; none was significantly upregulated.

CONCLUSIONS

Differences in gene expression between strabismic and normal human EOMs point to a relevant contribution of the peripheral oculomotor system to the strabismic condition. Decreases in expression of contractility genes and increases of extracellular matrix-associated genes indicate imbalances in EOM structure. We conclude that gene regulation of proteins fundamental to contractile mechanics and extracellular matrix structure is involved in pathogenesis and/or consequences of strabismus, suggesting potential novel therapeutic targets.

SNP markers retrieval for a non-model species: a practical approach

Background

SNP (Single Nucleotide Polymorphism) markers are rapidly becoming the markers of choice for applications in breeding because of next generation sequencing technology developments. For SNP development by NGS technologies, correct assembly of the huge amounts of sequence data generated is essential. Little is known about assembler's performance, especially when dealing with highly heterogeneous species that show a high genome complexity and what the possible consequences are of differences in assemblies on SNP retrieval. This study tested two assemblers (CAP3 and CLC) on 454 data from four lily genotypes and compared results with respect to SNP retrieval.

Results

CAP3 assembly resulted in higher numbers of contigs, lower numbers of reads per contig, and shorter average read lengths compared to CLC. Blast comparisons showed that CAP3 contigs were highly redundant. Contrastingly, CLC in rare cases combined paralogs in one contig. Redundant and chimeric contigs may lead to erroneous SNPs. Filtering for redundancy can be done by blasting selected SNP markers to the contigs and discarding all the SNP markers that show more than one blast hit. Results on chimeric contigs showed that only four out of 2,421 SNP markers were selected from chimeric contigs.

Conclusion

In practice, CLC performs better in assembling highly heterogeneous genome sequences compared to CAP3, and consequently SNP retrieval is more efficient. Additionally a simple flow scheme is suggested for SNP marker retrieval that can be valid for all non-model species.