DEAD-Box RNA Helicase Family in Physic Nut (Jatropha curcas L.): Structural Characterization and Response to Salinity.

Helicases, motor proteins present in both prokaryotes and eukaryotes, play a direct role in various steps of RNA metabolism. Specifically, SF2 RNA helicases, a subset of the DEAD-box family, are essential players in plant developmental processes and responses to biotic and abiotic stresses. Despite this, information on this family in the physic nut (Jatropha curcas L.) remains limited, spanning from structural patterns to stress responses. We identified 79 genes encoding DEAD-box RNA helicases (JcDHX) in the J. curcas genome. These genes were further categorized into three subfamilies: DEAD (42 genes), DEAH (30 genes), and DExH/D (seven genes). Characterization of the encoded proteins revealed a remarkable diversity, with observed patterns in domains, motifs, and exon-intron structures suggesting that the DEAH and DExH/D subfamilies in J. curcas likely contribute to the overall versatility of the family. Three-dimensional modeling of the candidates showed characteristic hallmarks, highlighting the expected functional performance of these enzymes. The promoter regions of the JcDHX genes revealed potential cis-elements such as Dof-type, BBR-BPC, and AP2-ERF, indicating their potential involvement in the response to abiotic stresses. Analysis of RNA-Seq data from the roots of physic nut accessions exposed to 150 mM of NaCl for 3 h showed most of the JcDHX candidates repressed. The protein-protein interaction network indicated that JcDHX proteins occupy central positions, connecting events associated with RNA metabolism. Quantitative PCR analysis validated the expression of nine DEAD-box RNA helicase transcripts, showing significant associations with key components of the stress response, including RNA turnover, ribosome biogenesis, DNA repair, clathrin-mediated vesicular transport, phosphatidyl 3,5-inositol synthesis, and mitochondrial translation. Furthermore, the induced expression of one transcript (JcDHX44) was confirmed, suggesting that it is a potential candidate for future functional analyses to better understand its role in salinity stress tolerance. This study represents the first global report on the DEAD-box family of RNA helicases in physic nuts and displays structural characteristics compatible with their functions, likely serving as a critical component of the plant's response pathways.

Exploiting DNA methylation in cassava under water deficit for crop improvement.

DNA methylation plays a key role in the development and plant responses to biotic and abiotic stresses. This work aimed to evaluate the DNA methylation in contrasting cassava genotypes for water deficit tolerance. The varieties BRS Formosa (bitter) and BRS Dourada (sweet) were grown under greenhouse conditions for 50 days, and afterwards, irrigation was suspended. The stressed (water deficit) and non-stressed plants (negative control) consisted the treatments with five plants per variety. The DNA samples of each variety and treatment provided 12 MethylRAD-Seq libraries (two cassava varieties, two treatments, and three replicates). The sequenced data revealed methylated sites covering 18 to 21% of the Manihot esculenta Crantz genome, depending on the variety and the treatment. The CCGG methylated sites mapped mostly in intergenic regions, exons, and introns, while the CCNGG sites mapped mostly intergenic, upstream, introns, and exons regions. In both cases, methylated sites in UTRs were less detected. The differentially methylated sites analysis indicated distinct methylation profiles since only 12% of the sites (CCGG and CCNGG) were methylated in both varieties. Enriched gene ontology terms highlighted the immediate response of the bitter variety to stress, while the sweet variety appears to suffer more potential stress-damages. The predicted protein-protein interaction networks reinforced such profiles. Additionally, the genomes of the BRS varieties uncovered SNPs/INDELs events covering genes stood out by the interactomes. Our data can be useful in deciphering the roles of DNA methylation in cassava drought-tolerance responses and adaptation to abiotic stresses.

The WRKY transcription factor family in cowpea: Genomic characterization and transcriptomic profiling under root dehydration.

