BACKGROUND
all orchids produce minute, endosperm- lacking seeds and are dependent on fungal colonization for germination and growth into an underground heterotrophic, achlorophyllous stage called a protocorm under natural conditions  <cit> .
after this fully mycotrophic stage, most species shift to autotrophy at adult stage with the development of photosynthetic organs, but approximately  <dig> species worldwide stay fully mycoheterotrophic   <cit> .
recently, some chlorophyllous species in orchidaceae were verified to obtain carbon at adult stage not only from their photosynthetic activity, but also from mycorrhizal fungi, which were called mixotrophs   <cit> .
among these, several terrestrial orchids such as c. lancifolium and c. goeringii were shown to be mx nutritional mode through isotopic analyses  <cit> .
the dependency of green orchids on mycorrhizal fungi throughout their life cycle prompts us to study these mycorrhizal interactions in more detail.
in nature, orchids largely associate with members of the ‘rhizoctonia’ complex, which contains three now taxonomically disparate agaricomycetes  taxa: sebacinales, ceratobasidiaceae and tulasnellaceae.
tulasnellaceae are the most frequently found ‘rhizoctonias’ in green orchids  <cit> .
besides orchid mycorrhizal fungi , recent studies also demonstrated that isolates of umbelopsis, as non-mycorrhizal beneficial endophytes, probably co-existed in the roots of wild orchids  <cit> .
c. hybridum collectively contains a large number of cultivars with large flowers.
over  <dig> species in cymbidium, including c. sinense, c. goeringii, c. faberi, c. ensifolium and c. kanran, were used in their breeding programs.
in our previous work, one omf isolate ml <dig> and one non-mycorrhizal beneficial isolate zh3a- <dig> were isolated from roots of wild c. spp.
adult plants and identified as e. repens  and u. nana , respectively.
here, we present that ml <dig> forms typical mycorrhiza with c. hybridum.
individual  or combinative  inoculation significantly improved the vegetative growth of c. hybridum plantlets compared with mock-inoculated controls.
although orchid mycorrhizal  symbiosis is an important feature in the life cycle of orchids , it has long been considered to be an atypical mycorrhizal association, with the fungus deriving little benefit from the orchid host  <cit> .
but cameron et al.
<cit>  reported that approximately  <dig> – <dig> % of the labeled carbon provided to the orchids was passed to the fungal partner  <cit> .
thus adult orchid mycorrhizas may represent a truly mutualistic interaction similar to ectomycorrhizae  and arbuscular mycorrhizae  associations  <cit> .
compared with the numerous transcriptional analyses to identify genes involved in am  and ecm symbiosis  , the molecular mechanisms behind om symbiosis are still poorly understood .
recently, different gene expression or proteome changes of orchid protocorms in response to mycorrhizal fungi were studied and mutualistic plant-fungus relationship was suggested .
as for orchids at the adult stage, valadares  <cit>  used rna-seq to determine the gene expression profile of wild oeceoclades maculata mycorrhizal roots in comparison to non-mycorrhizal controls  <cit> .
while this study provides new insights into our understanding of the om symbiosis for adult orchids, the symbiosis seems to be highly complex as wild orchids generally associate with multiple symbionts even within a single root  <cit> .
although genomic sequence resources currently available for orchids are limited, large expressed transcripts from phalaenopsis aphrodite, p. equestris, p. bellina, oncidium gower ramsey, erycina pusilla and c. sinense were generated in recent years with the rapid advances in dna sequencing technology known as next-generation sequencing  and the development of bioinformatic processes, especially for de novo assembly without reference genome information .
whereas no genetic information on orchid mycorrhizal fungi was available when our analyses began, large numbers of assembled transcripts and the genome information of t. calospora had been recently generated as part of a doe jgi community sequencing program.
this work has allowed us to assign some putative fungal genes to mycorrhizal fungi ml <dig> by blastn searches of these data.
in this paper, we used ngs approach and illumina digital gene expression  technology to obtain genome-wide insight into the transcriptomic responses of c. hybridum under specific interactions with mycorrhizal or/and other beneficial fungi.
to our knowledge, this represents the first genome-wide gene expression profile of orchid plants in response to mycorrhizal or non-mycorrhizal fungi.
the data enabled us to get new insight into the mechanisms underlying plant responses to mycorrhizal and non-mycorrhizal beneficial fungi.
the transcriptome for c. hybridum could also be used as an important resource for exploring other biological processes in orchids.
methods
biological materials and inoculation
the micro-propagated plantlets of c. hybridum ‘golden boy’ were provided by the guangdong key lab of ornamental plant germplasm innovation and utilization, china.
the propagation procedure for orchid plantlets was described in  <cit> .
the two fungal isolates ml <dig> and zh3a- <dig> were preserved in the guangdong key laboratory for innovative development and utilization of forest plant germplasm.
inocula were produced by growing each strain on 90 mm petri dishes containing potato-dextrose agar  for 10 days.
the plantlets with 3– <dig> roots about  <dig> -3 cm long were then transferred to de basal medium, an in vitro mycorrhizal induction medium in orchids  <cit> .
the co-cultivation de medium contained  <dig>  mm cacl <dig>   <dig>  mm mgs <dig>   <dig>  mm k2s <dig>   <dig>  mm kh2p <dig>  100 μm fes <dig>  25 μm h3bo <dig>  33 μm mncl <dig>   <dig>  μm znso <dig>   <dig>  μm namoo <dig>  140 μm na2edta, supplemented with yeast extract 1 g · l- <dig>  soluble starch 9 g · l- <dig> and agar 6 g · l- <dig>  each plantlet was inoculated with a 7 mm plug excised from an edge of an actively growing colony on pda of strain ml <dig>  zh3a- <dig> or ml01+ zh3a- <dig> the control plantlets were mock-inoculated with a plug excised from a pda plate without fungus.
all the treated and control plantlets were then incubated under a photoperiod of 16 hrs light at a temperature of 22°c.
symbiotic effects on vegetative growth of c. hybridum
after 45 days post- inoculation , ninety plantlets of each treatment were totally harvested and randomly divided into three groups, which represented three biological replicates.
the symbiotic effects on the vegetative growth of the hosts were assessed by total increased fresh weight and dry biomass of the plantlets.
the percentage of root growth  and mycorrhizal dependency  were also compared among different treatments.
md was calculated by the following formula:  × 100/shoot dry wt.
of mycorrhizal plants as described in  <cit> .
for statistical analysis of data, spss v <dig>  was employed for one-way analysis of variance  and comparisons among means were made using duncan’s multiple range test, calculated at p <  <dig> .
root staining and histology
the successful inoculation with strains ml <dig> and zh3a- <dig> and the absence of fungal contaminants in mock-treated roots were confirmed by trypan blue staining  <cit> .
to confirm their respective colonization pattern, semi-thin sections of root samples for each treatment were prepared and stained with a  <dig> % aqueous solution of toluidine blue for light microscopy .
reference cdna library construction
at  <dig>   <dig>   <dig> and  <dig> dpi, orchid roots of mz- treated or mock-inoculated control plantlets were harvested by carefully removing the substrate and mycelia, and the samples were frozen in liquid nitrogen and stored at -80°c.
total rna from each sample was extracted separately with column plant total rnaout kit according to the manufacturer’s protocol .
the rna was treated with dnase  and then quantified using a nanodrop 2000™ spectrophotometer .
equal amounts of total rna from mz- treated and mock-inoculated samples at each time point were pooled to produce minimum 20 μg total rna for constructing a cdna library for establishing the reference transcriptome.
mrna-sequencing and cdna de novoassembly
the library construction and sequencing were performed by the beijing genomics institute   and the general experimental pipeline was summarized by previously published papers  <cit> .
first, the cdna library was pair-end sequenced by illumina hiseq <dig> .
after removing adaptor sequences, empty reads and low quality sequences , we obtained the clean reads for de novo assembly.
three de novo assembly software were used including soapdenovo with optimized k-mer length of  <dig>  abyss  <cit>  with k-mer length of  <dig> and trinity  <cit>   with default k-mer length of  <dig>  sequence directions of the resulting unigenes were validated by blastx alignment  between unigenes and public protein databases with a priority order of nr , swiss-prot, kyoto encyclopedia of genes and genomes database , and cog.
when a unigene happens to be unaligned to none of the above databases, estscan will be introduced to predict its coding regions as well as to decide its sequence direction.
after sequencing and assembling, we obtained three de novo assembled transcriptome database with three different assembly methods.
the number of assembled transcripts, total bases of transcripts, average length of all unigenes, n <dig>  number of long-transcripts , and number of reads that could be mapped back to transcripts  were compared to evaluate the assembly quality with different methods and select the optimum reference transcriptome for further study.
the expression level of unigene was measured by the number of clean reads mapped to each annotated unigene normalized to rpkm   <cit>  and adjusted using an additional normalizing factor  <cit> .
functional annotation
the protein functional annotation of unigenes was given according to the known proteins with the highest sequence similarity by blastx to public protein databases mentioned above.
domain-based alignments were carried out against the cog database to predict and classify possible functions of unigenes.
we further used the blast2go program  to get go annotation of unigenes.
validation of the reference assembly
to validate the assembled reference transcripts,  <dig> unigenes were randomly chosen for further rt-pcr and sanger sequencing.
gene specific primers were designed according to the assembled transcripts with primer premier software  .
the reaction was performed in 25 μl, containing  <dig>  μl cdna template,  <dig>  μl ex taq mix , 1 μl of each 20 μm forward and reverse gene-specific primers and 10 μl of pcr-grade water.
touch-down pcr reactions were performed as follows: 95°c for 5 min before cycling  <dig> rounds of 94°c for 30 s, 65°c for 30 s  and 72°c for 1 min 30 s – followed by 15 cycles of 94°c for 30 s, 54°c for 30 s and 72°c for 1 min 30 s and finally a 72°c extension step for 10 min.
the specificity of products was inspected by 1% agarose gel electrophoresis before they were purified with ezgene™ gel/pcr extraction kit .
the purified pcr products were sanger-sequenced by the beijing genomics institute  .
identification of differentially expressed genes  using illumina short reads
to perform quantitative comparisons of the transcript variation with respect to different beneficial fungi, total rna samples generated from roots treated with ml <dig> , zh3a- <dig>  or ml01 + zh3a- <dig>  and mock-inoculated controls  at  <dig> dpi were used to construct un-normalised cdna libraries which were sequenced by illumina hiseq <dig> .
the cdna library construction, sequencing, transcript profiling and mapping back to the newly assembled reference transcriptome with trinity were carried out as a custom service .
the expression level of each assembled transcript in different samples was measured as rpkm values.
seqmap  <cit>  was used for read mapping and rseq  <cit>  was applied for rpkm based expression measurement.
abundance data from different samples were collected for each transcript.
in our work, the changes in relative abundance for unigenes between two samples were screened with the threshold: q-value ≤  <dig>  and log2ratio ≥  <dig> or ≤ - <dig>
quantitative rt-pcr  validation
qrt-pcr was further carried out to analyze the expression of twenty- six degs identified by rna-seq technology.
gene specific primers were designed according to the cdnas with primer premier software  .
for rna extraction of each treatment, three biological replicates were collected independently and immediately frozen in liquid nitrogen.
one microgram of total rna was reverse-transcribed in a 20 μl reaction mixture from a primescript ii 1st strand cdna synthesis kit .
qrt-pcr was performed in a 25 μl reaction mixture containing 2 × sybr master premix ex taq ii  <dig>  μl , 1 μl of cdna template , and 1 μl of each corresponding primer for the gene of interest.
qrt-pcr of three biological replicates for each sample was performed using a lightcycler  <dig> ii system with its relative quantification software  based on the delta-delta-ct method .
qrt-pcr was performed for 5 s at 95°c, 10 s at 56°c, and 20 s at 72°c.
the cdna samples were standardized to two reference genes: b <dig>  and b <dig> .
RESULTS
co-cultivation of c. hybridumplantlets with different fungi and their respective symbiotic effects on host vegetative growth
when co-cultivated with c. hybridum plantlets, the hyphae of the two isolates ml <dig> and zh3a- <dig> began to grow at  <dig> dpi and first contacted the hosts’ roots at  <dig> dpi.
after  <dig> dpi, the co-cultivation media of plantlets inoculated with ml <dig> changed from transparent to opaque, whereas the actively growing white fluffy hyphae of isolate zh3a- <dig> fully covered the plate.
the hyphal growth of isolate zh3a- <dig> was retarded by the co-inoculated isolate ml <dig>  which showed space at the colony margin .
after  <dig> dpi, ml <dig>  zh3a- <dig> and mz treatments all significantly improved the total fresh weight and dry biomass of c. hybridum plantlets compared to the control .table  <dig>
effect of inoculation with
epulorhiza repens
isolate ml01and/or
umbelopsis nana
isolate zh3a- <dig> on the growth parameters of
cymbidium hybridum
plantlets at  <dig> dpi
all the data were mean of three biological replicates, ± standard deviation.
: percentage of mycorrhizal dependency; prg : the elongated root number/total root number × 100).
values followed by different lower-case letters within a column are significantly different at p <  <dig>  according to dmrt.
detection the colonization feature of different beneficial fungi on roots of c. hybridumplantlets
in an initial observation, the hyphae of two isolates began to contact with orchids roots at  <dig> dpi.
after that time point, we stained the host's roots with each treatment every two days to monitor the colonization process of different fungi and the root staining results of each treatment at  <dig>   <dig>   <dig> and  <dig> dpi were presented in additional file 4: figure s <dig>  at  <dig> dpi, the first contacted roots were stained to confirm the hyphal penetration by isolates ml <dig> or zh3a- <dig>  after  <dig> dpi, orchid roots were extensively colonized by these fungi.
at  <dig> dpi, the symbiotic roots of orchid plantlets showed uneven staining owing to the rapid growth of roots.
to confirm the root staining results and characterize the exact colonization pattern of these beneficial fungi, we further analyzed the microscopic structure of the formed symbionts at  <dig> dpi.
root cross sections of c. hybridum plantlets generally included root hairs, velamen , exodermis , cortex , endodermis , pericycle , vascular cylinder and pith .
at  <dig> dpi, typical pelotons were formed by isolate ml <dig> in the host cortex cells  whereas hyphae of isolate zh3a- <dig> only colonized the velamen cells at regular intervals .
roots co-inoculated with isolates ml <dig> and zh3a- <dig> contained newly formed and degenerated pelotons, and also typical aggregated hyphae of isolate zh3a- <dig>  demonstrating the co-infection of two isolates .
none of these fungal structures were observed in control root tissue .figure  <dig>
microscopic structure of the symbionts formed by different
c. hybridum
- beneficial fungi interaction at  <dig> dpi.
mock-inoculated control  and inoculated with isolate zh3a- <dig> , ml01 + zh3a- <dig>  and ml <dig> .
ah, aggregated hyphae of isolate zh3a-3; co, cortex cells; dh, degenerated hyphae; ex, exodermis; mh, fungal hyphae of isolate ml01; pe, pelotons; rh, root hair; va, vascular cylinder; ve, velamen cells.
bar, 50 μm
root transcriptome profile of c. hybridiumgenerated by mrna-sequencing
up to date, none of genome sequences for orchidaceae was available in public database.
rna-seq, a cost-effective and highly efficient next-generation sequencing technology, provides high throughput and is therefore commonly used for de novo transcriptome assemblies in non-model species, including orchid species .
to analyze the transcript variation in c. hybridium colonized by beneficial fungi, pooled total rna from eight samples including mz-treated and mock-inoculated plantlets at each time point were used to generate an informative reference transcriptome database.
totally,  <dig> , <dig>  base pairs raw data were generated using the illumina hiseq <dig>  yielding  <dig> , <dig> clean reads that were 90 bp in length.
the raw reads data are available at the ncbi sra database with the accession number sra <dig>
three publicly available assemblers as soapdenovo, abyss and trinity were used to de novo assemble short-read rna-seq data into transcripts.
the assembly software tools as soapdenovo and abyss were chosen because of their previous wide application for the similar de novo assembly of rna-seq data sets  <cit> .trinity method , one of the newest de novo assembly methods, was reported to reconstruct of the majority of full-length transcripts in a sample from rna-seq reads directly, across a broad range of transcript levels  <cit> .
recent publications also showed that trinity had a consistently better performance than the other assemblers as soapdenovo, abyss, trans-abyss, oases with single k-mer value, when measured by transcript accuracy, integrity and completeness, and sensitivity to assemble transcripts from low to high expression levels  <cit> .
the outcomes of three assemblers in our study were summarized in table  <dig>  as shown in table  <dig>  de novo assembly using soapdenovo yielded a total of  <dig>  unique transcripts , the largest number of transcripts.
in contrast to the soapdenovo and abyss assembly, trinity assembled the longest mean transcript length , the largest n <dig> , the largest number of long-transcripts , and the highest percentage of rmbt .
we therefore elected to use the trinity assembly to retain maximum information for further transcript analysis and the trinity assembly is also available in additional file  <dig> table  <dig>
summary of
de novo
sequence assembly
annotation
to validate and annotate the assembled unigenes, the  <dig>  unigenes generated by trinity were subjected to blastx searches  against public protein databases.
we found  <dig>  ,  <dig>  ,  <dig>   sequences with matches in nr, swiss-prot, and kegg databases, respectively.
in total,  <dig>   sequences were annotated and another  <dig>   unigenes showed no homology to known sequences in these databases.
among the annotated sequences, approximately  <dig> transcripts shared the highest similarity to fungal or bacterial sequences in nr or swiss-prot database.
the unigenes homologous to known sequences in public databases were further annotated with go terms using blast2go.
a total of  <dig>   unigenes were assigned  <dig>  go term annotations in three categories: biological process , molecular function , and cellular component  .
metabolic process, cellular process and response to stimulus rank among the first three categories among the biological processes.
the molecular functions most represented were binding, catalytic activity and transporter activity.
cell, cell part and organelle are the most represented categories in cellular components.figure  <dig>
go classification of
c. hybridum
unigenes.
the global transcript abundance could be deduced from a non-normalized cdna library.
the transcripts highly abundant in the cymbidium transcriptome are listed in additional file 6: table s <dig>  there were  <dig>  unigenes mapped into  <dig> kegg pathways.
of the  <dig>  unigenes,  <dig>   were related to metabolic pathways,  <dig>  to biosynthesis of secondary metabolites,  <dig>  to spliceosome,  <dig>  to plant-pathogen interaction,  <dig>  to protein processing in endoplasmic reticulum, and  <dig>  to endocytosis .
from the c. hybridum transcript database,  <dig> of the assembled sequences were assigned to  <dig>  cog annotations, which could be grouped into  <dig> categories .
the largest category was “general function prediction only” , followed by “transcription” , “replication, recombination and repair” , “signal transduction mechanisms” , “posttranslational modification, protein turnover, chaperones”  and “carbohydrate transport and metabolism” .figure  <dig>
cog function classification of
c. hybridum
unigenes.
validation of assembled transcripts
to experimentally confirm that the unigenes obtained from sequencing and computational analysis were indeed expressed,  <dig> unigenes were chosen for further rt-pcr and sanger-sequencing.
blastn results showed that identities of the selected reference assemblies with sanger-sequencing results were between  <dig> -100% , confirming the reliability of assembled reference transcriptome.
these unigenes were also included in the next qrt-pcr analysis.
global analysis of differential gene expression in the symbiotic interaction
having generated a c. hybridum root reference transcriptome, our next goal was to perform comparative gene expression profiles in orchid roots inoculated with different beneficial fungi.
a total of four independent un-normalised cdna libraries were sequenced using the illumina hiseq <dig> 100 bp single read module .
for each library, the total number of counts for each read was determined and all the reads were mapped back to the reference transcriptome .
in this study,  <dig>  million and  <dig>  million reads were uniquely mapped, representing between  <dig> % - <dig> % of total reads, thus, providing good coverage for the transcript profiles.
the expression level of each assembled transcript sequence in different samples was measured through rpkm values, and the significant changes in the transcript abundance were identified applying a q-value ≤  <dig>  and a two-fold change threshold.
based on these criteria, at  <dig> dpi,  <dig>   <dig> or <dig> genes were significantly regulated in c. hybridum root inoculated with isolates ml <dig>  zh3a- <dig> or mz, respectively.
the general differential expressed profiles in the three interactions are presented in figure  <dig> and additional file 10: tables s7-s <dig>  the top  <dig> abundant genes in symbiotic roots with positive fold change and top  <dig> abundant genes in non-symbiotic roots with negative fold change in response to symbiotic fungi were presented in additional file 11: tables s <dig> and s <dig>  om formation specifically induced 57% of the top  <dig> abundant genes in symbiotic roots with up-regulation pattern.
unigene2357_all , encoding metallothionein-like protein type  <dig>  which limited oxidative damage through binding to and detoxifying excess copper and other metals, was the most abundant gene in om symbiotic roots with more than  <dig> -fold up-regulation.
other om-induced genes in this category were largely involved in inhibition of reactive oxygen species  , protein metabolism and processing , signaling , cell wall modification  and cell defense and recue .
twenty-three genes in this category were co-induced by three treatments, involving in signaling , defense  and general metabolism .
however, 63% of the top  <dig> abundant genes in non-symbiotic roots with negative regulation pattern in response to inoculants were shared by three treatments and another 29% of these genes were specifically down-regulated by om fungi.
the most representative group abundant in non-symbiotic roots and down- regulated by three treatments was attributed to genes encoding main components in photosynthesis system, whereas a mannose-specific lectin, with  <dig>  rpkm in control roots, was specifically down-regulated by om fungi.figure  <dig>
venn diagram of differentially expressed genes in the
c. hybridum
- beneficial fungi interaction.
mock-inoculated , ml01-treated , zh3a-3-treated  or ml01 + zh3a- <dig> -treated .
as shown in the venn diagram,  <dig> degs were shared by ml01-, zh3a-3- and mz- treatments and  <dig> transcripts were collectively regulated in cyex <dig> and cyex <dig> but not in cyex <dig> which might specifically functioned in om mycorrhizal formation.
individual inoculation with non-mycorrhizal fungus zh3a- <dig> specifically regulated  <dig> genes.
another  <dig> degs with special expression pattern in cyex <dig> probably functioned in the combinative symbiotic process.
all the degs were grouped according to hierarchical cluster analysis  and twenty different expression sub-clusters based on k-means were identified .
the lists of all the degs in  <dig> sub-clusters were also presented in additional file  <dig>  the genes with the same regulation pattern in three treatments were assigned to subcluster  <dig>   <dig>   <dig>   <dig>   <dig>   <dig>   <dig>   <dig>   <dig>   <dig>   <dig>   <dig> and  <dig>  the genes in subcluster  <dig>   <dig>   <dig>   <dig> were specifically regulated by om formation.
genes drastically up-regulated by om fungi were included in subcluster  <dig>  the functions of these genes were related to signaling , reactive oxygen species homeostasis and stress-related , protein metabolism and processing  and transport .figure  <dig>
twenty regulatory patterns of all the degs expressed in different
c. hybridum
- beneficial fungi interactions based on k-means method.
go annotation was also used to compare the functional distribution of degs in different treatments .
among the represented biological processes, the most represented items were “ion transport” for om formation and “cellular biosynthetic process” for individual zh3a- <dig> treatment.
among the molecular functions, most represented degs were categorized as “oxidoreductase activity” for om formation and “cellulose synthase activity” for individual zh3a- <dig> treatment.
among the cellular component, the most represented degs were assigned to “extracellular region” for om formation and “plastid” for individual zh3a- <dig> treatment.
to identify the biological pathways regulated during orchid- beneficial fungi interaction, we assigned the degs to the reference kegg pathways .
the top  <dig> pathways with the most representation by unique sequences for each treatment were included in figure s7-s <dig> .
a large number of unigenes related to primary and secondary metabolism were collectively regulated by three treatments as 'nitrogen metabolism', 'flavonoid biosynthesis', 'cyanoamino acid metabolism', 'pentose phosphate pathway', 'starch and sucrose metabolism', 'alanine, aspartate and glutamate metabolism', 'stilbenoid, diarylheptanoid and gingerol biosynthesis', 'porphyrin and chlorophyll metabolism', 'glyoxylate and dicarboxylate metabolism', 'metabolic pathways' and 'biosynthesis of secondary metabolites'.
substantial numbers of genes/proteins involved in signaling pathway, protein turnover, nutrient transport, hormone action, cell rescue and defense were identified in om symbiosis specifically.
validation of deg profiling with qrt-pcr
using additional biological samples of ml01-, zh3a-3- and mz- colonized roots in comparison to mock-inoculated control roots, we performed real-time rt-pcr to verify the digital expression profiles of  <dig> genes corresponding to a range of functional categories and regulation patterns in the c. hybridum- beneficial fungi interactions.
the pearson correlation coefficient was calculated by spss to assess the correlation between different platforms.
overall, the qrt-pcr values were highly correlated with the rna-seq results, which confirmed the general expression pattern for all the selected genes in three treatments .table  <dig>
comparison of the expression profiles of  <dig> genes by rna-seq and quantitative real-time rt-pcr
fold changes of selected c. hybridum trasncripts  were logarithmic transformed which were determined by rna-seq and qrt-pcr in root tissues treated with isolates ml <dig>  , zh3a- <dig>  and ml01+ zh3a- <dig>  at  <dig> dpi as compared to mock-inoculated control .
base on the qrt-pcr results, we further analyzed the functional distribution of all the confirmed genes.
the function of co-regulated genes in three treatments included signaling , membrane proliferation , scavenging of reactive oxygen species , general stress- responsive , nutrient transport , protein processing , transcription factor  and hormone transport .
degs derived exclusively from om roots were assigned to signaling , transport , transcription factor , cell defense  and hormone metabolism .two genes showed special expression pattern in response to zh3a- <dig> treatment, including unigene 4239_all  and unigene 13219_all .
discussion
some different effects of mycorrhizal and non-mycorrhizal fungi on the growth of c. hybridumplantlets
this study demonstrated that the vegetative growth of c. hybridum plantlets could be improved by orchid mycorrhizal fungus ml <dig> and/or non-mycorrhizal isolate zh3a- <dig>  it is widely accepted that tulasnella-like rhizoctonia formed typical mycorrhizae with orchids and improved their germination and their seedlings’ growth  <cit> .
consistent with these studies, ml <dig> formed typical pelotons in host root cortex and significantly improved the total fresh weight  and dry biomass  of c. hybridum plantlets compared to the control.
however, we observed that the roots of ml01- treated plantlets did not elongate obviously as the co-cultivation proceeded, which might be due to rapid acidification of the co-cultivation media .
it was thus inferred ml <dig> increased the total dry biomass of orchid plantlets through improving the absorption of mineral nutrient.
non-mycorrhizal fungus zh3a- <dig> alone showed similar promoting effects on total fw and dw of cymbidium plantlets to that of mz treatment at  <dig> dpi , but we also observed that individual zh3a-3- treated orchids initiated less new roots than that of dual- inoculated plantlets, and their leaves turned into light green after a prolonged co-cultivation.
isolate zh3a- <dig> always co-existed with tulasnella-like rhizoctonia in wild orchid roots and the hyphae of isolate zh3a- <dig> mainly spread in velamen cells and could not form pelotons in the cortex cells of cymbidium roots.
it was also revealed that the effects of isolate zh3a- <dig> and another omf isolate  on hosts’ uptake of mineral nutrient were complementary in our recent report  <cit> .
thus, we designated zh3a- <dig> as mycorrhizal-associated fungi.
functional analyses of differentially regulated genes by om symbiosis
the establishment of om symbiosis involves the mutual recognition of host and fungi, hypha penetrating the cell walls of root hairs or epidermis and then entering the cells of the cortical parenchma, forming characteristic and complex hyphal coils   <cit> .
as in other endomycorrhizal systems, the plant cell membranes in the infected tissues remain intact at all times, excluding the hyphae from the cytoplasmic compartment, and an interfacial matrix is formed between the invading hyphae and the invaginating plant plasmalemma  <cit> .
the fungi further spread either as a result of repeated colonization or by hyphae penetrating from cell to cell in the root cortex, and the cortex cells may be uncolonized or contain active pelotons or clumps of degenerating hyphae in different proportions.
so, rasmussen et al.
<cit>  proposed that there co-existed necrotrophic and biotrophic phase in orchid mycorrhiza, and the former referred to peloton degradation and unilateral nutrient transfer by endocytosis, and the latter, mutual nutritional exchange  <cit> .
a recent labeling study in spiranthes sinensis clearly indicated that in the symbiotic protocorm, hyphal degeneration played a significant role in the mass transfer of the elements to plant cells, and fungal carbon and nitrogen was also detected to be transferred from live hyphae through the interface between the symbionts  <cit> .
owing to the non- synchronous process, the symbiotic root system often contains fungal structures in different developmental stages , and the functions of degs for om symbiosis were widely distributed.
genes related to signalingpathway
lysm-domain containing receptor-like kinases  were identified as receptors of acylated chitin  or chitin produced by plant interacting microbes and confirmed to be involved in signaling during nodule and am development  <cit> .
in this context, it is noteworthy that two genes  predicted to encode lysm domain receptor-like kinases showed evident up-regulation in om symbiosis.
one of them, unigene 1869_all, encodes a protein with 56% identity to lysm type receptor kinase lys <dig> from l. japonicus .
its induction was validated by qrt-pcr, which measured 16-fold up-regulation in the mycorrhizal roots.
this was consistent with the notion that symbiosis signaling seems to occur throughout the process of root colonization by a symbiont and occurs in parallel to additional signaling mechanisms that require the lysm receptor-like kinases  <cit> .
ca2+ spiking in the nucleus and perinuclear region of root hair cells has been documented as one of the earliest cellular responses after the perception of symbionts by host plants  <cit> .
our digital analysis also revealed many transcripts associated with ca2+ signals, as annexin and calcium-binding ef hand family proteins were co-regulated by three treatments.
but we also detected some transcripts in this group were specifically regulated by om fungi, including five genes encoding calcium-binding proteins , two genes encoding homologues of calcium-dependent protein kinase  , one cyclic nucleotide-gated ion channel protein , and one gene , encoding a homologue of dmi <dig> nodulation receptor kinase from m. truncatula.
another interesting finding was that two transcripts  encoding putative gras family protein with equal 44% identity to ljnsp <dig>  from l. japonicas, and one gene , which encoded gras family transcription factor with 56% identity to nsp <dig>  from m. truncatula, were also specifically regulated by om symbiosis.
nsp <dig> and nsp <dig> are required for rhizobium-specific gene expression whereas only nsp <dig> might be required for mycorrhizal signaling and nsp <dig> is replaced with required for arbuscular mycorrhization  <dig>  to drive expression of mycorrhiza-specific genes  <cit> .
genes involved in cellular organization and protein processing
a substantial number of genes related to cell wall modification were found to be co-regulated in three treatments, which encoded probable xyloglucan endotransglucosylase, beta-d-xylosidase or cellulose synthase.
two of these were confirmed by qrt-pcr.
unigene 8895_all, encoding a protein with 75% identity to the xyloglucan endotransglucosylase  of medicago trancatula, was co-induced more than 4-fold.
another gene  was strongly up-regulated more than 32-fold upon symbiosis .
the translated protein shares 92% identity to cellulose synthase a  catalytic subunit of vitis vinifera.
cesa proteins are part of the cellulose synthase complex in higher plants located in the plasma membrane  <cit> .
guether et al.
<cit>  also confirmed that two transcripts of lotus, ljcel <dig> and ljcesa, accumulated specifically in the arbusculated cells and functioned in cell expansion during arbuscule development  <cit> .
the strongly up-regulated transcripts of putative cellulose synthase might be involved in plant cell wall reinforcement during fungal colonization.
similar to am fungi, orchid mycorrhizal fungi penetrate the plant cell wall but not the plasma membrane.
as the colonization proceeds, the host plasma membrane invaginates and the infected cells often contain numerous mitochondria, ribosomes and profusely developed er.
a large number of sequences that were specifically induced by om formation were assigned to protein metabolism and turnover, membrane dynamics and cell wall synthesis.
this group of om specifically induced genes in cymbidium were thought to be homologues associated with am and ectomycorrhizal symbioses  <cit> .
besides the co-regulated degs among three treatments, om symbiosis also specifically induced a high number of transcripts encoding cell wall-degrading enzymes, including five glucan endo- <dig> -beta-glucosidase, one pectinesterase  <dig> and one xyloglucan endotransglucosylase/hydrolase protein.
concomitantly, three extensin- like proteins , and three putative glycosyltransferases  involved in cell wall reinforcement were also found in our dataset.
overall, this indicates an importance of wall- strength maintenance in om symbiosis.
with regard to membrane proliferation, two homologues of a. thaliana syntaxin syp  <dig>  unigene10103_all and unigene7281_all, were exclusively induced 16- and 6- fold by mycorrhizal symbiosis.
syntaxin belong to the superfamily of the snare proteins that play key roles in membrane fusion events during vesicle trafficking.
in am symbiosis, a gene  encoding a protein with a 60% identity to the a. thaliana syntaxin syp  <dig> also showed seven-fold up-regulation and was suggested to facilitate plant/fungal compatibility in mycorrhizal lotus roots  <cit> .
protein metabolism and processing seems to be a major process in mycorrhizal symbiosis  <cit> .
most drastically up om- regulated genes in the present study were involved in protein synthesis and processing .
concerning protein synthesis, the strong induction of different ribosomal proteins in response to om fungi provided evidence.
a number of transcripts associated with protein turnover were also detected in our analysis, including three subtilisin-like proteases , three serine carboxypeptidase , and also some peptidase .
some transcripts associated with correct protein folding  and degradation of mis-folded protein  were also identified as specifically up-regulated in om symbiosis.
at the same time, a transcript of putative serine-type endopeptidase inhibitor, unigene16319_all, accumulated more than 40-fold in the mycorrhizal roots compared with controls, which might fine-tune protein processing.
transport-related genes
phosphate transfer to the adult orchids via mycorrhizal fungi was well documented, just as in other mycorrhizal systems.
recently, labeling experiments also directly demonstrated that the fully autotrophic orchid goodyera repens acquires carbon, nitrogen and phosphorous from its fungal partner  <cit> .
notably, g. repens also transfers significant amounts of photosynthate  back to its ceratobasidium mycobiont – the first direct demonstration of a net carbon flow from orchid to fungus  <cit> .
tulasnella-like rhizoctonia were also confirmed to improve the plant growth of c. spp.
by enhancing nutrient acquisition  <cit> .
in our global profiling analysis,  <dig> genes associated with membrane transport were identified as co-regulated in three treatments.
among them,  <dig> genes involved in nutrient transport were co-induced in symbiotic roots, including two inorganic phosphate transporters , a plasma membrane atpase  <dig> , a bcs <dig> protein precursor-like  a sorbitol transporter , an oligopeptide transporter 1-like  and a magnesium transporter , whereas  <dig> genes were co-repressed at this time point which encoded two putative copper-transporting atpase, a copper transporter, a probable nitrite transporter, a sulfate transporter  <dig> , an inorganic phosphate transporter 2– <dig>  chloroplastic, a probable sodium-coupled neutral amino acid transporter  <dig> and a member of major facilitator family.
of these, the expression pattern of five genes was further verified by qrt-pcr.
the first of these, unigene6086_all, was co-induced more than 5-fold, which encodes an inorganic phosphate transporter with 81% similarity to inorganic phosphate transporter 1– <dig> of arabidopsis thaliana and a shared 78% identity to inorganic phosphate transporter from lotus japonicus  and m. truncatula .
both deduced amino acid sequences of two up-regulated inorganic phosphate transporter genes  in our digital profiling have similarity to reported mycorrhiza-specific or -induced phosphate transporters  <cit>  .
the second gene, unigene2386_all, which was up-regulated between 4- and 16- fold in the symbiotic process, encodes a protein with 63% identity to v. vinifera oligopeptide transporter 1-like.
another two co-induced genes  were involved in carbohydrate transmembrane transport and further identified by qrt-pcr .
unigene15955_all encodes a putative sorbitol transporter with 71% similarity to sorbitol transporter of zea mays and unigene4130_all encodes a transporter with 60% identity to sugar transporter erd6-like  <dig> of a. thaliana.
finally, a member of a major facilitator family  was confirmed to be repressed more than 4-fold in symbiotic roots and its translated sequence shared 68% similarity to a major facilitator family protein of a. thaliana.
in general, the regulation of a number of transporter genes supported the notion that extensive nutrient exchange was involved in symbiosis of the host and beneficial fungi.figure  <dig>
phylogenetic tree for the amino acid sequences of plant phosphate transporters  described as mycorrhiza-specific or -induced transporters, and two putative pits of
c. hybridum
in this study were highlighted in red color.
the dendrogram was generated by mega  <dig>  software using clustlw for the alignment and the neighbor-joining method for the construction of phylogeny.
the plant pi-transporters accession numbers could be obtained in  <cit> .
we also identified a substantial proportion of genes involved in transport which were specifically regulated in mycorrhizal roots.
these included three transcripts  encoding putative plant plasma membrane proton atpases.
the amino acid sequence of the encoded protein by unigene11510_all shares 77% identity with the am-specific mtha1in m. truncatula.
the regulation of another transcript  was validated by qrt-pcr  and its translated sequence shared 70% similarity to h+-atpase  <dig> from a. thaliana.
in om, like in am and ecm, nutrients from symbiotic hyphae separated at the apoplastic space have to go through several membrane barriers before being assimilated by the partner's cells, so the activated plant plasma membrane proton atpase might be a sign of accelerated membrane transport.
the other two genes  were also identified as om-specific which encoded a putative copper-exporting atpase and a plasma membrane calcium-transporting atpase, respectively.
another group of mycorrhizal-induced genes identified was associated with nitrogen transport.
eight genes of peptide transporter  were found to be exclusively up-regulated between two-fold and more than 300-fold compared with controls.
the strongest up-regulated ptr , encoding a homologue of peptide transporter ptr <dig> of m. truncatula with 62% identity, showed none expression level in the control and cyex <dig> samples.
the second putative ptr  was 300-fold up regulated in mycorrhizal roots at  <dig> dpi and had a 76% similarity to a member of oligopeptide transporter family  from p. trichocarpa.
quantitative rt-pcr analysis confirmed that one transcript  was up-regulated more than 7-fold in om symbiosis whereas down regulated two-fold in zh3a-3-treated orchid roots.
besides these ptrs, two putative amino acid transporters were also detected to be induced by mycorrhizal formation.
a homology search revealed 70% identity between the translated amino acid sequence of unigene3142_all and cationic amino acid transporter  <dig>  of a. thaliana.
the second gene, unigene404_all, encodes a protein with 81% identity to lht <dig>  of a. thaliana , which was regulated by diverse factors.
our data supported the notion that omf played an important role in improving n acquisition of host plants  <cit> .
two putative monosaccharide transporter genes  were also found to be induced in om mycorrhizal roots.
the closest homologue of translated sequence by unigene12841_all in a. thaliana is sugar-proton symporter plt <dig> , with 71% identity at amino acid level.
another gene, unigene4472_all, encoded a homologue with 70% similarity to a tonoplast monosaccharide transporter  <dig>  in a. thaliana and tmt1/ <dig> probably represents a family of proton-coupled antiporters capable of high-capacity loading of glucose and sucrose into the vacuole.
gaude et al.
<cit>  also demonstrated that two putative sucrose transporters and one hexose transporter were induced in non-colonized cortical cells of mycorrhizal roots compared to the cortical cells of non-mycorrhizal roots, and suggested they might be involved in the uptake of carbohydrate from the apoplast into cells in the vicinity of hyphae in arbuscular mycorrhizal roots or functioned in the proton-coupled uptake of carbohydrate into the cytoplasm from the vacuole  <cit> .
this raises a compelling question of what is the exact function of these two putative monosaccharide transporters in our study.
finally, a gene encoding a putative potassium transporter  was also detected to be 6-fold up-regulated in mycorrhizal roots and shared 50% identity at the amino acid level with the potassium transporter  <dig> of a. thaliana.
in addition to a number of up-regulated genes probably involved in nutrient transport, atp binding cassette  transporter genes, which are unlikely associated with direct nutrient transport, were also identified to be regulated in mycorrhizal roots.
the strongest up-regulated gene of abc transporter family  encoded a protein with 66% identity to pdr <dig> abc transporter of oryza sativa, which was regulated by many abiotic and biotic factors.
recently, two plant half-size atp binding cassette  transporter proteins  were identified that are essential for am symbiosis  <cit>  and one full-size abc transporter gene mtabcb <dig>  showed strongly increased transcript levels in arbuscule-containing  cells and also in adjacent cells  <cit> .
these abc transporter proteins may be involved in the transport of signaling between plant cells during am symbiosis  <cit> .
here, these identified abc transporter genes showing different expression patterns in orchid mycorrhizal roots may be important for om symbiosis.
om-activated transcripts also included components associated with vesicle-mediated transport, which might involve in the membrane biogenesis for an intracellular accommodation of fungal structures.
the translated sequence of unigene7513_all  shares 70% identity with vamp7b  in a. lyrata subsp.
lyrata.
another gene, unigene11358_all, encoding a putative clathrin assembly protein, was up-regulated more than 8-fold in mycorrhizal roots.
defense -related and phytohormone regulation
oxidative burst belongs to early events of plant responses to compatible and incompatible fungi.
the success of beneficial microbes to colonize plant roots depends on their ability to manipulate immune responses of the host.
this was verified by the findings of transcripts associated with scavenging of reactive oxygen species , such as glutathione peroxidase, probable glutathione s-transferase  and germin-like protein , accumulating in the symbiotic roots.
one of these, unigene13078_all, encodes a protein with 77% identity to the glutathione peroxidase of populus trichocarpa and was further validated by qrt-pcr.
concomitantly, many transcripts encoding general 'stress- responsive' proteins ; bibenzyl synthase  were identified to be induced in these symbiotic process.
the enzymes in the pathway towards fungicidal phytoalexin synthesis are produced as a general feature in orchids.
reinecke  <cit>  have reported that the activity of benzyl synthase, catalyzing the biosynthesis of bibenzyls and  <dig>  10-dihydrophenanthrenes of bletilla striata rhizomes, increased substantially upon fungal infection  <cit> .
another interesting finding was that two transcripts encoding programmed cell death protein  <dig>  were down regulated collectively in three treatments.
pdcd <dig>  a suppressor of gene transcription and translation, plays a crucial inhibitory role in several types of human tumors through inhibiting autophagy in multiple cell types both in vitro and in vivo
<cit> .
down- regulation of these homologues of pdcd <dig> may promote the symbiotic process through suppressing plant defense.
furthermore, four transcripts , encoding homologues to key molecular regulators in two-component systems  which play important roles in responses to environmental stress stimuli  <cit> , were found to be down regulated regardless of the inoculated fungi species, highlighting the shared strategy by different beneficial fungi in the symbiosis.
apart from these shared degs, three chitinase genes  were exclusively induced by om fungi and the specific expression pattern of unigene 15958_all was confirmed by qrt-pcr.
the translated sequence of unigene 15958_all shared 59% identity with the chitinase 4-like of fragaria vesca.
chitinases hydrolyze β- <dig> -glycosidic bonds and can be induced in pathogenic interactions.
but the expression of unigene 15958_all was not induced by a pathogenic fusarium isolate .
unigene8905_all encoded a protein with 55% identity to class iii chitinase of m. truncatula and class iii chitinases were found to be expressed predominantly in the arbuscule-containing cells  <cit> , where they may suppress plant defense by the reduction of chitin-like elicitors during the formation of functional symbiotic interfaces  <cit> .
the similar role of class iii chitinase svchit <dig> in om was recently suggested by perotto  <cit> .
thus, these up-regulated chitinases are probably involved in the lysis of fungal cell wall and balancing the symbiosis.
another notable om-regulated group was mannose-specific lectins .
among these, unigene4318_all and unigene529_all showed high-level expression levels in mycorrhizal roots with little expression in the un-inoculated control and zh3a- <dig> treatments.
antifungal activity against a range of phytopathogenic fungi has been demonstrated for some of mannose-specific lectins both in vitro and in transgenic plants, and they were thus proposed to be part of defense mechanisms against the invading fungus.
the strong up-regulation of svlect <dig> and svlect <dig> in response to mycorrhizal fungi t. calospora in s. vomeracea protocorm was recently confirmed by qrt-pcr and this was suggested to be part of a mechanism to limit fungal growth in the mycorrhizal protocorm  <cit> .
however, another mannose-specific lectin gene, unigene14222_all, with  <dig>  rpkm in control roots, was specifically down-regulated by om fungi.
owing to the multi-functions of plant lectins, the role of these lectins in the orchid mycorrhizal process remains to be established.
among plant hormones, ethylene , salicylic acid , abscisic acid , and jasmonic acid  are known to be key elements in fine-tuning the plant defense response during interaction with other organisms  <cit> .
the inhibitory effect of exogenous ethylene in am development has been documented and the colonization of host tissues by mutualistic am fungi such as rhizophagus irregularis , ectomycorrhizal fungus laccaria bicolor, or endophytic fungus piriformospora indica was inhibited by et  <cit> .
however, et signaling is still essential for nodule formation in sesbania rostrata and in the initial stages of p. indica colonization  <cit> .
because ethylene-responsive transcription factors  can function as transcriptional activators and repressors, they are candidates for establishing a balanced defense response to the fungus without preventing growth and development in am symbiosis  <cit> .
valadares  <cit>  also suggested the role of ethylene in om regulation when compared the gene expression pattern of o. maculata mycorrhizal roots in comparison to non-mycorrhizal controls  <cit> .
in our global expression profiling, eight transcription factors  related to ethylene signaling  were found to accumulate in all the symbiotic roots and four genes among these were specifically induced in om roots.
among these, the induction of two transcripts  was further validated by by qrt-pcr.
the translated sequence of unigene1724_all showed 63% similarity to the ethylene receptor homolog  <dig>  of zea mays with two-fold induction.
unigene12197_all, encoding a protein with 45% similarity to ethylene-responsive element binding factor  <dig>  of g. max, was induced more than 5- fold in the symbiotic roots.
as with plett et al.
<cit> , the induction of these tfs is probably an adaptive response by the plant so that its growth and vigor are not compromised by the fungus in the late stage of symbiosis  <cit> .table  <dig>
transcription factors among beneficial fungi- induced genes
lists of tfs in this table included specifically regulated tfs in om formation and * co-regulated genes in three treatments; ⊕qpcr validation; number in parenthesis referred to the gene number regulated by om formation.
aba plays a crucial role in biotic and abiotic stress responses.
analysis of am colonization in the sitiens aba-deficient tomato mutant showed that aba was necessary in order to complete arbuscule formation and its functionality and to promote sustained colonization of the plant root  <cit> .
subsequent work supported the notion that aba deficiency negatively regulates mycorrhizal intensity/hyphal colonization directly or indirectly through increasing ethylene levels .
accumulation of the transcripts related to aba de novo biosynthesis was documented in am roots  <cit>  and om green protocorms  <cit> .
at the same time, genes associated with aba catabolism also showed increased transcript levels in am roots <cit>  and om roots  <cit> .
in our study,  <dig> transcripts related to aba biosynthesis: 9-cis-epoxycarotenoid dioxygenase   and zeaxanthin epoxidase  , were co-suppressed in symbiotic roots.
aba deficiency might be involved in restriction of fungal spread in orchid symbiotic roots.
jasmonic acid has long been implicated in the plant’s systemic response to pathogen attack  <cit> .
while experimental results sometimes appear contradictory, most authors agree that jasmonate signalling is important for am colonization and development  <cit> .
plett  <cit>  revealed that in the interaction between l. bicolor and poplar, ethylene and jasmonic acid showed extensive transcriptional cross-talk and functioned as negative modulators during mutualistic symbiosis  <cit> .
jasmonates belong to a diverse class of lipid metabolites known as oxylipins and the key genes in the biosynthesis of oxylipins were commonly induced in mycorrhizal roots  <cit> .
consistent with these studies, the induction of  <dig> lipoxygenase  genes in our study suggest the involvement of the oxylipin pathway in regulation of om symbiosis.
in plants, there are two main branches of the oxylipin pathway, determined by two different types of lipoxygenases , 9-lox and 13-loxs.
the 13-lox pathway leads to the biosynthesis of ja and derivatives, and the relevance of the 9-lox pathway in plant interactions with nematodes and pathogens has also been demonstrated recently  <cit> .
further research is required to elucidate the exact role of these up-regulated loxs in long-term balance of the om symbiosis.
apart from the potential roles for et, aba and ja in om symbiosis, three transcripts associated with gibberellin  biosynthesis were also increased in abundance in orchid mycorrhizal roots: two genes  encoding gibberellin 3-oxidase  and another  encoding ga 20-oxidase .
the translated sequence of unigene13170_all shares 77% identity to gibberellin 20-oxidase of triticum aestivum and its induction was confirmed by qrt-pcr.
the product of ga3ox and ga20ox activates ga, suggesting that ga synthesis is essential for om formation.
in barley- p. indica interaction, schäfer  <cit>  also observed that two barley mutants impaired in ga synthesis or perception, reduced colonization by p. indica.
similarly, homologues of genes involved in the ga biosynthesis were already identified as mycorrhizal-responsive in tomato plants at the mature  phase of the interaction  <cit> .
foo  <cit>  reported that biologically active gibberellins suppress arbuscule formation in pea roots, and gibberellin signaling della proteins are essential for this response.
owing to extensive transcriptional cross-talk between ga signaling and sa/ja responses, the ga-regulated defense response might be critical for establishing these long-term mutualisms.
transcription factors
as detailed above, three treatments initiated common stress responses for symbiotic roots, so the majority of co-regulated tfs encoded ap <dig> domain tfs and ethylene-responsive tfs , wrky tfs, myb- family tfs, basic helix-loop-helix  tfs and zinc-finger tfs .
dna-binding wrky, ap2-domain and erf tfs are generally involved during pathogen-related stress responses  <cit> .
the co-regulated tfs may also be involved in root developmental process, as are myb family, erf and nac family  <cit> .
we also detected gras family and arf tfs specifically induced by om formation.
four transcripts  encoding for gras members were putative scarecrow  protein encoding genes, considered to be genetic determinants of root identity.
previous studies also confirmed that members of the gras family are essential for nodule development  <cit>  and they may be important for the regulation of gene expression in am roots of m. truncatula, l. japonicus and s. lycopersicum
<cit> .
two putative arf transcripts were also identified to be up- regulated in orchid mycorrhizal roots.
the induction of putative arfs has also been found during am symbiosis in maize, rice and m. truncatula, but not in l. japonicus
<cit> .
and in am interaction, mycorrhizal roots exhibit morphological characteristics such as an increase in the number of lateral/fine roots during early growth phases, similar to auxin-treated roots  <cit> .
recently, sukumar  and fusconi  <cit>  summarized a compelling supportive role for involvement of auxin pathways in beneficial plant–microbe interactions  <cit> .
chutima  <cit>  reported that orchid-associated fungi tulasnella sp.
produced high levels of indole-3- acetic acid  with all biological activities in culture medium supplemented with 2 mg/ml of l-tryptophan  <cit> .
the iaa production by isolates ml <dig> was also identified in our research group .furthermore, unigene11039_all, predicted to translate a protein with 66% similarity to an auxin efflux facilitator slpin <dig> of solanum lycopersicum, showed 15-fold up-regulation in the symbiotic roots.
all these data point to the probable involvement of auxin signaling pathway in om functionality.
degs specifically regulated by non-mycorrhizal fungus zh3a-3
the hyphae of non-mycorrhizal fungus zh3a- <dig> colonized the velamen cells of orchid roots and individual inoculation significantly improved the root length and fw of c. hybridum plantlets.
as detailed above,  <dig> genes were significantly regulated in c. hybridum roots inoculated with zh3a- <dig> alone and most of these degs were shared with those of ml <dig> and/or mz treatment.
apart from these shared degs, we also identified  <dig> genes were specifically regulated by zh3a- <dig> treatment.
the largest number of up-regulated transcripts by zh3a- <dig> treatment was associated with orchid root development.
these genes are mainly related to hormonal balance and transcription regulation.
among plant phytohormones, auxin, cytokinin, and ethylene play important role for root development.
consistent with this, we detected the up-regulation of many corresponding genes, which encoded one cytokinin dehydrogenase , two erf , and three abc transporter b family members .
auxin plays a pivotal role in many plant- beneficial microbe interactions  <cit> , and it is well established that higher auxin/cytokinin ratio determines the root initiation.
a recent study also clearly demonstrated that osckx <dig>  a cytokinin oxidase/dehydrogenase  family gene, integrates cytokinin and auxin signaling to control rice crown root formation positively  <cit> .
so, the up-regulation of unigene14465_all, which encodes a homolog with 69% identity to ckx of a. thaliana, may suggest its positive role in orchid- zh3a- <dig> interaction.
polar auxin transport  is essential for establishment of auxin gradient that spatially defines the quiescent center  and consequently the fate of neighboring stem  <cit> .
abc transporter b family members were probably involved in auxin acropetal or basipetal transport  <cit> , thus, the up-regulation of three abc transporter b family members appears to converge on the auxin pathway to modulate root development.
as shown in table  <dig>  the putative myb-tfs ranked the third largest number of the regulated tfs in the mycorrhizal symbiosis, and most were shared by three treatments, which suggested their essential role in response to symbiotic fungi.
recently, a convincing research supported that a putative myb-like tf of l. japonicus  which was am-responsive may also have non-symbiotic functions: i.e.
root growth  <cit> .
in our study, two transcripts ,encoding homologues of myb family tfs of o. sativa japonica group and a. thaliana with 60%- 80% identity, were also detected to be up-regulated 4-fold by zh3a- <dig> the translated sequence of unigene2947_all also shared 46% similarity to ljmami at amino acid level.
another up-regulated tf related to root development was unigene489_all , which encodes a protein with 45% similarity to scarecrow-like protein of z. mays.
we also identified a substantial number of genes related to general disease resistance processes which were specifically up-regulated by zh3a- <dig> treatment.
these included  <dig> transcripts related to signaling ,  <dig> transcripts associated with balance of reactive oxygen , and  <dig> transcripts involved in secondary metabolism .
the increased activity of chitinases were generally found in nonspecific, broad-spectrum defenses, but we did not detect the transcript accumulation of chitinase genes in cyex <dig> the fungal cell wall of isolate zh3a- <dig> always ruptured as the cultivation proceeded , which might be the reason.
the activation of basal defense mechanisms of orchid plants by beneficial fungus zh3a- <dig> did not affect the symbiotic effects, as described in  <cit> .
as zh3a- <dig> colonization induced dramatic structural changes in root cells and promoted the root growth of orchids, transcripts accumulated in the symbiotic roots were also related to carbohydrate and amino acid metabolism, protein synthesis and processing, and structural re-organization.
furthermore, a transcript  involved in water transport was also up-regulated, which encodes a homolog of aquaporin pip1- <dig> of z. mays with 87% identity.
taken together, we confirmed the findings which suggested that the development and metabolism of plant symbionts are generally driven by differential regulation of transcriptional regulators, signal transduction, and metabolic pathways, rather than by expression of symbiosis-specific genes  <cit> .
detection of putative fungal genes in symbiotic roots
owing to the distinct phylogenetic relationship between isolate ml <dig> and zh3a- <dig>  we could detect  <dig> possible fungal genes combining the annotation information of de novo assembled transcripts with the comparative expressed level in different symbiotic roots .
of these  <dig> genes,  <dig> were probably genes from isolate zh3a- <dig>  and the rest, from ml <dig>  the origin of putative ml <dig> genes was further confirmed by blastn search against the t. calospora draft genome or t. calospora mycelium ests with an e value > 1e- <dig> .
transcripts from ml <dig> in the mycorrhizal roots were largely assigned to ribosomal proteins, enzymes that might be involved in plant cell wall degradation  and remodeling the fungal cell wall during growth and symbiosis .
we also identified numerous transcripts encoding different proteins associated with the nutrient transport and assimilation, such as h+-transporting atpase, phosphate transporter , ammonium transporter , amino acid and protein transporters, zip-like iron-zinc transporter, monosaccharide transporters  and mfs-type transporters .
abc-type fe3+-siderophore transporter and other transporters most probably implicated in the detoxification process were also detected in the mycorrhizal roots.
the translated amino acid sequence of unigene5906_all shared 59-60% identity to high-affinity pi- transporters from g. versiforme  and r. intraradices  .
the function of rinpt  has been suggested by tisserant et al.
<cit>  to be involved in stabilizing cooperation in the mycorrhizal symbiosis  <cit> .
the transcript abundance of unigene5906_all in cyex <dig> was higher than that in cyex <dig>  so, this gene might be responsible for the better growth- promoting effect of mz treatment than that of ml <dig> alone.
besides the transcripts in n and p assimilation pathway, iron-zinc transporter and abc-type fe3+-siderophore transporter co-existed in the mycorrhizal roots, which suggest an important role of some metals in om functionality.
our previous study showed that isolate cl <dig>  tulasnella-like rhizoctonia, formed typical mycorrhizae with orchids and enhanced the n, p, zn, fe contents of mycorrhzial plants  <cit> .
the its sequence of ml <dig> shared 84% similarity to that of isolate cl <dig> and both functioned similarly in improving the growth and nutritional absorption of orchid plantlets .
however, there is a need to compare the differences in gene expression between the intraradical and extraradical mycelium so as to elucidate the role of these transcripts in om symbiosis.
CONCLUSIONS
using a combination of rna-seq and real-time rt-pcr, we characterized the root transcriptomic reprogramming of c. hybridum with respect to different beneficial fungi and identified specific changes in relative abundance for transcripts of genes in om symbiosis.
genes involved in cell wall modification, reactive oxygen species detoxification, defense-related phytohormone and phosphate transport were co- induced in all the symbiotic interactions.
the transcripts increased specifically by om were related to signaling , cell wall degradation or reinforcement, protein metabolism and processing, defense  and n, fe transport.
among these degs, genes potentially functioned in c, n, p and fe transport, hormone metabolism and signaling are of particular interest.
aside from these orchid transcripts, we also identified some putative fungal genes in symbiotic roots which were associated with plant cell wall degradation, remodeling the fungal cell wall and nutrient transport.
this comprehensive dataset provides a basis for future research in exploring the molecular mechanisms underlying om formation and functionality.
electronic supplementary material
additional file 1: table s1: identities of the selected reference assemblies with sanger-sequencing results of corresponding touch-down pcr products.
additional file 2: table s2: gene-specific primers used for qrt-pcr.
additional file 3: figure s1: typical co-cultivation condition of cymbidium hybridum plantlets with different fungi at  <dig> dpi.
mock-inoculated control  and inoculated with isolate ml <dig> , zh3a- <dig>  and ml01 + zh3a- <dig> .
additional file 4: figure s2: root staining results of each treatment at different symbiotic time points.
roots inoculated with e. repens isolate ml <dig> at 6 days post-inoculation ;  roots inoculated with u. nana isolate zh3a- <dig> at  <dig> dpi;  roots inoculated with isolates ml <dig> and zh3a- <dig> at  <dig> dpi;  mock-inoculated roots at  <dig> dpi;  roots inoculated with e. repens isolate ml <dig> at  <dig> dpi;  roots inoculated with u. nana isolate zh3a- <dig> at  <dig> dpi;  roots inoculated with isolates ml <dig> and zh3a- <dig> at  <dig> dpi;  mock-inoculated roots at  <dig> dpi;  roots inoculated with e. repens isolate ml <dig> at  <dig> dpi;  roots inoculated with u. nana isolate zh3a- <dig> at  <dig> dpi;  roots inoculated with isolates ml <dig> and zh3a- <dig> at  <dig> dpi;  mock-inoculated roots at  <dig> dpi;  roots inoculated with e. repens isolate ml <dig> at  <dig> dpi;  roots inoculated with u. nana isolate zh3a- <dig> at  <dig> dpi;  roots inoculated with isolates ml <dig> and zh3a- <dig> at  <dig> dpi;  mock-inoculated roots at  <dig> dpi.
additional file 5:
cymbidium hybridum
root transcriptome assembly- illmina 90 bp reads.
additional file 6: table s3: highly abundant transcripts in cytr-unigenes .
additional file 7: table s4: pathway annotation of cytr-unigenes.
additional file 8: table s5: statistics of dge sequencing.
additional file 9: table s6: statistics of uniquely mapped reads for each sample.
additional file 10: table s7-s9: lists of differentially expressed genes in different c. hybridum- beneficial fungi interactions.
additional file 11: table s10: top  <dig> abundant genes in symbiotic roots with positive fold change in response to symbiotic fungi.
table s <dig>  top  <dig> abundant genes in non-symbiotic roots with negative fold change in response to symbiotic fungi.
additional file 12: figure s3: hierarchical clustering of degs from c. hybridum expressed in the roots in response to isolates ml <dig> , zh3a- <dig> , ml01 + zh3a- <dig> or mock inoculation .
additional file 13:
lists of degs in  <dig> sub-clusters based on k-means.
additional file 14: figure s4-s6: go enrichment bar graphs of degs in the interactions between c. hybridum and ml <dig> , zh3a- <dig>  or ml01 + zh3a- <dig> .
additional file 15: table s12-s14: lists of enriched pathways in different c. hybridum- beneficial fungi interactions.
additional file 16: figure s7-s9: scatterplot for top  <dig> enriched kegg pathways of degs in the interactions between c. hybridum and ml <dig> , zh3a- <dig>  or ml01 + zh3a- <dig> .
additional file 17: figure s10: comparison of expression profiles of random selected  <dig> genes by rna-seq and qrt-pcr.
additional file 18: table s15: expression data for putative genes from e. repens isolate ml <dig> in symbiotic roots.
table s <dig>  expression data for putative genes from u. nana isolate zh3a- <dig> in symbiotic roots.
additional file 19: figure s11: phylogenetic tree for the amino acid sequences of symbiosis- associated phosphate transporters from plant and fungi.the dendrogram was generated by mega  <dig>  software using clustlw for the alignment and the neighbor-joining method for the construction of phylogeny.
the plant pi-transporters accession numbers could be obtained in  <cit> .
fungal phosphate transporters accession numbers: glomus versiforme: gvept, aac <dig> ; gvept <dig>  gi|74654973; rhizophagus intraradices: rinpt, aal <dig> ; yeast: yepho <dig>  gi|1346710; sesbania rostrata: serpt, cac <dig> ; piriformospora indica: pipt <dig>  abi <dig> ; pipt <dig>  cca <dig> ; pipt <dig>  cca <dig> ; pipt <dig>  cca <dig> .
the nucleotide sequences of putative phosphate transporters in our study were available in additional file 6: table s <dig>
abbreviations
abcatp-binding cassette
amsarbuscular mycorrhizas
ap2/erfapetala2/et response factor
arfauxin response factor
bhlhbasic helix-loop-helix
cat5cationic amino acid transporter 5
cdpkcalcium-dependent protein kinase
ckxcytokinin oxidase/dehydrogenase
degsdifferentially expressed genes
dgedigital gene expression
ecmsectomycorrhizas
etrethylene receptor homolog
lht1lysine histidine transporter 1
lysm-rlklysm-domain containing receptor-like kinases
mafmycorrhizal- associated fungi
mdmycorrhizal dependency
mhheterotrophy
mxmixotrophs
naataanicotianamine aminotransferase a
ngsnext generation sequencing
nsp2nodulation-signaling pathway  <dig> protein
omorchid mycorrhizae
omforchid mycorrhizal fungi
optoligopeptide transporter family
patpolar auxin transport
pdapotato-dextrose agar
pdcdprogrammed cell death protein 4
ptphosphate transporter
ptrpeptide transporter
prgpercentage of root growth
qcquiescent center
ram1protein required for arbuscular mycorrhization 1
rmbtreads that could be mapped back to transcripts
rpkmreads per kilo base per million reads
scrscarecrow
spasuppressor of phya-105
tfstranscription factors
tmt2tonoplast monosaccharide transporter 2
vamp7bvesicle-associated membrane protein 7b
xetxyloglucan endotransglucosylase.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
zxl and zjx carried out the experiments, performed the bioinformatics analyses, and drafted the manuscript.
ccl cultured the experimental material and extracted all the rna samples.
zhy performed the symbiotic experiments.
yjz and lm participated in the qrt-pcr experiment.
lfb designed the study and revised the manuscript.
all authors read and approved the final manuscript.
we are grateful to xiaoyong liu  for his help with the identification of fungal isolates in our research.
we thank professor ronald r. sederoff  for his critical reading and revision of the manuscript.
we also thank for the constructive comments of the anonymous reviewers to improve our manuscript.
this research was financially supported by national natural science foundation of p. r. china , guangdong province science and technology project , and the cooperation project in industry, education and research of guangdong province and ministry of education of p. r. china , and the science and technology fund from south china agricultural university .

 all orchids produce minute, endosperm- lacking seeds and are dependent on fungal colonization for germination and growth into an underground heterotrophic, achlorophyllous stage called a protocorm under natural conditions. after this fully mycotrophic stage, most species shift to autotrophy at adult stage with the development of photosynthetic organs, but approximately species worldwide stay fully mycoheterotrophic. recently, some chlorophyllous species in orchidaceae were verified to obtain carbon at adult stage not only from their photosynthetic activity, but also from mycorrhizal fungi, which were called mixotrophs. among these, several terrestrial orchids such as c. lancifolium and c. goeringii were shown to be mx nutritional mode through isotopic analyses. the dependency of green orchids on mycorrhizal fungi throughout their life cycle prompts us to study these mycorrhizal interactions in more detail. in nature, orchids largely associate with members of the ‘rhizoctonia’ complex, which contains three now taxonomically disparate agaricomycetes taxa: sebacinales, ceratobasidiaceae and tulasnellaceae. tulasnellaceae are the most frequently found ‘rhizoctonias’ in green orchids. besides orchid mycorrhizal fungi, recent studies also demonstrated that isolates of umbelopsis, as non-mycorrhizal beneficial endophytes, probably co-existed in the roots of wild orchids. c. hybridum collectively contains a large number of cultivars with large flowers. over species in cymbidium, including c. sinense, c. goeringii, c. faberi, c. ensifolium and c. kanran, were used in their breeding programs. in our previous work, one omf isolate ml and one non-mycorrhizal beneficial isolate zh3a- were isolated from roots of wild c. spp. adult plants and identified as e. repens and u. nana, respectively. here, we present that ml forms typical mycorrhiza with c. hybridum. individual or combinative inoculation significantly improved the vegetative growth of c. hybridum plantlets compared with mock-inoculated controls. although orchid mycorrhizal symbiosis is an important feature in the life cycle of orchids, it has long been considered to be an atypical mycorrhizal association, with the fungus deriving little benefit from the orchid host. but cameron et al. reported that approximately – % of the labeled carbon provided to the orchids was passed to the fungal partner. thus adult orchid mycorrhizas may represent a truly mutualistic interaction similar to ectomycorrhizae and arbuscular mycorrhizae associations. compared with the numerous transcriptional analyses to identify genes involved in am and ecm symbiosis, the molecular mechanisms behind om symbiosis are still poorly understood. recently, different gene expression or proteome changes of orchid protocorms in response to mycorrhizal fungi were studied and mutualistic plant-fungus relationship was suggested. as for orchids at the adult stage, valadares used rna-seq to determine the gene expression profile of wild oeceoclades maculata mycorrhizal roots in comparison to non-mycorrhizal controls. while this study provides new insights into our understanding of the om symbiosis for adult orchids, the symbiosis seems to be highly complex as wild orchids generally associate with multiple symbionts even within a single root. although genomic sequence resources currently available for orchids are limited, large expressed transcripts from phalaenopsis aphrodite, p. equestris, p. bellina, oncidium gower ramsey, erycina pusilla and c. sinense were generated in recent years with the rapid advances in dna sequencing technology known as next-generation sequencing and the development of bioinformatic processes, especially for de novo assembly without reference genome information. whereas no genetic information on orchid mycorrhizal fungi was available when our analyses began, large numbers of assembled transcripts and the genome information of t. calospora had been recently generated as part of a doe jgi community sequencing program. this work has allowed us to assign some putative fungal genes to mycorrhizal fungi ml by blastn searches of these data. in this paper, we used ngs approach and illumina digital gene expression technology to obtain genome-wide insight into the transcriptomic responses of c. hybridum under specific interactions with mycorrhizal or/and other beneficial fungi. to our knowledge, this represents the first genome-wide gene expression profile of orchid plants in response to mycorrhizal or non-mycorrhizal fungi. the data enabled us to get new insight into the mechanisms underlying plant responses to mycorrhizal and non-mycorrhizal beneficial fungi. the transcriptome for c. hybridum could also be used as an important resource for exploring other biological processes in orchids. methods biological materials and inoculation the micro-propagated plantlets of c. hybridum ‘golden boy’ were provided by the guangdong key lab of ornamental plant germplasm innovation and utilization, china. the propagation procedure for orchid plantlets was described in. the two fungal isolates ml and zh3a- were preserved in the guangdong key laboratory for innovative development and utilization of forest plant germplasm. inocula were produced by growing each strain on 90 mm petri dishes containing potato-dextrose agar for 10 days. the plantlets with 3– roots about -3 cm long were then transferred to de basal medium, an in vitro mycorrhizal induction medium in orchids. the co-cultivation de medium contained  mm cacl  mm mgs  mm k2s  mm kh2p 100 μm fes 25 μm h3bo 33 μm mncl  μm znso  μm namoo 140 μm na2edta, supplemented with yeast extract 1 g · l- soluble starch 9 g · l- and agar 6 g · l- each plantlet was inoculated with a 7 mm plug excised from an edge of an actively growing colony on pda of strain ml zh3a- or ml01+ zh3a- the control plantlets were mock-inoculated with a plug excised from a pda plate without fungus. all the treated and control plantlets were then incubated under a photoperiod of 16 hrs light at a temperature of 22°c. symbiotic effects on vegetative growth of c. hybridum after 45 days post- inoculation, ninety plantlets of each treatment were totally harvested and randomly divided into three groups, which represented three biological replicates. the symbiotic effects on the vegetative growth of the hosts were assessed by total increased fresh weight and dry biomass of the plantlets. the percentage of root growth and mycorrhizal dependency were also compared among different treatments. md was calculated by the following formula:  × 100/shoot dry wt. of mycorrhizal plants as described in. for statistical analysis of data, spss v was employed for one-way analysis of variance and comparisons among means were made using duncan’s multiple range test, calculated at p < . root staining and histology the successful inoculation with strains ml and zh3a- and the absence of fungal contaminants in mock-treated roots were confirmed by trypan blue staining. to confirm their respective colonization pattern, semi-thin sections of root samples for each treatment were prepared and stained with a % aqueous solution of toluidine blue for light microscopy. reference cdna library construction at and dpi, orchid roots of mz- treated or mock-inoculated control plantlets were harvested by carefully removing the substrate and mycelia, and the samples were frozen in liquid nitrogen and stored at -80°c. total rna from each sample was extracted separately with column plant total rnaout kit according to the manufacturer’s protocol. the rna was treated with dnase and then quantified using a nanodrop 2000™ spectrophotometer. equal amounts of total rna from mz- treated and mock-inoculated samples at each time point were pooled to produce minimum 20 μg total rna for constructing a cdna library for establishing the reference transcriptome. mrna-sequencing and cdna de novoassembly the library construction and sequencing were performed by the beijing genomics institute and the general experimental pipeline was summarized by previously published papers. first, the cdna library was pair-end sequenced by illumina hiseq. after removing adaptor sequences, empty reads and low quality sequences, we obtained the clean reads for de novo assembly. three de novo assembly software were used including soapdenovo with optimized k-mer length of abyss with k-mer length of and trinity with default k-mer length of sequence directions of the resulting unigenes were validated by blastx alignment between unigenes and public protein databases with a priority order of nr, swiss-prot, kyoto encyclopedia of genes and genomes database, and cog. when a unigene happens to be unaligned to none of the above databases, estscan will be introduced to predict its coding regions as well as to decide its sequence direction. after sequencing and assembling, we obtained three de novo assembled transcriptome database with three different assembly methods. the number of assembled transcripts, total bases of transcripts, average length of all unigenes, n number of long-transcripts, and number of reads that could be mapped back to transcripts were compared to evaluate the assembly quality with different methods and select the optimum reference transcriptome for further study. the expression level of unigene was measured by the number of clean reads mapped to each annotated unigene normalized to rpkm and adjusted using an additional normalizing factor. functional annotation the protein functional annotation of unigenes was given according to the known proteins with the highest sequence similarity by blastx to public protein databases mentioned above. domain-based alignments were carried out against the cog database to predict and classify possible functions of unigenes. we further used the blast2go program to get go annotation of unigenes. validation of the reference assembly to validate the assembled reference transcripts, unigenes were randomly chosen for further rt-pcr and sanger sequencing. gene specific primers were designed according to the assembled transcripts with primer premier software. the reaction was performed in 25 μl, containing  μl cdna template,  μl ex taq mix, 1 μl of each 20 μm forward and reverse gene-specific primers and 10 μl of pcr-grade water. touch-down pcr reactions were performed as follows: 95°c for 5 min before cycling rounds of 94°c for 30 s, 65°c for 30 s and 72°c for 1 min 30 s – followed by 15 cycles of 94°c for 30 s, 54°c for 30 s and 72°c for 1 min 30 s and finally a 72°c extension step for 10 min. the specificity of products was inspected by 1% agarose gel electrophoresis before they were purified with ezgene™ gel/pcr extraction kit. the purified pcr products were sanger-sequenced by the beijing genomics institute. identification of differentially expressed genes using illumina short reads to perform quantitative comparisons of the transcript variation with respect to different beneficial fungi, total rna samples generated from roots treated with ml, zh3a- or ml01 + zh3a- and mock-inoculated controls at dpi were used to construct un-normalised cdna libraries which were sequenced by illumina hiseq. the cdna library construction, sequencing, transcript profiling and mapping back to the newly assembled reference transcriptome with trinity were carried out as a custom service. the expression level of each assembled transcript in different samples was measured as rpkm values. seqmap was used for read mapping and rseq was applied for rpkm based expression measurement. abundance data from different samples were collected for each transcript. in our work, the changes in relative abundance for unigenes between two samples were screened with the threshold: q-value ≤  and log2ratio ≥  or ≤ - quantitative rt-pcr validation qrt-pcr was further carried out to analyze the expression of twenty- six degs identified by rna-seq technology. gene specific primers were designed according to the cdnas with primer premier software. for rna extraction of each treatment, three biological replicates were collected independently and immediately frozen in liquid nitrogen. one microgram of total rna was reverse-transcribed in a 20 μl reaction mixture from a primescript ii 1st strand cdna synthesis kit. qrt-pcr was performed in a 25 μl reaction mixture containing 2 × sybr master premix ex taq ii  μl, 1 μl of cdna template, and 1 μl of each corresponding primer for the gene of interest. qrt-pcr of three biological replicates for each sample was performed using a lightcycler ii system with its relative quantification software based on the delta-delta-ct method. qrt-pcr was performed for 5 s at 95°c, 10 s at 56°c, and 20 s at 72°c. the cdna samples were standardized to two reference genes: b and b.  co-cultivation of c. hybridumplantlets with different fungi and their respective symbiotic effects on host vegetative growth when co-cultivated with c. hybridum plantlets, the hyphae of the two isolates ml and zh3a- began to grow at dpi and first contacted the hosts’ roots at dpi. after dpi, the co-cultivation media of plantlets inoculated with ml changed from transparent to opaque, whereas the actively growing white fluffy hyphae of isolate zh3a- fully covered the plate. the hyphal growth of isolate zh3a- was retarded by the co-inoculated isolate ml which showed space at the colony margin. after dpi, ml zh3a- and mz treatments all significantly improved the total fresh weight and dry biomass of c. hybridum plantlets compared to the control.table effect of inoculation with epulorhiza repens isolate ml01and/or umbelopsis nana isolate zh3a- on the growth parameters of cymbidium hybridum plantlets at dpi all the data were mean of three biological replicates, ± standard deviation.: percentage of mycorrhizal dependency; prg: the elongated root number/total root number × 100). values followed by different lower-case letters within a column are significantly different at p < according to dmrt. detection the colonization feature of different beneficial fungi on roots of c. hybridumplantlets in an initial observation, the hyphae of two isolates began to contact with orchids roots at dpi. after that time point, we stained the host's roots with each treatment every two days to monitor the colonization process of different fungi and the root staining results of each treatment at and dpi were presented in additional file 4: figure s at dpi, the first contacted roots were stained to confirm the hyphal penetration by isolates ml or zh3a- after dpi, orchid roots were extensively colonized by these fungi. at dpi, the symbiotic roots of orchid plantlets showed uneven staining owing to the rapid growth of roots. to confirm the root staining results and characterize the exact colonization pattern of these beneficial fungi, we further analyzed the microscopic structure of the formed symbionts at dpi. root cross sections of c. hybridum plantlets generally included root hairs, velamen, exodermis, cortex, endodermis, pericycle, vascular cylinder and pith. at dpi, typical pelotons were formed by isolate ml in the host cortex cells whereas hyphae of isolate zh3a- only colonized the velamen cells at regular intervals. roots co-inoculated with isolates ml and zh3a- contained newly formed and degenerated pelotons, and also typical aggregated hyphae of isolate zh3a- demonstrating the co-infection of two isolates. none of these fungal structures were observed in control root tissue.figure microscopic structure of the symbionts formed by different c. hybridum - beneficial fungi interaction at dpi. mock-inoculated control and inoculated with isolate zh3a-, ml01 + zh3a- and ml. ah, aggregated hyphae of isolate zh3a-3; co, cortex cells; dh, degenerated hyphae; ex, exodermis; mh, fungal hyphae of isolate ml01; pe, pelotons; rh, root hair; va, vascular cylinder; ve, velamen cells. bar, 50 μm root transcriptome profile of c. hybridiumgenerated by mrna-sequencing up to date, none of genome sequences for orchidaceae was available in public database. rna-seq, a cost-effective and highly efficient next-generation sequencing technology, provides high throughput and is therefore commonly used for de novo transcriptome assemblies in non-model species, including orchid species. to analyze the transcript variation in c. hybridium colonized by beneficial fungi, pooled total rna from eight samples including mz-treated and mock-inoculated plantlets at each time point were used to generate an informative reference transcriptome database. totally,, base pairs raw data were generated using the illumina hiseq yielding, clean reads that were 90 bp in length. the raw reads data are available at the ncbi sra database with the accession number sra three publicly available assemblers as soapdenovo, abyss and trinity were used to de novo assemble short-read rna-seq data into transcripts. the assembly software tools as soapdenovo and abyss were chosen because of their previous wide application for the similar de novo assembly of rna-seq data sets.trinity method, one of the newest de novo assembly methods, was reported to reconstruct of the majority of full-length transcripts in a sample from rna-seq reads directly, across a broad range of transcript levels. recent publications also showed that trinity had a consistently better performance than the other assemblers as soapdenovo, abyss, trans-abyss, oases with single k-mer value, when measured by transcript accuracy, integrity and completeness, and sensitivity to assemble transcripts from low to high expression levels. the outcomes of three assemblers in our study were summarized in table  as shown in table  de novo assembly using soapdenovo yielded a total of unique transcripts, the largest number of transcripts. in contrast to the soapdenovo and abyss assembly, trinity assembled the longest mean transcript length, the largest n, the largest number of long-transcripts, and the highest percentage of rmbt. we therefore elected to use the trinity assembly to retain maximum information for further transcript analysis and the trinity assembly is also available in additional file table summary of de novo sequence assembly annotation to validate and annotate the assembled unigenes, the unigenes generated by trinity were subjected to blastx searches against public protein databases. we found,, sequences with matches in nr, swiss-prot, and kegg databases, respectively. in total, sequences were annotated and another unigenes showed no homology to known sequences in these databases. among the annotated sequences, approximately transcripts shared the highest similarity to fungal or bacterial sequences in nr or swiss-prot database. the unigenes homologous to known sequences in public databases were further annotated with go terms using blast2go. a total of unigenes were assigned go term annotations in three categories: biological process, molecular function, and cellular component. metabolic process, cellular process and response to stimulus rank among the first three categories among the biological processes. the molecular functions most represented were binding, catalytic activity and transporter activity. cell, cell part and organelle are the most represented categories in cellular components.figure go classification of c. hybridum unigenes. the global transcript abundance could be deduced from a non-normalized cdna library. the transcripts highly abundant in the cymbidium transcriptome are listed in additional file 6: table s there were unigenes mapped into kegg pathways. of the unigenes, were related to metabolic pathways, to biosynthesis of secondary metabolites, to spliceosome, to plant-pathogen interaction, to protein processing in endoplasmic reticulum, and to endocytosis. from the c. hybridum transcript database, of the assembled sequences were assigned to cog annotations, which could be grouped into categories. the largest category was “general function prediction only”, followed by “transcription”, “replication, recombination and repair”, “signal transduction mechanisms”, “posttranslational modification, protein turnover, chaperones” and “carbohydrate transport and metabolism”.figure cog function classification of c. hybridum unigenes. validation of assembled transcripts to experimentally confirm that the unigenes obtained from sequencing and computational analysis were indeed expressed, unigenes were chosen for further rt-pcr and sanger-sequencing. blastn results showed that identities of the selected reference assemblies with sanger-sequencing results were between -100%, confirming the reliability of assembled reference transcriptome. these unigenes were also included in the next qrt-pcr analysis. global analysis of differential gene expression in the symbiotic interaction having generated a c. hybridum root reference transcriptome, our next goal was to perform comparative gene expression profiles in orchid roots inoculated with different beneficial fungi. a total of four independent un-normalised cdna libraries were sequenced using the illumina hiseq 100 bp single read module. for each library, the total number of counts for each read was determined and all the reads were mapped back to the reference transcriptome. in this study, million and million reads were uniquely mapped, representing between % - % of total reads, thus, providing good coverage for the transcript profiles. the expression level of each assembled transcript sequence in different samples was measured through rpkm values, and the significant changes in the transcript abundance were identified applying a q-value ≤  and a two-fold change threshold. based on these criteria, at dpi, or genes were significantly regulated in c. hybridum root inoculated with isolates ml zh3a- or mz, respectively. the general differential expressed profiles in the three interactions are presented in figure  and additional file 10: tables s7-s the top abundant genes in symbiotic roots with positive fold change and top abundant genes in non-symbiotic roots with negative fold change in response to symbiotic fungi were presented in additional file 11: tables s and s om formation specifically induced 57% of the top abundant genes in symbiotic roots with up-regulation pattern. unigene2357_all, encoding metallothionein-like protein type which limited oxidative damage through binding to and detoxifying excess copper and other metals, was the most abundant gene in om symbiotic roots with more than -fold up-regulation. other om-induced genes in this category were largely involved in inhibition of reactive oxygen species, protein metabolism and processing, signaling, cell wall modification and cell defense and recue. twenty-three genes in this category were co-induced by three treatments, involving in signaling, defense and general metabolism. however, 63% of the top abundant genes in non-symbiotic roots with negative regulation pattern in response to inoculants were shared by three treatments and another 29% of these genes were specifically down-regulated by om fungi. the most representative group abundant in non-symbiotic roots and down- regulated by three treatments was attributed to genes encoding main components in photosynthesis system, whereas a mannose-specific lectin, with rpkm in control roots, was specifically down-regulated by om fungi.figure venn diagram of differentially expressed genes in the c. hybridum - beneficial fungi interaction. mock-inoculated, ml01-treated, zh3a-3-treated or ml01 + zh3a- -treated. as shown in the venn diagram, degs were shared by ml01-, zh3a-3- and mz- treatments and transcripts were collectively regulated in cyex and cyex but not in cyex which might specifically functioned in om mycorrhizal formation. individual inoculation with non-mycorrhizal fungus zh3a- specifically regulated genes. another degs with special expression pattern in cyex probably functioned in the combinative symbiotic process. all the degs were grouped according to hierarchical cluster analysis and twenty different expression sub-clusters based on k-means were identified. the lists of all the degs in sub-clusters were also presented in additional file the genes with the same regulation pattern in three treatments were assigned to subcluster and the genes in subcluster were specifically regulated by om formation. genes drastically up-regulated by om fungi were included in subcluster the functions of these genes were related to signaling, reactive oxygen species homeostasis and stress-related, protein metabolism and processing and transport.figure twenty regulatory patterns of all the degs expressed in different c. hybridum - beneficial fungi interactions based on k-means method. go annotation was also used to compare the functional distribution of degs in different treatments. among the represented biological processes, the most represented items were “ion transport” for om formation and “cellular biosynthetic process” for individual zh3a- treatment. among the molecular functions, most represented degs were categorized as “oxidoreductase activity” for om formation and “cellulose synthase activity” for individual zh3a- treatment. among the cellular component, the most represented degs were assigned to “extracellular region” for om formation and “plastid” for individual zh3a- treatment. to identify the biological pathways regulated during orchid- beneficial fungi interaction, we assigned the degs to the reference kegg pathways. the top pathways with the most representation by unique sequences for each treatment were included in figure s7-s. a large number of unigenes related to primary and secondary metabolism were collectively regulated by three treatments as'nitrogen metabolism','flavonoid biosynthesis','cyanoamino acid metabolism','pentose phosphate pathway','starch and sucrose metabolism','alanine, aspartate and glutamate metabolism','stilbenoid, diarylheptanoid and gingerol biosynthesis','porphyrin and chlorophyll metabolism','glyoxylate and dicarboxylate metabolism','metabolic pathways' and'biosynthesis of secondary metabolites'. substantial numbers of genes/proteins involved in signaling pathway, protein turnover, nutrient transport, hormone action, cell rescue and defense were identified in om symbiosis specifically. validation of deg profiling with qrt-pcr using additional biological samples of ml01-, zh3a-3- and mz- colonized roots in comparison to mock-inoculated control roots, we performed real-time rt-pcr to verify the digital expression profiles of genes corresponding to a range of functional categories and regulation patterns in the c. hybridum- beneficial fungi interactions. the pearson correlation coefficient was calculated by spss to assess the correlation between different platforms. overall, the qrt-pcr values were highly correlated with the rna-seq results, which confirmed the general expression pattern for all the selected genes in three treatments.table comparison of the expression profiles of genes by rna-seq and quantitative real-time rt-pcr fold changes of selected c. hybridum trasncripts were logarithmic transformed which were determined by rna-seq and qrt-pcr in root tissues treated with isolates ml, zh3a- and ml01+ zh3a- at dpi as compared to mock-inoculated control. base on the qrt-pcr results, we further analyzed the functional distribution of all the confirmed genes. the function of co-regulated genes in three treatments included signaling, membrane proliferation, scavenging of reactive oxygen species, general stress- responsive, nutrient transport, protein processing, transcription factor and hormone transport. degs derived exclusively from om roots were assigned to signaling, transport, transcription factor, cell defense and hormone metabolism.two genes showed special expression pattern in response to zh3a- treatment, including unigene 4239_all and unigene 13219_all. discussion some different effects of mycorrhizal and non-mycorrhizal fungi on the growth of c. hybridumplantlets this study demonstrated that the vegetative growth of c. hybridum plantlets could be improved by orchid mycorrhizal fungus ml and/or non-mycorrhizal isolate zh3a- it is widely accepted that tulasnella-like rhizoctonia formed typical mycorrhizae with orchids and improved their germination and their seedlings’ growth. consistent with these studies, ml formed typical pelotons in host root cortex and significantly improved the total fresh weight and dry biomass of c. hybridum plantlets compared to the control. however, we observed that the roots of ml01- treated plantlets did not elongate obviously as the co-cultivation proceeded, which might be due to rapid acidification of the co-cultivation media. it was thus inferred ml increased the total dry biomass of orchid plantlets through improving the absorption of mineral nutrient. non-mycorrhizal fungus zh3a- alone showed similar promoting effects on total fw and dw of cymbidium plantlets to that of mz treatment at dpi, but we also observed that individual zh3a-3- treated orchids initiated less new roots than that of dual- inoculated plantlets, and their leaves turned into light green after a prolonged co-cultivation. isolate zh3a- always co-existed with tulasnella-like rhizoctonia in wild orchid roots and the hyphae of isolate zh3a- mainly spread in velamen cells and could not form pelotons in the cortex cells of cymbidium roots. it was also revealed that the effects of isolate zh3a- and another omf isolate on hosts’ uptake of mineral nutrient were complementary in our recent report. thus, we designated zh3a- as mycorrhizal-associated fungi. functional analyses of differentially regulated genes by om symbiosis the establishment of om symbiosis involves the mutual recognition of host and fungi, hypha penetrating the cell walls of root hairs or epidermis and then entering the cells of the cortical parenchma, forming characteristic and complex hyphal coils. as in other endomycorrhizal systems, the plant cell membranes in the infected tissues remain intact at all times, excluding the hyphae from the cytoplasmic compartment, and an interfacial matrix is formed between the invading hyphae and the invaginating plant plasmalemma. the fungi further spread either as a result of repeated colonization or by hyphae penetrating from cell to cell in the root cortex, and the cortex cells may be uncolonized or contain active pelotons or clumps of degenerating hyphae in different proportions. so, rasmussen et al. proposed that there co-existed necrotrophic and biotrophic phase in orchid mycorrhiza, and the former referred to peloton degradation and unilateral nutrient transfer by endocytosis, and the latter, mutual nutritional exchange. a recent labeling study in spiranthes sinensis clearly indicated that in the symbiotic protocorm, hyphal degeneration played a significant role in the mass transfer of the elements to plant cells, and fungal carbon and nitrogen was also detected to be transferred from live hyphae through the interface between the symbionts. owing to the non- synchronous process, the symbiotic root system often contains fungal structures in different developmental stages, and the functions of degs for om symbiosis were widely distributed. genes related to signalingpathway lysm-domain containing receptor-like kinases were identified as receptors of acylated chitin or chitin produced by plant interacting microbes and confirmed to be involved in signaling during nodule and am development. in this context, it is noteworthy that two genes predicted to encode lysm domain receptor-like kinases showed evident up-regulation in om symbiosis. one of them, unigene 1869_all, encodes a protein with 56% identity to lysm type receptor kinase lys from l. japonicus. its induction was validated by qrt-pcr, which measured 16-fold up-regulation in the mycorrhizal roots. this was consistent with the notion that symbiosis signaling seems to occur throughout the process of root colonization by a symbiont and occurs in parallel to additional signaling mechanisms that require the lysm receptor-like kinases. ca2+ spiking in the nucleus and perinuclear region of root hair cells has been documented as one of the earliest cellular responses after the perception of symbionts by host plants. our digital analysis also revealed many transcripts associated with ca2+ signals, as annexin and calcium-binding ef hand family proteins were co-regulated by three treatments. but we also detected some transcripts in this group were specifically regulated by om fungi, including five genes encoding calcium-binding proteins, two genes encoding homologues of calcium-dependent protein kinase, one cyclic nucleotide-gated ion channel protein, and one gene, encoding a homologue of dmi nodulation receptor kinase from m. truncatula. another interesting finding was that two transcripts encoding putative gras family protein with equal 44% identity to ljnsp from l. japonicas, and one gene, which encoded gras family transcription factor with 56% identity to nsp from m. truncatula, were also specifically regulated by om symbiosis. nsp and nsp are required for rhizobium-specific gene expression whereas only nsp might be required for mycorrhizal signaling and nsp is replaced with required for arbuscular mycorrhization to drive expression of mycorrhiza-specific genes. genes involved in cellular organization and protein processing a substantial number of genes related to cell wall modification were found to be co-regulated in three treatments, which encoded probable xyloglucan endotransglucosylase, beta-d-xylosidase or cellulose synthase. two of these were confirmed by qrt-pcr. unigene 8895_all, encoding a protein with 75% identity to the xyloglucan endotransglucosylase of medicago trancatula, was co-induced more than 4-fold. another gene was strongly up-regulated more than 32-fold upon symbiosis. the translated protein shares 92% identity to cellulose synthase a catalytic subunit of vitis vinifera. cesa proteins are part of the cellulose synthase complex in higher plants located in the plasma membrane. guether et al. also confirmed that two transcripts of lotus, ljcel and ljcesa, accumulated specifically in the arbusculated cells and functioned in cell expansion during arbuscule development. the strongly up-regulated transcripts of putative cellulose synthase might be involved in plant cell wall reinforcement during fungal colonization. similar to am fungi, orchid mycorrhizal fungi penetrate the plant cell wall but not the plasma membrane. as the colonization proceeds, the host plasma membrane invaginates and the infected cells often contain numerous mitochondria, ribosomes and profusely developed er. a large number of sequences that were specifically induced by om formation were assigned to protein metabolism and turnover, membrane dynamics and cell wall synthesis. this group of om specifically induced genes in cymbidium were thought to be homologues associated with am and ectomycorrhizal symbioses. besides the co-regulated degs among three treatments, om symbiosis also specifically induced a high number of transcripts encoding cell wall-degrading enzymes, including five glucan endo- -beta-glucosidase, one pectinesterase and one xyloglucan endotransglucosylase/hydrolase protein. concomitantly, three extensin- like proteins, and three putative glycosyltransferases involved in cell wall reinforcement were also found in our dataset. overall, this indicates an importance of wall- strength maintenance in om symbiosis. with regard to membrane proliferation, two homologues of a. thaliana syntaxin syp unigene10103_all and unigene7281_all, were exclusively induced 16- and 6- fold by mycorrhizal symbiosis. syntaxin belong to the superfamily of the snare proteins that play key roles in membrane fusion events during vesicle trafficking. in am symbiosis, a gene encoding a protein with a 60% identity to the a. thaliana syntaxin syp also showed seven-fold up-regulation and was suggested to facilitate plant/fungal compatibility in mycorrhizal lotus roots. protein metabolism and processing seems to be a major process in mycorrhizal symbiosis. most drastically up om- regulated genes in the present study were involved in protein synthesis and processing. concerning protein synthesis, the strong induction of different ribosomal proteins in response to om fungi provided evidence. a number of transcripts associated with protein turnover were also detected in our analysis, including three subtilisin-like proteases, three serine carboxypeptidase, and also some peptidase. some transcripts associated with correct protein folding and degradation of mis-folded protein were also identified as specifically up-regulated in om symbiosis. at the same time, a transcript of putative serine-type endopeptidase inhibitor, unigene16319_all, accumulated more than 40-fold in the mycorrhizal roots compared with controls, which might fine-tune protein processing. transport-related genes phosphate transfer to the adult orchids via mycorrhizal fungi was well documented, just as in other mycorrhizal systems. recently, labeling experiments also directly demonstrated that the fully autotrophic orchid goodyera repens acquires carbon, nitrogen and phosphorous from its fungal partner. notably, g. repens also transfers significant amounts of photosynthate back to its ceratobasidium mycobiont – the first direct demonstration of a net carbon flow from orchid to fungus. tulasnella-like rhizoctonia were also confirmed to improve the plant growth of c. spp. by enhancing nutrient acquisition. in our global profiling analysis, genes associated with membrane transport were identified as co-regulated in three treatments. among them, genes involved in nutrient transport were co-induced in symbiotic roots, including two inorganic phosphate transporters, a plasma membrane atpase, a bcs protein precursor-like a sorbitol transporter, an oligopeptide transporter 1-like and a magnesium transporter, whereas genes were co-repressed at this time point which encoded two putative copper-transporting atpase, a copper transporter, a probable nitrite transporter, a sulfate transporter, an inorganic phosphate transporter 2– chloroplastic, a probable sodium-coupled neutral amino acid transporter and a member of major facilitator family. of these, the expression pattern of five genes was further verified by qrt-pcr. the first of these, unigene6086_all, was co-induced more than 5-fold, which encodes an inorganic phosphate transporter with 81% similarity to inorganic phosphate transporter 1– of arabidopsis thaliana and a shared 78% identity to inorganic phosphate transporter from lotus japonicus and m. truncatula. both deduced amino acid sequences of two up-regulated inorganic phosphate transporter genes in our digital profiling have similarity to reported mycorrhiza-specific or -induced phosphate transporters. the second gene, unigene2386_all, which was up-regulated between 4- and 16- fold in the symbiotic process, encodes a protein with 63% identity to v. vinifera oligopeptide transporter 1-like. another two co-induced genes were involved in carbohydrate transmembrane transport and further identified by qrt-pcr. unigene15955_all encodes a putative sorbitol transporter with 71% similarity to sorbitol transporter of zea mays and unigene4130_all encodes a transporter with 60% identity to sugar transporter erd6-like of a. thaliana. finally, a member of a major facilitator family was confirmed to be repressed more than 4-fold in symbiotic roots and its translated sequence shared 68% similarity to a major facilitator family protein of a. thaliana. in general, the regulation of a number of transporter genes supported the notion that extensive nutrient exchange was involved in symbiosis of the host and beneficial fungi.figure phylogenetic tree for the amino acid sequences of plant phosphate transporters described as mycorrhiza-specific or -induced transporters, and two putative pits of c. hybridum in this study were highlighted in red color. the dendrogram was generated by mega software using clustlw for the alignment and the neighbor-joining method for the construction of phylogeny. the plant pi-transporters accession numbers could be obtained in. we also identified a substantial proportion of genes involved in transport which were specifically regulated in mycorrhizal roots. these included three transcripts encoding putative plant plasma membrane proton atpases. the amino acid sequence of the encoded protein by unigene11510_all shares 77% identity with the am-specific mtha1in m. truncatula. the regulation of another transcript was validated by qrt-pcr and its translated sequence shared 70% similarity to h+-atpase from a. thaliana. in om, like in am and ecm, nutrients from symbiotic hyphae separated at the apoplastic space have to go through several membrane barriers before being assimilated by the partner's cells, so the activated plant plasma membrane proton atpase might be a sign of accelerated membrane transport. the other two genes were also identified as om-specific which encoded a putative copper-exporting atpase and a plasma membrane calcium-transporting atpase, respectively. another group of mycorrhizal-induced genes identified was associated with nitrogen transport. eight genes of peptide transporter were found to be exclusively up-regulated between two-fold and more than 300-fold compared with controls. the strongest up-regulated ptr, encoding a homologue of peptide transporter ptr of m. truncatula with 62% identity, showed none expression level in the control and cyex samples. the second putative ptr was 300-fold up regulated in mycorrhizal roots at dpi and had a 76% similarity to a member of oligopeptide transporter family from p. trichocarpa. quantitative rt-pcr analysis confirmed that one transcript was up-regulated more than 7-fold in om symbiosis whereas down regulated two-fold in zh3a-3-treated orchid roots. besides these ptrs, two putative amino acid transporters were also detected to be induced by mycorrhizal formation. a homology search revealed 70% identity between the translated amino acid sequence of unigene3142_all and cationic amino acid transporter of a. thaliana. the second gene, unigene404_all, encodes a protein with 81% identity to lht of a. thaliana, which was regulated by diverse factors. our data supported the notion that omf played an important role in improving n acquisition of host plants. two putative monosaccharide transporter genes were also found to be induced in om mycorrhizal roots. the closest homologue of translated sequence by unigene12841_all in a. thaliana is sugar-proton symporter plt, with 71% identity at amino acid level. another gene, unigene4472_all, encoded a homologue with 70% similarity to a tonoplast monosaccharide transporter in a. thaliana and tmt1/ probably represents a family of proton-coupled antiporters capable of high-capacity loading of glucose and sucrose into the vacuole. gaude et al. also demonstrated that two putative sucrose transporters and one hexose transporter were induced in non-colonized cortical cells of mycorrhizal roots compared to the cortical cells of non-mycorrhizal roots, and suggested they might be involved in the uptake of carbohydrate from the apoplast into cells in the vicinity of hyphae in arbuscular mycorrhizal roots or functioned in the proton-coupled uptake of carbohydrate into the cytoplasm from the vacuole. this raises a compelling question of what is the exact function of these two putative monosaccharide transporters in our study. finally, a gene encoding a putative potassium transporter was also detected to be 6-fold up-regulated in mycorrhizal roots and shared 50% identity at the amino acid level with the potassium transporter of a. thaliana. in addition to a number of up-regulated genes probably involved in nutrient transport, atp binding cassette transporter genes, which are unlikely associated with direct nutrient transport, were also identified to be regulated in mycorrhizal roots. the strongest up-regulated gene of abc transporter family encoded a protein with 66% identity to pdr abc transporter of oryza sativa, which was regulated by many abiotic and biotic factors. recently, two plant half-size atp binding cassette transporter proteins were identified that are essential for am symbiosis and one full-size abc transporter gene mtabcb showed strongly increased transcript levels in arbuscule-containing cells and also in adjacent cells. these abc transporter proteins may be involved in the transport of signaling between plant cells during am symbiosis. here, these identified abc transporter genes showing different expression patterns in orchid mycorrhizal roots may be important for om symbiosis. om-activated transcripts also included components associated with vesicle-mediated transport, which might involve in the membrane biogenesis for an intracellular accommodation of fungal structures. the translated sequence of unigene7513_all shares 70% identity with vamp7b in a. lyrata subsp. lyrata. another gene, unigene11358_all, encoding a putative clathrin assembly protein, was up-regulated more than 8-fold in mycorrhizal roots. defense -related and phytohormone regulation oxidative burst belongs to early events of plant responses to compatible and incompatible fungi. the success of beneficial microbes to colonize plant roots depends on their ability to manipulate immune responses of the host. this was verified by the findings of transcripts associated with scavenging of reactive oxygen species, such as glutathione peroxidase, probable glutathione s-transferase and germin-like protein, accumulating in the symbiotic roots. one of these, unigene13078_all, encodes a protein with 77% identity to the glutathione peroxidase of populus trichocarpa and was further validated by qrt-pcr. concomitantly, many transcripts encoding general'stress- responsive' proteins; bibenzyl synthase were identified to be induced in these symbiotic process. the enzymes in the pathway towards fungicidal phytoalexin synthesis are produced as a general feature in orchids. reinecke have reported that the activity of benzyl synthase, catalyzing the biosynthesis of bibenzyls and 10-dihydrophenanthrenes of bletilla striata rhizomes, increased substantially upon fungal infection. another interesting finding was that two transcripts encoding programmed cell death protein were down regulated collectively in three treatments. pdcd a suppressor of gene transcription and translation, plays a crucial inhibitory role in several types of human tumors through inhibiting autophagy in multiple cell types both in vitro and in vivo. down- regulation of these homologues of pdcd may promote the symbiotic process through suppressing plant defense. furthermore, four transcripts, encoding homologues to key molecular regulators in two-component systems which play important roles in responses to environmental stress stimuli, were found to be down regulated regardless of the inoculated fungi species, highlighting the shared strategy by different beneficial fungi in the symbiosis. apart from these shared degs, three chitinase genes were exclusively induced by om fungi and the specific expression pattern of unigene 15958_all was confirmed by qrt-pcr. the translated sequence of unigene 15958_all shared 59% identity with the chitinase 4-like of fragaria vesca. chitinases hydrolyze β- -glycosidic bonds and can be induced in pathogenic interactions. but the expression of unigene 15958_all was not induced by a pathogenic fusarium isolate. unigene8905_all encoded a protein with 55% identity to class iii chitinase of m. truncatula and class iii chitinases were found to be expressed predominantly in the arbuscule-containing cells, where they may suppress plant defense by the reduction of chitin-like elicitors during the formation of functional symbiotic interfaces. the similar role of class iii chitinase svchit in om was recently suggested by perotto. thus, these up-regulated chitinases are probably involved in the lysis of fungal cell wall and balancing the symbiosis. another notable om-regulated group was mannose-specific lectins. among these, unigene4318_all and unigene529_all showed high-level expression levels in mycorrhizal roots with little expression in the un-inoculated control and zh3a- treatments. antifungal activity against a range of phytopathogenic fungi has been demonstrated for some of mannose-specific lectins both in vitro and in transgenic plants, and they were thus proposed to be part of defense mechanisms against the invading fungus. the strong up-regulation of svlect and svlect in response to mycorrhizal fungi t. calospora in s. vomeracea protocorm was recently confirmed by qrt-pcr and this was suggested to be part of a mechanism to limit fungal growth in the mycorrhizal protocorm. however, another mannose-specific lectin gene, unigene14222_all, with rpkm in control roots, was specifically down-regulated by om fungi. owing to the multi-functions of plant lectins, the role of these lectins in the orchid mycorrhizal process remains to be established. among plant hormones, ethylene, salicylic acid, abscisic acid, and jasmonic acid are known to be key elements in fine-tuning the plant defense response during interaction with other organisms. the inhibitory effect of exogenous ethylene in am development has been documented and the colonization of host tissues by mutualistic am fungi such as rhizophagus irregularis, ectomycorrhizal fungus laccaria bicolor, or endophytic fungus piriformospora indica was inhibited by et. however, et signaling is still essential for nodule formation in sesbania rostrata and in the initial stages of p. indica colonization. because ethylene-responsive transcription factors can function as transcriptional activators and repressors, they are candidates for establishing a balanced defense response to the fungus without preventing growth and development in am symbiosis. valadares also suggested the role of ethylene in om regulation when compared the gene expression pattern of o. maculata mycorrhizal roots in comparison to non-mycorrhizal controls. in our global expression profiling, eight transcription factors related to ethylene signaling were found to accumulate in all the symbiotic roots and four genes among these were specifically induced in om roots. among these, the induction of two transcripts was further validated by by qrt-pcr. the translated sequence of unigene1724_all showed 63% similarity to the ethylene receptor homolog of zea mays with two-fold induction. unigene12197_all, encoding a protein with 45% similarity to ethylene-responsive element binding factor of g. max, was induced more than 5- fold in the symbiotic roots. as with plett et al., the induction of these tfs is probably an adaptive response by the plant so that its growth and vigor are not compromised by the fungus in the late stage of symbiosis.table transcription factors among beneficial fungi- induced genes lists of tfs in this table included specifically regulated tfs in om formation and * co-regulated genes in three treatments; ⊕qpcr validation; number in parenthesis referred to the gene number regulated by om formation. aba plays a crucial role in biotic and abiotic stress responses. analysis of am colonization in the sitiens aba-deficient tomato mutant showed that aba was necessary in order to complete arbuscule formation and its functionality and to promote sustained colonization of the plant root. subsequent work supported the notion that aba deficiency negatively regulates mycorrhizal intensity/hyphal colonization directly or indirectly through increasing ethylene levels. accumulation of the transcripts related to aba de novo biosynthesis was documented in am roots and om green protocorms. at the same time, genes associated with aba catabolism also showed increased transcript levels in am roots and om roots. in our study, transcripts related to aba biosynthesis: 9-cis-epoxycarotenoid dioxygenase and zeaxanthin epoxidase, were co-suppressed in symbiotic roots. aba deficiency might be involved in restriction of fungal spread in orchid symbiotic roots. jasmonic acid has long been implicated in the plant’s systemic response to pathogen attack. while experimental results sometimes appear contradictory, most authors agree that jasmonate signalling is important for am colonization and development. plett revealed that in the interaction between l. bicolor and poplar, ethylene and jasmonic acid showed extensive transcriptional cross-talk and functioned as negative modulators during mutualistic symbiosis. jasmonates belong to a diverse class of lipid metabolites known as oxylipins and the key genes in the biosynthesis of oxylipins were commonly induced in mycorrhizal roots. consistent with these studies, the induction of lipoxygenase genes in our study suggest the involvement of the oxylipin pathway in regulation of om symbiosis. in plants, there are two main branches of the oxylipin pathway, determined by two different types of lipoxygenases, 9-lox and 13-loxs. the 13-lox pathway leads to the biosynthesis of ja and derivatives, and the relevance of the 9-lox pathway in plant interactions with nematodes and pathogens has also been demonstrated recently. further research is required to elucidate the exact role of these up-regulated loxs in long-term balance of the om symbiosis. apart from the potential roles for et, aba and ja in om symbiosis, three transcripts associated with gibberellin biosynthesis were also increased in abundance in orchid mycorrhizal roots: two genes encoding gibberellin 3-oxidase and another encoding ga 20-oxidase. the translated sequence of unigene13170_all shares 77% identity to gibberellin 20-oxidase of triticum aestivum and its induction was confirmed by qrt-pcr. the product of ga3ox and ga20ox activates ga, suggesting that ga synthesis is essential for om formation. in barley- p. indica interaction, schäfer also observed that two barley mutants impaired in ga synthesis or perception, reduced colonization by p. indica. similarly, homologues of genes involved in the ga biosynthesis were already identified as mycorrhizal-responsive in tomato plants at the mature phase of the interaction. foo reported that biologically active gibberellins suppress arbuscule formation in pea roots, and gibberellin signaling della proteins are essential for this response. owing to extensive transcriptional cross-talk between ga signaling and sa/ja responses, the ga-regulated defense response might be critical for establishing these long-term mutualisms. transcription factors as detailed above, three treatments initiated common stress responses for symbiotic roots, so the majority of co-regulated tfs encoded ap domain tfs and ethylene-responsive tfs, wrky tfs, myb- family tfs, basic helix-loop-helix tfs and zinc-finger tfs. dna-binding wrky, ap2-domain and erf tfs are generally involved during pathogen-related stress responses. the co-regulated tfs may also be involved in root developmental process, as are myb family, erf and nac family. we also detected gras family and arf tfs specifically induced by om formation. four transcripts encoding for gras members were putative scarecrow protein encoding genes, considered to be genetic determinants of root identity. previous studies also confirmed that members of the gras family are essential for nodule development and they may be important for the regulation of gene expression in am roots of m. truncatula, l. japonicus and s. lycopersicum. two putative arf transcripts were also identified to be up- regulated in orchid mycorrhizal roots. the induction of putative arfs has also been found during am symbiosis in maize, rice and m. truncatula, but not in l. japonicus. and in am interaction, mycorrhizal roots exhibit morphological characteristics such as an increase in the number of lateral/fine roots during early growth phases, similar to auxin-treated roots. recently, sukumar and fusconi summarized a compelling supportive role for involvement of auxin pathways in beneficial plant–microbe interactions. chutima reported that orchid-associated fungi tulasnella sp. produced high levels of indole-3- acetic acid with all biological activities in culture medium supplemented with 2 mg/ml of l-tryptophan. the iaa production by isolates ml was also identified in our research group.furthermore, unigene11039_all, predicted to translate a protein with 66% similarity to an auxin efflux facilitator slpin of solanum lycopersicum, showed 15-fold up-regulation in the symbiotic roots. all these data point to the probable involvement of auxin signaling pathway in om functionality. degs specifically regulated by non-mycorrhizal fungus zh3a-3 the hyphae of non-mycorrhizal fungus zh3a- colonized the velamen cells of orchid roots and individual inoculation significantly improved the root length and fw of c. hybridum plantlets. as detailed above, genes were significantly regulated in c. hybridum roots inoculated with zh3a- alone and most of these degs were shared with those of ml and/or mz treatment. apart from these shared degs, we also identified genes were specifically regulated by zh3a- treatment. the largest number of up-regulated transcripts by zh3a- treatment was associated with orchid root development. these genes are mainly related to hormonal balance and transcription regulation. among plant phytohormones, auxin, cytokinin, and ethylene play important role for root development. consistent with this, we detected the up-regulation of many corresponding genes, which encoded one cytokinin dehydrogenase, two erf, and three abc transporter b family members. auxin plays a pivotal role in many plant- beneficial microbe interactions, and it is well established that higher auxin/cytokinin ratio determines the root initiation. a recent study also clearly demonstrated that osckx a cytokinin oxidase/dehydrogenase family gene, integrates cytokinin and auxin signaling to control rice crown root formation positively. so, the up-regulation of unigene14465_all, which encodes a homolog with 69% identity to ckx of a. thaliana, may suggest its positive role in orchid- zh3a- interaction. polar auxin transport is essential for establishment of auxin gradient that spatially defines the quiescent center and consequently the fate of neighboring stem. abc transporter b family members were probably involved in auxin acropetal or basipetal transport, thus, the up-regulation of three abc transporter b family members appears to converge on the auxin pathway to modulate root development. as shown in table  the putative myb-tfs ranked the third largest number of the regulated tfs in the mycorrhizal symbiosis, and most were shared by three treatments, which suggested their essential role in response to symbiotic fungi. recently, a convincing research supported that a putative myb-like tf of l. japonicus which was am-responsive may also have non-symbiotic functions: i.e. root growth. in our study, two transcripts,encoding homologues of myb family tfs of o. sativa japonica group and a. thaliana with 60%- 80% identity, were also detected to be up-regulated 4-fold by zh3a- the translated sequence of unigene2947_all also shared 46% similarity to ljmami at amino acid level. another up-regulated tf related to root development was unigene489_all, which encodes a protein with 45% similarity to scarecrow-like protein of z. mays. we also identified a substantial number of genes related to general disease resistance processes which were specifically up-regulated by zh3a- treatment. these included transcripts related to signaling, transcripts associated with balance of reactive oxygen, and transcripts involved in secondary metabolism. the increased activity of chitinases were generally found in nonspecific, broad-spectrum defenses, but we did not detect the transcript accumulation of chitinase genes in cyex the fungal cell wall of isolate zh3a- always ruptured as the cultivation proceeded, which might be the reason. the activation of basal defense mechanisms of orchid plants by beneficial fungus zh3a- did not affect the symbiotic effects, as described in. as zh3a- colonization induced dramatic structural changes in root cells and promoted the root growth of orchids, transcripts accumulated in the symbiotic roots were also related to carbohydrate and amino acid metabolism, protein synthesis and processing, and structural re-organization. furthermore, a transcript involved in water transport was also up-regulated, which encodes a homolog of aquaporin pip1- of z. mays with 87% identity. taken together, we confirmed the findings which suggested that the development and metabolism of plant symbionts are generally driven by differential regulation of transcriptional regulators, signal transduction, and metabolic pathways, rather than by expression of symbiosis-specific genes. detection of putative fungal genes in symbiotic roots owing to the distinct phylogenetic relationship between isolate ml and zh3a- we could detect possible fungal genes combining the annotation information of de novo assembled transcripts with the comparative expressed level in different symbiotic roots. of these genes, were probably genes from isolate zh3a- and the rest, from ml the origin of putative ml genes was further confirmed by blastn search against the t. calospora draft genome or t. calospora mycelium ests with an e value > 1e-. transcripts from ml in the mycorrhizal roots were largely assigned to ribosomal proteins, enzymes that might be involved in plant cell wall degradation and remodeling the fungal cell wall during growth and symbiosis. we also identified numerous transcripts encoding different proteins associated with the nutrient transport and assimilation, such as h+-transporting atpase, phosphate transporter, ammonium transporter, amino acid and protein transporters, zip-like iron-zinc transporter, monosaccharide transporters and mfs-type transporters. abc-type fe3+-siderophore transporter and other transporters most probably implicated in the detoxification process were also detected in the mycorrhizal roots. the translated amino acid sequence of unigene5906_all shared 59-60% identity to high-affinity pi- transporters from g. versiforme and r. intraradices. the function of rinpt has been suggested by tisserant et al. to be involved in stabilizing cooperation in the mycorrhizal symbiosis. the transcript abundance of unigene5906_all in cyex was higher than that in cyex so, this gene might be responsible for the better growth- promoting effect of mz treatment than that of ml alone. besides the transcripts in n and p assimilation pathway, iron-zinc transporter and abc-type fe3+-siderophore transporter co-existed in the mycorrhizal roots, which suggest an important role of some metals in om functionality. our previous study showed that isolate cl tulasnella-like rhizoctonia, formed typical mycorrhizae with orchids and enhanced the n, p, zn, fe contents of mycorrhzial plants. the its sequence of ml shared 84% similarity to that of isolate cl and both functioned similarly in improving the growth and nutritional absorption of orchid plantlets. however, there is a need to compare the differences in gene expression between the intraradical and extraradical mycelium so as to elucidate the role of these transcripts in om symbiosis.  using a combination of rna-seq and real-time rt-pcr, we characterized the root transcriptomic reprogramming of c. hybridum with respect to different beneficial fungi and identified specific changes in relative abundance for transcripts of genes in om symbiosis. genes involved in cell wall modification, reactive oxygen species detoxification, defense-related phytohormone and phosphate transport were co- induced in all the symbiotic interactions. the transcripts increased specifically by om were related to signaling, cell wall degradation or reinforcement, protein metabolism and processing, defense and n, fe transport. among these degs, genes potentially functioned in c, n, p and fe transport, hormone metabolism and signaling are of particular interest. aside from these orchid transcripts, we also identified some putative fungal genes in symbiotic roots which were associated with plant cell wall degradation, remodeling the fungal cell wall and nutrient transport. this comprehensive dataset provides a basis for future research in exploring the molecular mechanisms underlying om formation and functionality. electronic supplementary material additional file 1: table s1: identities of the selected reference assemblies with sanger-sequencing results of corresponding touch-down pcr products. additional file 2: table s2: gene-specific primers used for qrt-pcr. additional file 3: figure s1: typical co-cultivation condition of cymbidium hybridum plantlets with different fungi at dpi. mock-inoculated control and inoculated with isolate ml, zh3a- and ml01 + zh3a-. additional file 4: figure s2: root staining results of each treatment at different symbiotic time points. roots inoculated with e. repens isolate ml at 6 days post-inoculation; roots inoculated with u. nana isolate zh3a- at dpi; roots inoculated with isolates ml and zh3a- at dpi; mock-inoculated roots at dpi; roots inoculated with e. repens isolate ml at dpi; roots inoculated with u. nana isolate zh3a- at dpi; roots inoculated with isolates ml and zh3a- at dpi; mock-inoculated roots at dpi; roots inoculated with e. repens isolate ml at dpi; roots inoculated with u. nana isolate zh3a- at dpi; roots inoculated with isolates ml and zh3a- at dpi; mock-inoculated roots at dpi; roots inoculated with e. repens isolate ml at dpi; roots inoculated with u. nana isolate zh3a- at dpi; roots inoculated with isolates ml and zh3a- at dpi; mock-inoculated roots at dpi. additional file 5: cymbidium hybridum root transcriptome assembly- illmina 90 bp reads. additional file 6: table s3: highly abundant transcripts in cytr-unigenes. additional file 7: table s4: pathway annotation of cytr-unigenes. additional file 8: table s5: statistics of dge sequencing. additional file 9: table s6: statistics of uniquely mapped reads for each sample. additional file 10: table s7-s9: lists of differentially expressed genes in different c. hybridum- beneficial fungi interactions. additional file 11: table s10: top abundant genes in symbiotic roots with positive fold change in response to symbiotic fungi. table s top abundant genes in non-symbiotic roots with negative fold change in response to symbiotic fungi. additional file 12: figure s3: hierarchical clustering of degs from c. hybridum expressed in the roots in response to isolates ml, zh3a-, ml01 + zh3a- or mock inoculation. additional file 13: lists of degs in sub-clusters based on k-means. additional file 14: figure s4-s6: go enrichment bar graphs of degs in the interactions between c. hybridum and ml, zh3a- or ml01 + zh3a-. additional file 15: table s12-s14: lists of enriched pathways in different c. hybridum- beneficial fungi interactions. additional file 16: figure s7-s9: scatterplot for top enriched kegg pathways of degs in the interactions between c. hybridum and ml, zh3a- or ml01 + zh3a-. additional file 17: figure s10: comparison of expression profiles of random selected genes by rna-seq and qrt-pcr. additional file 18: table s15: expression data for putative genes from e. repens isolate ml in symbiotic roots. table s expression data for putative genes from u. nana isolate zh3a- in symbiotic roots. additional file 19: figure s11: phylogenetic tree for the amino acid sequences of symbiosis- associated phosphate transporters from plant and fungi.the dendrogram was generated by mega software using clustlw for the alignment and the neighbor-joining method for the construction of phylogeny. the plant pi-transporters accession numbers could be obtained in. fungal phosphate transporters accession numbers: glomus versiforme: gvept, aac; gvept gi|74654973; rhizophagus intraradices: rinpt, aal; yeast: yepho gi|1346710; sesbania rostrata: serpt, cac; piriformospora indica: pipt abi; pipt cca; pipt cca; pipt cca. the nucleotide sequences of putative phosphate transporters in our study were available in additional file 6: table s abbreviations abcatp-binding cassette amsarbuscular mycorrhizas ap2/erfapetala2/et response factor arfauxin response factor bhlhbasic helix-loop-helix cat5cationic amino acid transporter 5 cdpkcalcium-dependent protein kinase ckxcytokinin oxidase/dehydrogenase degsdifferentially expressed genes dgedigital gene expression ecmsectomycorrhizas etrethylene receptor homolog lht1lysine histidine transporter 1 lysm-rlklysm-domain containing receptor-like kinases mafmycorrhizal- associated fungi mdmycorrhizal dependency mhheterotrophy mxmixotrophs naataanicotianamine aminotransferase a ngsnext generation sequencing nsp2nodulation-signaling pathway protein omorchid mycorrhizae omforchid mycorrhizal fungi optoligopeptide transporter family patpolar auxin transport pdapotato-dextrose agar pdcdprogrammed cell death protein 4 ptphosphate transporter ptrpeptide transporter prgpercentage of root growth qcquiescent center ram1protein required for arbuscular mycorrhization 1 rmbtreads that could be mapped back to transcripts rpkmreads per kilo base per million reads scrscarecrow spasuppressor of phya-105 tfstranscription factors tmt2tonoplast monosaccharide transporter 2 vamp7bvesicle-associated membrane protein 7b xetxyloglucan endotransglucosylase. competing interests the authors declare that they have no competing interests. authors’ contributions zxl and zjx carried out the experiments, performed the bioinformatics analyses, and drafted the manuscript. ccl cultured the experimental material and extracted all the rna samples. zhy performed the symbiotic experiments. yjz and lm participated in the qrt-pcr experiment. lfb designed the study and revised the manuscript. all authors read and approved the final manuscript. we are grateful to xiaoyong liu for his help with the identification of fungal isolates in our research. we thank professor ronald r. sederoff for his critical reading and revision of the manuscript. we also thank for the constructive comments of the anonymous reviewers to improve our manuscript. this research was financially supported by national natural science foundation of p. r. china, guangdong province science and technology project, and the cooperation project in industry, education and research of guangdong province and ministry of education of p. r. china, and the science and technology fund from south china agricultural university.
