Optimization of cDNA-AFLP experiments using genomic sequence data

Emerging tools, concepts and ideas to track the modulator genes underlying plant drought adaptive traits: An overview

Crop vulnerability to multiple abiotic stresses is increasing at an alarming rate in the current global climate change scenario, especially drought. Crop improvement for adaptive adjustments to accomplish stress tolerance requires a comprehensive understanding of the key contributory processes. This requires the identification and careful analysis of the critical morpho-physiological plant attributes and their genetic control. In this review we try to discuss the crucial traits underlying drought tolerance and the various modes followed to understand their molecular level regulation. Plant stress biology is progressing into new dimensions and a conscious attempt has been made to traverse through the various approaches and checkpoints that would be relevant to tackle drought stress limitations for sustainable crop production.

Differential gene expression in tea (Camellia sinensis L.) calli with different morphologies and catechin contents

Tea (Camellia sinensis) is a commercially important crop that contains valuable secondary metabolites. To understand the molecular regulation of secondary metabolism in tea, we selected and analyzed two cell lines of tea callus (Yunjing63Y and Yunjing63X) that showed different morphological characteristics and catechin contents. Yunjing63Y callus was yellow and tight, while yunjing63X callus was white and loose. HPLC analyses showed that Yunjing63Y contained 3.71 times higher levels of catechins than Yunjing63X. Using cDNA amplified fragment-length polymorphism (cDNA-AFLP) we identified 68 genes that were differentially expressed between the two lines. Of the 68 differentially expressed ESTs, 40 showed higher expressions in Yunjing63Y and 28 showed higher expressions in Yunjing63X. BLASTX comparisons classified these ESTs into seven functional groups; phenylpropanoid metabolism (2.9%), UDPG-dependent glucosyl transferase (8.8%), transcription factors (11.8%), transporters (13.2%), signal transduction (19.1%), other metabolism (26.5%), and unknown (17.7%). We used qRT-PCR to validate the expression of genes and ESTs, and found that genes associated with flavan-3-ols biosynthesis and metabolism were expressed at higher levels in Yunjing63Y than in Yunjing63X. In addition, the expression of ESTs associated with flavonoid biosynthesis, regulation and transport were higher in Yunjing63Y than in Yunjing63X. The full-length cDNA of a EST coding for a putative MYB transcription factor was amplified using rapid amplification of cDNA ends (RACE). The resulting 1270 bp long cDNA, named CsMYB1, contained a 933-bp ORF encoding a 310-amino acid protein with a predicted molecular weight of 105.27 kDa and a predicted isoelectric point of 4.85 and showed highest homology to plant MYBs likely involved in stress signaling.

Isolation and expression analysis of partial sequences of heavy metal transporters from Brassica juncea by coupling high throughput cloning with a molecular fingerprinting technique

Heavy metal transporters play a key role in regulating metal accumulation and transport in plants. These are important candidate genes to study in metal tolerant and accumulator plants for their potential use in environmental clean up. We coupled a degenerate primer-based RT-PCR approach with a molecular fingerprinting technique based on amplified rDNA restriction analysis (ARDRA) to identify novel ESTs corresponding to heavy metal transporters from metal accumulator Brassica juncea. We utilized this technique to clone several family members of natural resistance-associated macrophage proteins (NRAMP) and yellow stripe-like proteins (YSL) in a high throughput manner to distinguish between closely related isoforms and/or allelic variants from the allopolyploid B. juncea. Partial clones of 23 Brassica juncea NRAMPs and 27 YSLs were obtained with similarity to known Arabidopsis thaliana and Noccaea (Thlaspi) caerulescens NRAMP and YSL genes. The cloned transporters showed Brassica-specific changes in domains, which can have important functional consequences. Semi-quantitative RT-PCR-based expression analysis of chosen members indicated that even closely related isoforms/allelic variants of BjNRAMP and BjYSL have distinct tissue-specific and metal-dependent expressions which might be essential for adaptive fitness and heavy metal tolerance. Consistent to this, BjYSL6.1 and BjYSL5.8 were found to show elevated expressions specifically in cadmium-treated shoots and lead-treated roots of B. juncea, respectively.