Cowpea [Vigna unguiculata (L.) Walp.] is one of the most tolerant legume crops to drought and salt stresses. WRKY transcription factor (TF) family members stand out among plant transcriptional regulators related to abiotic stress tolerance. However, little information is currently available on the expression of the cowpea WRKY gene family (VuWRKY) in response to water deficit. Thus, we analyzed genomic and transcriptomic data from cowpea to identify VuWRKY members and characterize their structure and transcriptional response under root dehydration stress. Ninety-two complete VuWRKY genes were found in the cowpea genome based on their domain characteristics. They were clustered into three groups: I (15 members), II (58), and III (16), while three genes were unclassified. Domain analysis of the encoded proteins identified four major variants of the conserved heptapeptide motif WRKYGQK. In silico analysis of VuWRKY gene promoters identified eight candidate binding motifs of cis-regulatory elements, regulated mainly by six TF families associated with abiotic stress responses. Ninety-seven VuWRKY modulated splicing variants associated with 55 VuWRKY genes were identified via RNA-Seq analysis available at the Cowpea Genomics Consortium (CpGC) database. qPCR analyses showed that 22 genes are induced under root dehydration, with VuWRKY18, 21, and 75 exhibiting the most significant induction levels. Given their central role in activating signal transduction cascades in abiotic stress response, the data provide a foundation for the targeted modification of specific VuWRKY family members to improve drought tolerance in this important climate-resilient legume in the developing world and beyond.

The Cowpea Kinome: Genomic and Transcriptomic Analysis Under Biotic and Abiotic Stresses.

The present work represents a pioneering effort, being the first to analyze genomic and transcriptomic data from Vigna unguiculata (cowpea) kinases. We evaluated the cowpea kinome considering its genome-wide distribution and structural characteristics (at the gene and protein levels), sequence evolution, conservation among Viridiplantae species, and gene expression in three cowpea genotypes under different stress situations, including biotic (injury followed by virus inoculation-CABMV or CPSMV) and abiotic (root dehydration). The structural features of cowpea kinases (VuPKs) indicated that 1,293 bona fide VuPKs covered 20 groups and 118 different families. The RLK-Pelle was the largest group, with 908 members. Insights on the mechanisms of VuPK genomic expansion and conservation among Viridiplantae species indicated dispersed and tandem duplications as major forces for VuPKs' distribution pattern and high orthology indexes and synteny with other legume species, respectively. K a /K s ratios showed that almost all (91%) of the tandem duplication events were under purifying selection. Candidate cis-regulatory elements were associated with different transcription factors (TFs) in the promoter regions of the RLK-Pelle group. C2H2 TFs were closely associated with the promoter regions of almost all scrutinized families for the mentioned group. At the transcriptional level, it was suggested that VuPK up-regulation was stress, genotype, or tissue dependent (or a combination of them). The most prominent families in responding (up-regulation) to all the analyzed stresses were RLK-Pelle_DLSV and CAMK_CAMKL-CHK1. Concerning root dehydration, it was suggested that the up-regulated VuPKs are associated with ABA hormone signaling, auxin hormone transport, and potassium ion metabolism. Additionally, up-regulated VuPKs under root dehydration potentially assist in a critical physiological strategy of the studied cowpea genotype in this assay, with activation of defense mechanisms against biotic stress while responding to root dehydration. This study provides the foundation for further studies on the evolution and molecular function of VuPKs.

Reference genes for quantitative real-time PCR normalization of Cenostigma pyramidale roots under salt stress and mycorrhizal association.

Cenostigma pyramidale is a native legume of the Brazilian semiarid region which performs symbiotic association with arbuscular mycorrhizal fungi (AMF), being an excellent model for studying genes associated with tolerance against abiotic and biotic stresses. In RT-qPCR approach, the use of reference genes is mandatory to avoid incorrect interpretation of the relative expression. This study evaluated the stability of ten candidate reference genes (CRGs) from C. pyramidale root tissues under salt stress (three collection times) and associated with AMF (three different times of salinity). The de novo transcriptome was obtained via RNA-Seq sequencing. Three algorithms were used to calculate the stability of CRGs under different conditions: (i) global (Salt, Salt+AMF, AMF and Control, and collection times), (ii) only non-inoculated plants, and (iii) AMF (only inoculated plants). HAG2, SAC1, aRP3 were the most stable CRGs for global and AMF assays, whereas HAG2, SAC1, RHS1 were the best for salt stress assay. This CRGs were used to validate the relative expression of two up-regulated transcripts in Salt2h (RAP2-3 and PIN8). Our study provides the first set of reference genes for C. pyramidale under salinity and AMF, supporting future researches on gene expression with this species.

Plant Antimicrobial Peptides: State of the Art, In Silico Prediction and Perspectives in the Omics Era.

