Background: Aegilops variabilis No.1 is highly resistant to cereal cyst nematode (CCN). However, a lack of genomic information has restricted studies on CCN resistance genes in Ae. variabilis and has limited genetic applications in wheat breeding. Results: Using RNA-Seq technology, we generated a root transcriptome at a sequencing depth of 4.69 gigabases of Ae. variabilis No. 1 from a pooled RNA sample. The sample contained equal amounts of RNA extracted from CCN-infected and untreated control plants at three time-points. Using the Trinity method, nearly 52,081,238 high-quality trimmed reads were assembled into a non-redundant set of 118,064 unigenes with an average length of 500 bp and an N50 of 599 bp. The total assembly was 59.09 Mb of unique transcriptome sequences with average read-depth coverage of 33.25×. In BLAST searches of our database against public databases, 66.46% (78,467) of the unigenes were annotated with gene descriptions, conserved protein domains, or gene ontology terms. Functional categorization further revealed 7,408 individual unigenes and three pathways related to plant stress resistance. Conclusions: We conducted high-resolution transcriptome profiling related to root development and the response to CCN infection in Ae. variabilis No.1. This research facilitates further studies on gene discovery and on the molecular mechanisms related to CCN resistance.

The proteins containing the TIFY domain belong to a plant-specific family of putative transcription factors and could be divided into four subfamilies: ZML, TIFY, PPD and JAZ. They not only function as key regulators of jasmonate hormonal response, but are also involved in responding to abiotic stress. In this study, we identified 24 TIFY genes (PeTIFYs) in Moso bamboo (Phyllostachys edulis) of Poaceae by analyzing the whole genome sequence. OnePeTIFYbelongs to TIFY subfamily, 18 and five belong to JAZ and ZML subfamilies, respectively. Two equivocal gene models were re-predicted and a putative retrotransposition event was found in a ZML protein. The distribution and conservation of domain or motif, and gene structure were also analyzed. Phylogenetic analysis with TIFY proteins ofArabidopsisandOryza sativaindicated that JAZ subfamily could be further divided to four groups. Evolutionary analysis revealed intragenomic duplication and orthologous relationship betweenP. edulis,O. sativa, andB. distachyon. Calculation of the non-synonymous (Ka) and synonymous (Ks) substitution rates and their ratios indicated that the duplication ofPeTIFYmay have occurred around 16.7 million years ago (MYA), the divergence time of TIFY family among theP. edulis-O. sativa,P. edulis-B. distachyon,andO. sativa-B. distachyonwas approximately 39 MYA, 39 MYA, and 45 MYA, respectively. They appear to have undergone extensive purifying selection during evolution. Transcriptome sequencing revealed that more than 50% ofPeTIFYgenes could be up-regulated by cold and dehydration stresses, and somePeTIFYsalso share homology to know TIFYs involved in abiotic stress tolerance. Our results made insights into TIFY family of Moso bamboo, an economically important non-timber forest resource, and provided candidates for further identification of genes involved in regulating responses to abiotic stress.

Late embryogenesis abundant (LEA) proteins have been identified in a wide range of organisms and are believed to play a role in the adaptation of plants to stress conditions. In this study, we performed genome-wide identification of LEA proteins and their coding genes in Moso bamboo (Phyllostachys edulis) of Poaceae. A total of 23 genes encoding LEA proteins (PeLEAs) were found in P. edulis that could be classified to six groups based on Pfam protein family and homologous analysis. Further in silico analyses of the structures, gene amount, and biochemical characteristics were conducted and compared with those of O. sativa (OsLEAs), B. distachyon (BdLEAs), Z. mays (ZmLEAs), S. bicolor (SbLEAs), Arabidopsis, and Populus trichocarpa. The less number of PeLEAs was found. Evolutionary analysis revealed orthologous relationship and colinearity between P. edulis, O. sativa, B. distachyon, Z. mays, and S. bicolor. Analyses of the non-synonymous (Ka) and synonymous (Ks)substitution rates and their ratios indicated that the duplication of PeLEAs may have occurred around 18.8 million years ago (MYA), and divergence time of LEA family among the P. edulis-O. sativa and P. edulis–B. distachyon, P. edulis-S. bicolor, and P. edulis-Z. mays was approximately 30 MYA, 36 MYA, 48 MYA, and 53 MYA, respectively. Almost all PeLEAs contain ABA- and (or) stress-responsive regulatory elements. Further RNA-seq analysis revealed approximately 78% of PeLEAs could be up-regulated by dehydration and cold stresses. The present study makes insights into the LEA family in P. edulis and provides inventory of stress-responsive genes for further functional validation and transgenic research aiming to plant genetic improvement of abiotic stress tolerance.