Identification of genes involved in stem rust resistance from wheat mutant D51 with the cDNA-AFLP technique

Wheat (Triticum aestivum L.) stem rust caused by Puccinia graminis f. sp. tritici is one of the main diseases of wheat worldwide. Wheat mutant line D51, which was derived from the highly susceptible cultivar L6239, shows resistance to the prevailing races 21C3CPH, 21C3CKH, and 21C3CTR of P. graminis f. sp. tritici in China. In this study, we used the cDNA-AFLP technology to identify the genes that are likely involved in the stem rust resistance. EcoRI/MseI selective primers were used to generate approximately 1920 DNA fragments. Seventy five differentially transcribed fragments (3.91%) were identified by comparing the samples of 21C3CPH infected D51 with infected L6239 or uninfected D51. Eleven amplified cDNA fragments were sequenced. Eight showed significant similarity to known genes, including TaLr1 (leaf rust resistance gene), wlm24 (wheat powdery mildew resistance gene), stress response genes and ESTs of environment stress of tall fescue. These identified genes are involved in plant defense response and stem rust resistance and need further research to determine their usefulness in breeding new resistance cultivars.

Eukaryotic transcriptomics in silico: optimizing cDNA-AFLP efficiency

Background

Complementary-DNA based amplified fragment length polymorphism (cDNA-AFLP) is a commonly used tool for assessing the genetic regulation of traits through the correlation of trait expression with cDNA expression profiles. In spite of the frequent application of this method, studies on the optimization of the cDNA-AFLP assay design are rare and have typically been taxonomically restricted. Here, we model cDNA-AFLPs on all 92 eukaryotic species for which cDNA pools are currently available, using all combinations of eight restriction enzymes standard in cDNA-AFLP screens.

Results

In silco simulations reveal that cDNA pool coverage is largely determined by the choice of individual restriction enzymes and that, through the choice of optimal enzyme combinations, coverage can be increased from <40% to 75% without changing the underlying experimental design. We find evidence of phylogenetic signal in the coverage data, which is largely mediated by organismal GC content. There is nonetheless a high degree of consistency in cDNA pool coverage for particular enzyme combinations, indicating that our recommendations should be applicable to most eukaryotic systems. We also explore the relationship between the average observed fragment number per selective AFLP-PCR reaction and the size of the underlying cDNA pool, and show how AFLP experiments can be used to estimate the number of genes expressed in a target tissue.

Conclusion

The insights gained from in silico screening of cDNA-AFLPs from a broad sampling of eukaryotes provide a set of guidelines that should help to substantially increase the efficiency of future cDNA-AFLP experiments in eukaryotes. In silico simulations also suggest a novel use of cDNA-AFLP screens to determine the number of transcripts expressed in a target tissue, an application that should be invaluable as next-generation sequencing technologies are adapted for differential display.

Identification and functional validation of a unique set of drought induced genes preferentially expressed in response to gradual water stress in peanut

Peanut, found to be relatively drought tolerant crop, has been the choice of study to characterize the genes expressed under gradual water deficit stress. Nearly 700 genes were identified to be enriched in subtractive cDNA library from gradual process of drought stress adaptation. Further, expression of the drought inducible genes related to various signaling components and gene sets involved in protecting cellular function has been described based on dot blot experiments. Fifty genes (25 regulators and 25 functional related genes) selected based on dot blot experiments were tested for their stress responsiveness using northern blot analysis and confirmed their nature of differential regulation under different field capacity of drought stress treatments. ESTs generated from this subtracted cDNA library offered a rich source of stress-related genes including signaling components. Additional 50% uncharacterized sequences are noteworthy. Insights gained from this study would provide the foundation for further studies to understand the question of how peanut plants are able to adapt to naturally occurring harsh drought conditions. At present functional validation cannot be deemed in peanut, hence as a proof of concept seven orthologues of drought induced genes of peanut have been silenced in heterologous N. benthamiana system, using virus induced gene silencing method. These results point out the functional importance for HSP70 gene and key regulators such as Jumonji in drought stress response.

Improved coverage of cDNA-AFLP by sequential digestion of immobilized cDNA

BACKGROUND