Even before the perception or interaction with pathogens, plants rely on constitutively guardian molecules, often specific to tissue or stage, with further expression after contact with the pathogen. These guardians include small molecules as antimicrobial peptides (AMPs), generally cysteine-rich, functioning to prevent pathogen establishment. Some of these AMPs are shared among eukaryotes (eg, defensins and cyclotides), others are plant specific (eg, snakins), while some are specific to certain plant families (such as heveins). When compared with other organisms, plants tend to present a higher amount of AMP isoforms due to gene duplications or polyploidy, an occurrence possibly also associated with the sessile habit of plants, which prevents them from evading biotic and environmental stresses. Therefore, plants arise as a rich resource for new AMPs. As these molecules are difficult to retrieve from databases using simple sequence alignments, a description of their characteristics and in silico (bioinformatics) approaches used to retrieve them is provided, considering resources and databases available. The possibilities and applications based on tools versus database approaches are considerable and have been so far underestimated.

Plant Thaumatin-like Proteins: Function, Evolution and Biotechnological Applications.

Thaumatin-like proteins (TLPs) are a highly complex protein family associated with host defense and developmental processes in plants, animals, and fungi. They are highly diverse in angiosperms, for which they are classified as the PR-5 (Pathogenesis-Related-5) protein family. In plants, TLPs have a variety of properties associated with their structural diversity. They are mostly associated with responses to biotic stresses, in addition to some predicted activities under drought and osmotic stresses. The present review covers aspects related to the structure, evolution, gene expression, and biotechnological potential of TLPs. The efficiency of the discovery of new TLPs is below its potential, considering the availability of omics data. Furthermore, we present an exemplary bioinformatics annotation procedure that was applied to cowpea (Vigna unguiculata) transcriptome, including libraries of two tissues (root and leaf), and two stress types (biotic/abiotic) generated using different sequencing approaches. Even without using genomic sequences, the pipeline uncovered 56 TLP candidates in both tissues and stresses. Interestingly, abiotic stress (root dehydration) was associated with a high number of modulated TLP isoforms. The nomenclature used so far for TLPs was also evaluated, considering TLP structure and possible functions identified to date. It is clear that plant TLPs are promising candidates for breeding purposes and for plant transformation aiming a better performance under biotic and abiotic stresses. The development of new therapeutic drugs against human fungal pathogens also deserves attention. Despite that, applications derived from TLP molecules are still below their potential, as it is evident in our review.

Cowpea and abiotic stresses: identification of reference genes for transcriptional profiling by qPCR.

BACKGROUND: Due to cowpea ability to fix nitrogen in poor soils and relative tolerance to drought and salt stresses, efforts have been directed to identifying genes and pathways that confer stress tolerance in this species. Real-time quantitative PCR (qPCR) has been widely used as the most reliable method to measure gene expression, due to its high accuracy and specificity. In the present study, nine candidate reference genes were rigorously tested for their application in normalization of qPCR data onto roots of four distinct cowpea accessions under two abiotic stresses: root dehydration and salt (NaCl, 100 mM). In addition, the regulation of four target transcripts, under the same referred conditions was also scrutinized. RESULTS: geNorm, NormFinder, BestKeeper, and ΔCt method results indicated a set of three statistically validated RGs for each stress condition: (I) root dehydration (actin, ubiquitin-conjugating enzyme E2 variant 1D, and a Phaseolus vulgaris unknown gene-UNK), and (II) salt (ubiquitin-conjugating enzyme E2 variant 1D, F-box protein, and UNK). The expression profile of the target transcripts suggests that flavonoids are important players in the cowpea response to the abiotic stresses analyzed, since chalcone isomerase and chalcone synthase were up-regulated in the tolerant and sensitive accessions. A lipid transfer protein also participates in the cowpea tolerance mechanisms to root dehydration and salt stress. The referred transcript was up-regulated in the two tolerant accessions and presented no differential expression in the sensitive counterparts. Chitinase B, in turn, generally related to plant defense, was an important target transcript under salt stress, being up-regulated at the tolerant, and down-regulated in the sensitive accession. CONCLUSIONS: Reference genes suitable for qPCR analyses in cowpea under root dehydration and salt stress were identified. This action will lead to a more accurate and reliable analysis of gene expression on this species. Additionally, the results obtained in this study may guide future research on gene expression in cowpea under other abiotic stress types that impose osmotic imbalance. The target genes analyzed, in turn, deserve functional evaluation due to their transcriptional regulation under stresses and biotechnological potential.