more than, clean reads were obtained from un-normalized cdna library prepared from fungal- and mock- treated cymbidium roots at four time points using rna-seq technology. a total of sequences were annotated based on blastx results and over ninety percent of which were assigned to plant origin. aside from these orchid transcripts, some putative fungal genes were also identified in symbiotic roots related to plant cell wall degradation, remodeling the fungal cell wall and nutrient transport. genes involved in phosphate transport and root morphogenesis were also detected to be up-regulated collectively. these reads were assembled into unique transcripts, with a mean length of bp. the orchidaceae is one of the largest families in the plant kingdom and orchid mycorrhizae are indispensable in the life cycle of all orchids under natural conditions. our previous work demonstrated that the non-mycorrhizal fungus umbelopsis nana zh3a- could improve the symbiotic effects of orchid mycorrhizal fungus epulorhiza repens ml by co-cultivation with cymbidium hybridum plantlets. among the om specifically induced transcripts, genes related to signaling, protein metabolism and processing, defense, transport and auxin response were identifed. the molecular basis of om formation and function will improve our knowledge of plant- mycorrhzial fungi interactions, and their effects on plant and fungal growth, development and differentiation. our analysis showed that general defense responses were co- induced by three treatments, including cell wall modification, reactive oxygen species detoxification, secondary biosynthesis and hormone balance. the comparative expression profiling revealed that the transcriptional reprogramming by om symbiosis generally overlapped that of arbuscular mycorrhizas and ectomycorrhizas. thus, we investigated the c. hybridum transcript profile associated with different beneficial fungi.
 the orchidaceae is one of the largest families in the plant kingdom and orchid mycorrhizae are indispensable in the life cycle of all orchids under natural conditions. in spite of this, little is known concerning the mechanisms underlying orchid- mycorrhizal fungi interactions. our previous work demonstrated that the non-mycorrhizal fungus umbelopsis nana zh3a- could improve the symbiotic effects of orchid mycorrhizal fungus epulorhiza repens ml by co-cultivation with cymbidium hybridum plantlets. thus, we investigated the c. hybridum transcript profile associated with different beneficial fungi.  more than, clean reads were obtained from un-normalized cdna library prepared from fungal- and mock- treated cymbidium roots at four time points using rna-seq technology. these reads were assembled into unique transcripts, with a mean length of bp. a total of sequences were annotated based on blastx results and over ninety percent of which were assigned to plant origin. the digital gene expression profiles in cymbidium root at days post inoculation revealed that and genes were sigificantly regulated in response to ml zh3a- and ml01+ zh3a- treatments, respectively. twenty-six genes in different regulation patterns were validated using quantitative rt-pcr. our analysis showed that general defense responses were co- induced by three treatments, including cell wall modification, reactive oxygen species detoxification, secondary biosynthesis and hormone balance. genes involved in phosphate transport and root morphogenesis were also detected to be up-regulated collectively. among the om specifically induced transcripts, genes related to signaling, protein metabolism and processing, defense, transport and auxin response were identifed. aside from these orchid transcripts, some putative fungal genes were also identified in symbiotic roots related to plant cell wall degradation, remodeling the fungal cell wall and nutrient transport.  the orchid root transcriptome will facilitate our understanding of orchid - associated biological mechanism. the comparative expression profiling revealed that the transcriptional reprogramming by om symbiosis generally overlapped that of arbuscular mycorrhizas and ectomycorrhizas. the molecular basis of om formation and function will improve our knowledge of plant- mycorrhzial fungi interactions, and their effects on plant and fungal growth, development and differentiation. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords root transcriptomedigital gene expressionplant mycorrhizal symbiosiscymbidium hybridumepulorhiza repens ml01umbelopsis nana zh3a-3issue-copyright-statement© the author 2014 