cDNA-AFLP is a transcriptomics technique which does not require prior sequence information and can therefore be used as a gene discovery tool. The method is based on selective amplification of cDNA fragments generated by restriction endonucleases, electrophoretic separation of the products and comparison of the band patterns between treated samples and controls. Unequal distribution of restriction sites used to generate cDNA fragments negatively affects the performance of cDNA-AFLP. Some transcripts are represented by more than one fragment while other escape detection, causing redundancy and reducing the coverage of the analysis, respectively.

RESULTS

With the goal of improving the coverage of cDNA-AFLP without increasing its redundancy, we designed a modified cDNA-AFLP protocol. Immobilized cDNA is sequentially digested with several restriction endonucleases and the released DNA fragments are collected in mutually exclusive pools. To investigate the performance of the protocol, software tool MECS (Multiple Enzyme cDNA-AFLP Simulation) was written in Perl. cDNA-AFLP protocols described in the literature and the new sequential digestion protocol were simulated on sets of cDNA sequences from mouse, human and Arabidopsis thaliana. The redundancy and coverage, the total number of PCR reactions, and the average fragment length were calculated for each protocol and cDNA set.

CONCLUSION

Simulation revealed that sequential digestion of immobilized cDNA followed by the partitioning of released fragments into mutually exclusive pools outperformed other cDNA-AFLP protocols in terms of coverage, redundancy, fragment length, and the total number of PCRs. Primers generating 30 to 70 amplicons per PCR provided the highest fraction of electrophoretically distinguishable fragments suitable for normalization. For A. thaliana, human and mice transcriptome, the use of two marking enzymes and three sequentially applied releasing enzymes for each of the marking enzymes is recommended.

AFLP-based transcript profiling (cDNA-AFLP) for genome-wide expression analysis

Although DNA microarrays are currently the standard tool for genome-wide expression analysis, their application is limited to organisms for which the complete genome sequence or large collections of known transcript sequences are available. Here, we describe a protocol for cDNA-AFLP, an AFLP-based transcript profiling method that allows genome-wide expression analysis in any species without the need for prior sequence knowledge. In essence, the cDNA-AFLP method involves reverse transcription of mRNA into double-stranded cDNA, followed by restriction digestion, ligation of specific adapters and fractionation of this mixture of cDNA fragments into smaller subsets by selective PCR amplification. The resulting cDNA-AFLP fragments are separated on high-resolution gels, and visualization of cDNA-AFLP fingerprints is described using either a conventional autoradiography platform or an automated LI-COR system. Observed differences in band intensities between samples provide a good measure of the relative differences in the gene expression levels. Identification of differentially expressed genes can be accomplished by purifying cDNA-AFLP fragments from sequence gels and subsequent sequencing. This method has found widespread use as an attractive technology for gene discovery on the basis of fragment detection and for temporal quantitative gene expression analysis. The protocol can be completed in 3-4 d.

Computer-assisted identification of multitrace electrophoretic patterns in differential display experiments

Modern multicapillary devices allow researchers to address increasingly complex biological questions involving comparisons of gene expression patterns across electrophoretic samples under various experimental conditions. As labor-intensive visual evaluation of the electrophoretic results is often the bottleneck of large-scale differential display (DD) studies, one way to further streamline this process is to focus only on a highly compressed list of the most potential patterns that are likely to provide reliable findings. To enable the identification of such candidate patterns, we present a computer-assisted method for objective ranking of multitrace peak patterns in DD experiments. The fundamental component of the multitrace pattern ranking method (MRANK) is the multiple alignment algorithm that allows for discovery of patterns involving sets of peak complexes from various electrophoretic samples. A score value is attached to each detected pattern which characterizes how accurately the pattern resembles the desired pattern query, freely defined by the researcher. The ranked pattern list produced by MRANK is validated against visual evaluation in terms of detecting and ranking a group of relevant patterns in a DD analysis of T-helper cell differentiation. We demonstrate high enrichment of the desired patterns on top of the score-ranked list (e.g., 90% of the visually selected patterns are discovered by looking through the first 3% of patterns in the ranked list of all patterns). The results suggest that a substantial amount of manual labor can be saved without compromising the accuracy of the findings by prioritizing the patterns according to MRANK output in the visual confirmation phase.