Genotype-dependent regulation of drought-responsive genes in tolerant and sensitive sugarcane cultivars.

Drought is the most damaging among the major abiotic stresses. Transcriptomic studies allow a global overview of expressed genes, providing the basis for molecular markers development. Here, the HT-SuperSAGE technique allowed the evaluation of four drought-tolerant cultivars and four-sensitive cultivars, after 24h of irrigation suppression. We identified 9831 induced unitags from roots of the tolerant cultivars with different regulations by the -sensitive cultivars after the applied stress. These unitags allowed a proposal of 15 genes, whose expressed profiles were validated by RT-qPCR, evaluating each cultivar independently. These genes covered broad metabolic processes: ethylene stress attenuation (ACCD); root growth (ß-EXP8); protein degradation [ubiquitination pathway (E2, 20SPß4); plant proteases (AP, C13)]; oxidative detoxification (TRX); fatty acid synthesis (ACC); amino acid transport (AAT), and carbohydrate metabolism [glycolysis (PFK, TPI, FBA); TCA cycle (LDP, MDH); pentose phosphate pathway (TKT)]. The expressed profiles showed a genotype-dependent regulation of the target genes. Two drought-tolerant cultivars (SP83-2847; CTC6) presented each one, nine of the induced genes. Among the -sensitive cultivars, CTC13 induced only one, while SP90-1636 induced two genes. These genes should help breeders to identify accessions managing drought stress tolerance responses, showing better ethylene stress attenuation, energy allocation, amino acid transport, and protein homeostasis.

Snakin: Structure, Roles and Applications of a Plant Antimicrobial Peptide.

Snakins are plant antimicrobial peptides (AMPs) of the Snakin/GASA family, formed by three distinct regions: an N-terminal signal peptide; a variable site; and the GASA domain in the Cterminal region composed by twelve conserved cysteine residues that contribute to the biochemical stability of the molecule. These peptides are known to play different roles in response to a variety of biotic (i.e., induced by bacteria, fungi and nematode pathogens) and abiotic (salinity, drought and ROS) stressors, as well as in crosstalk promoted by plant hormones, with emphasis on abscisic and salicylic acid (ABA and SA, respectively). Such properties make snakin/GASA members promising biotechnological sources for potential therapeutic and agricultural applications. However, information regarding their tertiary structure, mode of action and function are not yet completely elucidated. The present review presents aspects of snakin structure, expression, functional studies and perspectives about the potential applications for agricultural and medical purposes.

Resistance (R) Genes: Applications and Prospects for Plant Biotechnology and Breeding.

The discovery of novel plant resistance (R) genes (including their homologs and analogs) opened interesting possibilities for controlling plant diseases caused by several pathogens. However, due to environmental pressure and high selection operated by pathogens, several crop plants have lost specificity, broad-spectrum or durability of resistance. On the other hand, the advances in plant genome sequencing and biotechnological approaches, combined with the increasing knowledge on Rgenes have provided new insights on their applications for plant genetic breeding, allowing the identification and implementation of novel and efficient strategies that enhance or optimize their use for efficiently controlling plant diseases. The present review focuses on main perspectives of application of R-genes and its co-players for the acquisition of resistance to pathogens in cultivated plants, with emphasis on biotechnological inferences, including transgenesis, cisgenesis, directed mutagenesis and gene editing, with examples of success and challenges to be faced.

Plant Elite Squad: First Defense Line and Resistance Genes - Identification, Diversity and Functional Roles.

Plants exhibit sensitive mechanisms to respond to environmental stresses, presenting some specific and non-specific reactions when attacked by pathogens, including organisms from different classes and complexity, as viroids, viruses, bacteria, fungi and nematodes. A crucial step to define the fate of the plant facing an invading pathogen is the activation of a compatible Resistance (R) gene, the focus of the present review. Different aspects regarding R-genes and their products are discussed, including pathogen recognition mechanisms, signaling and effects on induced and constitutive defense processes, splicing and post transcriptional mechanisms involved. There are still countless challenges to the complete understanding of the mechanisms involving R-genes in plants, in particular those related to the interactions with other genes of the pathogen and of the host itself, their regulation, acting mechanisms at transcriptional and post-transcriptional levels, as well as the influence of other types of stress over their regulation. A magnification of knowledge is expected when considering the novel information from the omics and systems biology.

Epigenetic Signals on Plant Adaptation: a Biotic Stress Perspective.

For sessile organisms such as plants, regulatory mechanisms of gene expression are vital, since they remain exposed to climatic and biological threats. Thus, they have to face hazards with instantaneous reorganization of their internal environment. For this purpose, besides the use of transcription factors, the participation of chromatin as an active factor in the regulation of transcription is crucial. Chemical changes in chromatin structure affect the accessibility of the transcriptional machinery and acting in signaling, engaging/inhibiting factors that participate in the transcription processes. Mechanisms in which gene expression undergoes changes without the occurrence of DNA gene mutations in the monomers that make up DNA, are understood as epigenetic phenomena. These include (1) post-translational modifications of histones, which results in stimulation or repression of gene activity and (2) cytosine methylation in the promoter region of individual genes, both preventing access of transcriptional activators as well as signaling the recruitment of repressors. There is evidence that such modifications can pass on to subsequent generations of daughter cells and even generations of individuals. However, reports indicate that they persist only in the presence of a stressor factor (or an inductor of the above-mentioned modifications). In its absence, these modifications weaken or lose heritability, being eliminated in the next few generations. In this review, it is argued how epigenetic signals influence gene regulation, the mechanisms involved and their participation in processes of resistance to biotic stresses, controlling processes of the plant immune system.

Validation of novel reference genes for reverse transcription quantitative real-time PCR in drought-stressed sugarcane.

One of the most challenging aspects of RT-qPCR data analysis is the identification of reliable reference genes. Ideally, they should be neither induced nor repressed under different experimental conditions. To date, few reference genes have been adequately studied for sugarcane (Saccharum spp.) using statistical approaches. In this work, six candidate genes ( αTUB, GAPDH, H1, SAMDC, UBQ, and 25S rRNA) were tested for gene expression normalization of sugarcane root tissues from drought-tolerant and -sensitive accessions after continuous dehydration (24 h). By undergoing different approaches (GeNorm, NormFinder, and BestKeeper), it was shown that most of them could be used in combinations for normalization purposes, with the exception of SAMDC. Nevertheless three of them (H1, αTUB, and GAPDH) were considered the most reliable reference genes. Their suitability as reference genes validated the expression profiles of two targets (AS and PFP α1), related to SuperSAGE unitags, in agreement with results revealed by previous in silico analysis. The other two sugarcane unitags (ACC oxidase and PIP1-1), after salt stress (100 mM NaCl), presented their expressions validated in the same way. In conclusion, these reference genes will be useful for dissecting gene expression in sugarcane roots under abiotic stress, especially in transcriptomic studies using SuperSAGE or RNAseq approaches.

Early transcriptional response of soybean contrasting accessions to root dehydration.

Drought is a significant constraint to yield increase in soybean. The early perception of water deprivation is critical for recruitment of genes that promote plant tolerance. DeepSuperSAGE libraries, including one control and a bulk of six stress times imposed (from 25 to 150 min of root dehydration) for drought-tolerant and sensitive soybean accessions, allowed to identify new molecular targets for drought tolerance. The survey uncovered 120,770 unique transcripts expressed by the contrasting accessions. Of these, 57,610 aligned with known cDNA sequences, allowing the annotation of 32,373 unitags. A total of 1,127 unitags were up-regulated only in the tolerant accession, whereas 1,557 were up-regulated in both as compared to their controls. An expression profile concerning the most representative Gene Ontology (GO) categories for the tolerant accession revealed the expression "protein binding" as the most represented for "Molecular Function", whereas CDPK and CBL were the most up-regulated protein families in this category. Furthermore, particular genes expressed different isoforms according to the accession, showing the potential to operate in the distinction of physiological behaviors. Besides, heat maps comprising GO categories related to abiotic stress response and the unitags regulation observed in the expression contrasts covering tolerant and sensitive accessions, revealed the unitags potential for plant breeding. Candidate genes related to "hormone response" (LOX, ERF1b, XET), "water response" (PUB, BMY), "salt stress response" (WRKY, MYB) and "oxidative stress response" (PER) figured among the most promising molecular targets. Additionally, nine transcripts (HMGR, XET, WRKY20, RAP2-4, EREBP, NAC3, PER, GPX5 and BMY) validated by RT-qPCR (four different time points) confirmed their differential expression and pointed that already after 25 minutes a transcriptional reorganization started in response to the new condition, with important differences between both accessions.

Natural antisense transcripts in plants: a review and identification in soybean infected with Phakopsora pachyrhizi SuperSAGE library.

Natural antisense ranscripts (NAT) are RNA molecules complementary to other endogenous RNAs. They are capable of regulating the expression of target genes at different levels (transcription, mRNA stability, translation, etc.). Such a property makes them ideal for interventions in organisms' metabolism. The present study reviewed plant NAT aspects, including features, availability and genesis, conservation and distribution, coding capacity, NAT pair expression, and functions. Besides, an in silico identification of NATs pairs was presented, using deepSuperSAGE libraries of soybean infected or not with Phakopsora pachyrhizi. Results showed that around 1/3 of the 77,903 predicted trans-NATs (by PlantsNATsDB database) detected had unitags mapped in both sequences of each pair. The same 1/3 of the 436 foreseen cis-NATs showed unitags anchored in both sequences of the related pairs. For those unitags mapped in NAT pairs, a modulation expression was assigned as upregulated, downregulated, or constitutive, based on the statistical analysis (P < 0.05). As a result, the infected treatment promoted the expression of 2,313 trans-NATs pairs comprising unitags exclusively from that library (1,326 pairs had unitags only found in the mock library). To understand the regulation of these NAT pairs could be a key aspect in the ASR plant response.

SymGRASS: a database of sugarcane orthologous genes involved in arbuscular mycorrhiza and root nodule symbiosis.

BACKGROUND: The rationale for gathering information from plants procuring nitrogen through symbiotic interactions controlled by a common genetic program for a sustainable biofuel production is the high energy demanding application of synthetic nitrogen fertilizers. We curated sequence information publicly available for the biofuel plant sugarcane, performed an analysis of the common SYM pathway known to control symbiosis in other plants, and provide results, sequences and literature links as an online database. METHODS: Sugarcane sequences and informations were downloaded from the nucEST database, cleaned and trimmed with seqclean, assembled with TGICL plus translating mapping method, and annotated. The annotation is based on BLAST searches against a local formatted plant Uniprot90 generated with CD-HIT for functional assignment, rpsBLAST to CDD database for conserved domain analysis, and BLAST search to sorghum's for Gene Ontology (GO) assignment. Gene expression was normalized according the Unigene standard, presented as ESTs/100 kb. Protein sequences known in the SYM pathway were used as queries to search the SymGRASS sequence database. Additionally, antimicrobial peptides described in the PhytAMP database served as queries to retrieve and generate expression profiles of these defense genes in the libraries compared to the libraries obtained under symbiotic interactions. RESULTS: We describe the SymGRASS, a database of sugarcane orthologous genes involved in arbuscular mycorrhiza (AM) and root nodule (RN) symbiosis. The database aggregates knowledge about sequences, tissues, organ, developmental stages and experimental conditions, and provides annotation and level of gene expression for sugarcane transcripts and SYM orthologous genes in sugarcane through a web interface. Several candidate genes were found for all nodes in the pathway, and interestingly a set of symbiosis specific genes was found. CONCLUSIONS: The knowledge integrated in SymGRASS may guide studies on molecular, cellular and physiological mechanisms by which sugarcane controls the establishment and efficiency of endophytic associations. We believe that the candidate sequences for the SYM pathway together with the pool of exclusively expressed tentative consensus (TC) sequences are crucial for the design of molecular studies to unravel the mechanisms controlling the establishment of symbioses in sugarcane, ultimately serving as a basis for the improvement of grass crops.

Expression Analysis of Sugarcane Aquaporin Genes under Water Deficit.

The present work is a pioneer study specifically addressing the aquaporin transcripts in sugarcane transcriptomes. Representatives of the four aquaporin subfamilies (PIP, TIP, SIP, and NIP), already described for higher plants, were identified. Forty-two distinct aquaporin isoforms were expressed in four HT-SuperSAGE libraries from sugarcane roots of drought-tolerant and -sensitive genotypes, respectively. At least 10 different potential aquaporin isoform targets and their respective unitags were considered to be promising for future studies and especially for the development of molecular markers for plant breeding. From those 10 isoforms, four (SoPIP2-4, SoPIP2-6, OsPIP2-4, and SsPIP1-1) showed distinct responses towards drought, with divergent expressions between the bulks from tolerant and sensitive genotypes, when they were compared under normal and stress conditions. Two targets (SsPIP1-1 and SoPIP1-3/PIP1-4) were selected for validation via RT-qPCR and their expression patterns as detected by HT-SuperSAGE were confirmed. The employed validation strategy revealed that different genotypes share the same tolerant or sensitive phenotype, respectively, but may use different routes for stress acclimation, indicating the aquaporin transcription in sugarcane to be potentially genotype-specific.

New insights in the sugarcane transcriptome responding to drought stress as revealed by superSAGE.

In the scope of the present work, four SuperSAGE libraries have been generated, using bulked root tissues from four drought-tolerant accessions as compared with four bulked sensitive genotypes, aiming to generate a panel of differentially expressed stress-responsive genes. Both groups were submitted to 24 hours of water deficit stress. The SuperSAGE libraries produced 8,787,315 tags (26 bp) that, after exclusion of singlets, allowed the identification of 205,975 unitags. Most relevant BlastN matches comprised 567,420 tags, regarding 75,404 unitags with 164,860 different ESTs. To optimize the annotation efficiency, the Gene Ontology (GO) categorization was carried out for 186,191 ESTs (BlastN against Uniprot-SwissProt), permitting the categorization of 118,208 ESTs (63.5%). In an attempt to elect a group of the best tags to be validated by RTqPCR, the GO categorization of the tag-related ESTs allowed the in silico identification of 213 upregulated unitags responding basically to abiotic stresses, from which 145 presented no hits after BlastN analysis, probably concerning new genes still uncovered in previous studies. The present report analyzes the sugarcane transcriptome under drought stress, using a combination of high-throughput transcriptome profiling by SuperSAGE with the Solexa sequencing technology, allowing the identification of potential target genes during the stress response.

Identification of plant protein kinases in response to abiotic and biotic stresses using SuperSAGE.

Plants are sessile organisms subjected to many environmental adversities. For their survival they must sense and respond to biotic and abiotic stresses efficiently. During this process, protein kinases are essential in the perception of environmental stimuli, triggering signaling cascades. Kinases are among the largest and most important gene families for biotechnological purposes, bringing many challenges to the bioinformaticians due to the combination of conserved domains besides diversified regions. Cowpea [Vigna unguiculata (L.) Walp.] is an important legume that is adapted to different agroclimatic conditions, including drought, humidity and a range of temperatures. For this crop, the association of the SuperSAGE method with high-throughput sequencing technology would generate reliable transcriptome profiles with millions of tags counted and statistically analyzed. An approach evaluating biotic and abiotic stresses was carried out generating over 13 million cowpea SuperSAGE tags available from leaves/roots of plants under abiotic (mechanical injury and salinity) or biotic (CABMV, Cowpea aphid born mosaic virus) stresses. The annotation and identification of tags linked by BlastN to previously well described ESTs, allowed the posterior identification of kinases. The annotation efficiency depended on the database used, with the KEGG figuring as a good source for annotated ESTs especially when complemented by an independent Gene Ontology categorization, as well as the Gene Index using selected species. The use of different approaches allowed the identification of 1,350 kinase candidates considering biotic libraries and 2,268 regarding abiotic libraries, based on a combination of both, adequate descriptions and GO terms. Additional searches in kinase specific databases allowed the identification of a relatively low number of additional kinases, uncovering the lack of kinase databases for non-model organisms, especially plants. Concerning the kinase families, a total of 713 potential kinases were classified into 13 families of the CMGC and STE groups. Concerning the differentially expressed kinases, 169 of the 713 potential kinases were identified (p < 0.05), 100 up- and 69 down-regulated when comparing distinct libraries, allowing the generation of a comprehensive panel of the differentially expressed kinases under biotic and abiotic stresses in a non-model plant as cowpea.