Growth Promotion-Related miRNAs in Oncidium Orchid Roots Colonized by the Endophytic Fungus Piriformospora indica

Piriformospora indica, an endophytic fungus of Sebacinales, colonizes the roots of a wide range of host plants and establishes various benefits for the plants. In this work, we describe miRNAs which are upregulated in Oncidium orchid roots after colonization by the fungus. Growth promotion and vigorous root development were observed in Oncidium hybrid orchid, while seedlings were colonized by P. indica. We performed a genome-wide expression profiling of small RNAs in Oncidium orchid roots either colonized or not-colonized by P. indica. After sequencing, 24,570,250 and 24744,141 clean reads were obtained from two libraries. 13,736 from 17,036,953 unique sequences showed homology to either 86 miRNA families described in 41 plant species, or to 46 potential novel miRNAs, or to 51 corresponding miRNA precursors. The predicted target genes of these miRNAs are mainly involved in auxin signal perception and transduction, transcription, development and plant defense. The expression analysis of miRNAs and target genes demonstrated the regulatory functions they may participate in. This study revealed that growth stimulation of the Oncidium orchid after colonization by P. indica includes an intricate network of miRNAs and their targets. The symbiotic function of P. indica on Oncidium orchid resembles previous findings on Chinese cabbage. This is the first study on growth regulation and development of Oncidium orchid by miRNAs induced by the symbiotic fungus P. indica.


Introduction
Piriformospora indica, a root-colonizing endophyte with a broad host range, is intensively studied due to its diverse beneficial effects on the performance of both monocot and eudicot plants. P. indica improves nutrition uptake of the host plants [1,2], enhances their resistance to biotic and abiotic stress [3][4][5][6], promotes the biomass of the aerial and under-ground parts [7] as well as the accumulation of secondary metabolites followed by a cell deathdependent phase, and this lifestyle is critical for local and systemic resistance induced by P. indica [6,[8][9][10][11]. Since the fungus can easily be cultivated without any host, it became a model organism to study symbiotic root interactions [3]. Proteomics and transcriptomics were employed to study the interaction of P. indica with Arabidopsis thaliana and barley, and stage-specific up-or downregulated proteins and genes involved in phytohormone metabolism [11,12] and membrane protein synthesis [13] were detected. P. indica also promotes growth of Chinese cabbage (Brassica campestris subsp.chinensis). Analysis of a double-subtracted EST cDNA library identified genes associated with auxin biosynthesis and signaling indicating that this phytohormone might be crucial in the Chinese cabbage/P. indica symbiosis [14].
Oncidium orchid, also known as ''dancing lady orchid'', is one of the most important pot and cut flowers in Taiwan, Southeast Asia and China. Oncidium belongs to the epiphytic orchids inherently known for slow growth rate. They require 2-3 years from the seedling stage to flowering [15]. Besides, the most commercialized cultivar Onc. 'Gower Ramsey' (Onc. GR) is highly susceptible to the soft-rot (Erwinia carotovora) disease [16]. Since most of the commercial Oncidium cultivars are self-incompatible [17,18] or show high sterility which prevents breeding by cross-hybridization, it is necessary to improve new varieties to overcome these problems. Although gene transformation technics to improve the Oncidium cultivar property have been developed [19][20][21], this technology is highly dependent on the available information on the genetics of Oncidium. Recently, some progress has been made with regard to Oncidium flowering time [22][23][24][25][26] and pigment biosynthesis [19,20,27,28], but little is known about growth regulation and disease resistance in orchids. We observed that colonization of Oncidium roots by P. indica results in growth promotion, similar to results obtained for Chinese cabbage and other plant species. P. indica is phylogentically related to fungi isolated from the rhizosphere of orchids, and helpful for the germination of the orchid seeds in vitro [29]. However, the underlying molecular mechanisms establishing and maintaining the symbiosis are barely studied. The identification of genes responding to the colonization may help to elucidate the molecular basis of orchid growth and development.
Recent evidence indicates that micro RNAs (miRNAs) play an important role in the interaction between plants and soil microbes [30]. MiRNAs are a class of endogenous small non-coding RNAs, often 20-22 nt long, which regulate gene expression by mediating gene silencing at transcriptional (TGS) and post-transcriptional (PTGS) levels, including histone modification, DNA methylation, RNA silicing and translational repression [31,32]. They play important roles in numerous plant processes including development, differentiation and response to biotic and abiotic stresses [32][33][34][35][36][37][38][39]. Furthermore, miRNAs are rapidly upregulated in response to wounding and other stress stimuli [40]. The biogenesis and mode of action of miRNAs have been extensively reviewed [39,41].
MiRNAs recognize their mRNA targets based on near-perfect complementarity and suppress expression of the target genes by guiding degradation or translational repression of the cognate mRNA targets [42]. Furthermore, they tend to be largely acting as an ''early'' regulator of signal transduction in operating at the transcription factor (TF) level in various systems [43,44]. During symbiosis, miRNAs are involved in regulating plant nutritional balance [45,46], hormone homeostasis and signaling [47] and symbiotic nodule spatial and temporal development [48,49]. MiRNAs can respond quickly to infection by symbiotic bacteria. In soybean roots, a set of miRNAs were found to be intensively upor down-regulated by infection with the rhizobial bacterium Brodyrnia japonicum where they target a wide range of mRNAs [47,49].
Traditional cloning and studying of miRNAs rely on the cloning technology,which is labor-intense and low efficient. Due to the development of next generation sequence technology, genomewide detection of miRNAs became much easier. Moreover, the read numbers of unigenes among different libraries helps to estimate the expression level of an individual small RNA [50]. Here, we used high-throughput technology to study the miRNAs in roots of the Oncidium orchid colonized or not colonized by P. indica. A great number of conserved and novel miRNAs were identified. The predicted target genes of these miRNAs were mainly cataloged to phytohormone signal perception and transduction, transcription factors, secondary metabolites and plant defence mechanisms. This study was focused on the identification and characterization of conserved and novel miRNAs. The target transcripts of some Oncidium miRNAs were identified. The data revealed that the identified miRNAs may participate in establishing an intricate network which regulates plant growth, root development and defense. This study may help to understand the molecular basis of the symbiotic interaction of P. indica with orchid plants.

Results
The Effect of P. indica on Oncidium GR Growth and Root Development To understand the symbiotic interaction between Oncidium GR and P. indica, root tissue was analyzed after fungal colonization. P. indica penetrated into the root epidermal layers 24 h after inoculation ( Figure 1B-D). Five days after inoculation, the hyphae were widely distributed over the root surface, and no significant difference was observed between root tip, elongation zone and differentiation zone ( Figure 1E). This case is different from A. thaliana [6,51] and barley [10], where hyphae preferentially colonized the differentiation zone. Cross and longitudinal sections show that hyphae fully colonized the velamen ( Figure 1F, H). Since the velamen and exodermis do not contain plastids, the hyphae were easily detected after visualization by the carbohydrate binding lectin concanavalin A-AF633 (conA-633) ( Figure 1H). Moreover, unlike other endophytic fungi which penetrate into the cortex cells through the exodermis [52], no hyphae were detected in the exodermis and cortex layer of Oncidium, even after a relative long period of co-cultivation (8 weeks). This suggests that P. indica may preferentially colonize dead cells in Oncidium, consistent with findings in barley where massive development of P. indica takes place in dead host cells [10].
Eight weeks after inoculation, the stems and roots developed better than orchids without P. indica co-cultivation. Furthermore, no necrotic lesions were observed ( Figure 1A). With the colonization by P. indica, the fresh weight of the orchid seedlings is 2-fold higher than that of uninocculated control seedlings. Also the leaf number was more and the stem diameter was slightly thicker compared to the control. The number of the roots markedly increased approximately 2-fold and the diameter of the main root was 1,2-fold bigger. These results indicated that the establishment of the beneficial symbiosis is relatively slow and mainly caused by an increase in the root biomass.
Sequencing Small RNA from Oncidium Roots after Cocultivation with P. indica Roots colonized or not colonized by P. indica for 8 weeks were chosen for small RNAs sequencing. This time point was analyzed because a growth-promoting effect of the fungus is visible. After removal of the low quality contaminant and adapter reads, 24,570,250 and 24,744,141 clean read sequences were obtained from the libraries of control and colonized roots (Accession: SRP031471). The clean read sequences represented 97.15% and 99.66% of all reads, respectively. Their lengths range from 13 to 30 nucleotides (nt). Small RNAs of 20 to 24 nt represented 95.02% and 75.66%, respectively ( Figure S1), of all small RNAs, indicating that both libraries are of high quality and can be used for further miRNA studies.
For both libraries, 17,036,953 unique sequences were obtained from a total of 49,314,391 clean reads. Among them, 62.56% were specific for the control roots, 27.47% were specific for the library of P. indica-colonized roots, and 9.97% overlapped in both libraries. After alignment to the Rfam 10.1 and Genbank databases, rRNA, snRNA and other no-coding RNA sequences were removed. The remaining 43,190 and 33,619 unique sequences in the two libraries were considered as potential miRNAs. However, most of the remaining unique sequences were still un-annotated. They accounted for 99.09% and 97.40% in the two libraries, respectively ( Table 1). The high percentage of unannotated small RNAs is mostly due to the insufficient genomic information of the Oncidium orchid. Error bars represent SD for three independent experiments. *, P value,0.05; **, P,0.001. Hyphae were stained with chitin-specific WGA-AF488 (green). Penetrated sites (Ps) were stained with lectin-specific conA-AF633 (red) [88]. Samples were analyzed and photographed with an Olympus IX71 inverted microscope system (Japan). Ve, velamen; EX, exodermis; Cort: cortex; Ch: Chloroplast. Hy, hyphae; Ps, penetration site. doi:10.1371/journal.pone.0084920.g001 For identification of conserved miRNAs, those with 20 to 24 nt were aligned with all known mature miRNAs from plants in the miRBase 19.0 Databank (E value = 1). Sequences with more than 4 nt differences were removed. Afterward, 2,789 from 17,036,953 unique sequences from both libraries showed homology to 86 miRNA families distributed in 41 species (Table S2). They are mainly found in Cucumis melo, Malus domestica, Saccharum ssp., and Brachypodium distachyon ( Figure S2). Furthermore, it is well known, that miRNA* strands bound to Ago1 typically generate miRNAs starting with uridine, while Ago2-bound miRNA* usually starts with cytidine [53]. In our two libraries, most conserved miRNAs started with ''U'' ( Figure S3), indicating that the miRNAs in Oncidium roots are mainly regulated by Ago1.
The vast majority and most abundantly conserved miRNAs were present in both libraries, and the number of miRNAs detectable only in one of the two libraries was low (Table S2). 1605 miRNAs belonging to 62 families were detected in both libraries. 1083 miRNAs belonging to 56 families were detected only in the library from P. indica-colonized roots. 101 miRNAs belonging to 31 families were detected only in the library from control roots. The miRNAs specifically detected in the library from P. indicacolonized roots were mir158a (t0035494; 41 reads), mir166l (t0039545; 36 reads) and mir528 (t0043631; 33 reads). The read number of those miRNAs specifically expressed in P. indicacolonized roots suggests that their accumulation is low (Table S2).
Thus, miRNAs specifically expressed in the library of P. indicacolonized roots are difficult to detect and may not have significant impact on conserved miRNAs analysis. These results are similar to those reported from the plant pathogenic fungus Sclerotinia sclerotiorum [54], in which the number of low expressed conserved and specific miRNAs was also low. Furthermore, the number of conserved miRNAs in the control library was much lower than that in the library derived from P. indica-colonized roots. They were miRNA 397 (CK0328724; 6 reads), miRNA 171 (CK0347435; 6 reads) and miRNA 169 (CK0450655; 5 reads).
On the other hand, the number of reads of conserved miRNA families commonly occurring in both libraries showed significant differences between the two libraries. For example, the most abundant family mir528 generated 21287.27 TPM (times per million) in the P. indica-colonized library and only 3587.79 TPM in the control library. The top 20 miRNA families, including miR156 and miR528, were all significantly up-regulated in the P. indicacolonized library ( Table 2). In contrast, a few miRNA families such as miR169, miR396, miR319, mir160, and mir5648 were down-regulated in the P. indica-colonized library, when compared to the control library ( Table 2).

Identification of Novel miRNAs in Roots of 6 P. indicacolonized Oncidium Orchid
Since the whole genome sequence of Oncidium GR or any other orchid is not available, potential novel miRNAs were searched against the transcriptome database of Oncidium GR (http:// predictor.nchu.edu.tw/oogb/) [26]. The secondary structure of miRNA precursors and dicer cleavage sites were predicted using the Mireap software. 46 putative miRNAs and 51 precursors were detected. Among them, 7 putative miRNAs, including novel2, novel8, novel11, novel12, novel17, novel21 and novel44 were detected in both libraries. Interestingly, 3 putative miRNAs (novel9, novel10 and novel25) detected specifically in the P. indica-colonized library were identical in nucleotide sequence but different in their precursors. Similar results were found for novel11, novel12 and novel44 (Table  S3). These potential novel miRNAs were 20 to 23 nt long. Similar to conserved miRNAs, the majority (65.21%) of these novel miRNAs begins with ''U'' indicating that they were bound and processed by AGO1 [53]. The stem-loop structures of the potential precursors for the putative miRNAs were analyzed for RNA structure ( Figure  S4). The minimum free energy values ranged from 21.4 to 127.1 kcal/mol, and the average value was 49.99 kcal/mol.
The putative novel miRNAs showed much lower read numbers than those of the conserved miRNAs. 46 putative miRNAs generated only 175.09 TPM and 270.53 TPM in both libraries, respectively. The most abundantly detected putative miRNA was novel11, which was predicted to be generated from several different precursors and also named novel12 and novel44. In total, they were expressed with 50.88 TPM in the control library, but were downregulated with 3.71 TPM in the P. indica-colonized library. The predicted target gene was either CC-NBS-LRR (Unigene38669) or an rp1-like protein (JL929550.1) ( Table S3), suggesting that they may play a role in triggering plant immunity [55,56]. A second putative miRNA was novel9, which generated 49.43 TPM specifically in the P. indica-colonized library. It was also predicted to be generated from several different precursors, and also named novel10 and novel25. The predicted target gene coding for a protein kinase (JL915652.1 or JL916309.1) may be involved in a wide range of symbiotic processes [57] (Table S3).

Prediction of miRNA Target Genes
To understand the regulatory function of miRNAs during symbiosis, those miRNAs which were abundantly detected and significantly up-/down-regulated by P. indica were selected for further investigation. 1,834 conserved miRNAs belonging to 26 families were selected for target gene prediction. Due to the limited genome information, the information from the Oncidium mRNA database [26] was supplemented with mRNA database information from Oryza sativa and A. thaliana. After alignment in psRNAtarget, 702 best fit target candidates were obtained (Table  S4).
Subsequently, annotation and GO analysis were conducted by Blast 2 GO. 653 targets were annotated and distributed in 38 categories (Table S4). In the subcategory of biological process, targets genes were mainly involved in root development, defense, hormone-mediated signaling, cell death and cell cycle (Figure 2). In the subcategory of molecular function, target genes were mainly related to amino acid metabolism, secondary metabolites, kinases, epigenetic and posttranscriptional modifications, redox regulation and transporters (Figure 2). In the subcategory of cellular components, target genes belonging the mitochondrial respiratory chain complex I, heterotrimeric G-protein complex, integral membrane proteins of the endoplasmic reticulum and chloroplast proteins were detected (Figure 2). Combined with the number of reads (Table S5), the most abundant miRNAs were predicted to target the genes involved in root development, such as mir535, mir156, mir166, mir168, mir393, mir894 (Table 3). Others are related to auxin signaling, such as mir164 and mir167 (Table 3). Another group of miRNAs, such as mir528, mir397, mir408, may target genes involved in cell wall metabolic and redox regulation (Table 3). Subsequently, we chose some representative miRNAs to identify the predicted targets by surveying the Oncidium GR transcriptomic database (Table S6). The results indicated that most of the predicted target genes of the conserved miRNAs were consistent with the results using the O. sativa and A. thaliana mRNA databases, but some miRNAs such as mir528, which was only detected in monocotyledons, were predicted to target different genes in each of the three databases. For example, the targets of mir528 were predicted as FAR1 (NM_119978.3), DUF3444 (NM_001084639.1) and VTE1 (NM_119430.4) in the A. thaliana database; but they were predicted as plastocyanin-like protein (LOC_Os08g04310.1), a laccase precursor (LOC_Os01g44330.1), Cu/Zn-SOD protein (LOC_Os08g44770.1,2), multi-copper oxidase (LOC_Os01g03620.1) and L-ascorbate oxidase precursor (LOC_Os06g37150.1) [58][59][60][61][62][63] in the O. sativa database. However, a-lactate dehydrogenase (JL926279.1), pyruvate dehydrogenase e1 component subunit alpha (Unigene38362), F-box protein 13-like protein (Unigene13132) and a tetratricopeptide repeat-containing protein (Unigene39801) were predicted as target genes of mir528 in the Oncidium GR transcriptomic database. Similar results were found for mir390, which was predicted to target an LRR-kinase in Oncidium GR transcriptome database, but predicted as TSA3 in A. thaliana mRNA database.
In summary, the prediction of target genes revealed that the miRNAs in the P. indica-colonized roots of Oncidium may regulate a wide range of molecular processes. Especially, proteins involved in auxin signaling and development were identified (Table 3). Auxin has been shown to play a critical role in P. indica-mediated growth promotion and development in Chinese cabbage roots [14,64]. This implicates that miRNAs play important roles in the regulation of growth promotion during symbiosis of Oncidium and P. indica.

Expression Patterns of miRNAs and their Putative Targets in P. indica -colonized Roots
To elucidate the regulatory function of miRNAs on their putative targets, quantitative real time PCR (qPCR) was conducted to investigate the expression levels during different periods after colonization by P. indica. Most of the conserved miRNAs were detectable by qPCR ( Figure 3). Members in the same family showed similar expression pattern in qPCR, such as mir156s, mir535s and mir528s. Furthermore, in accordance with bioinformatic data of the deep sequencing, the most abundant miRNAs, such as mir528 and mir156b did not show higher expression level than other miRNAs. Conversely, mir894, mir535a, mir166a, mir166h, mir167, mir168, mir397 and mir2950 showed higher expression levels in P. indica infected Oncidium root tissue (Figure 3). A few conserved miRNAs with low read numbers can barely be detected in the qPCR assay. These results indicated that qPCR data are largely in accordance with the read number data from deep sequencing.
Subsequently, we used qPCR to study the expression pattern of miRNA target genes during the symbiotic process. Among the putative targets, mRNAs for auxin signal tranduction components showed significant changes (Figure 4). The mRNA for the auxin receptor protein TIR1, target of mir393, showed a gradual downregulation from 0 to 8 weeks. The mRNAs for auxin response factors (ARFs), targets of mir160, were up-regulated during the first week followed by down-regulation between the 3 rd and 8 th week. We could not detect the two ARFs, type 6 and 8, targets of mir167, in our qPCR assay, suggesting that the transcript abundance is low. ARFs are considered to be a key player in regulating the auxin signal pathway postively or negatively, and they are involved in auxin homeostasis by conjugating free auxin via GH3 [65].
Plants need a ''reprogram cycle'' to adjust growth and development during a symbiotic process [66]. The mRNA level involved in development, such as for the squamosa promoter binding protein, SPL, AGO1, and for the class III homeoboxleucine zipper transcription factor (PHV; HD-ZIP III transcription factors) were regulated by mir156/mir529, mir168 and mir166, respectively [67]. Their expression was down-regulated from the 1 st to the 3 rd week followed by up-regulation at the 8 th week ( Figure 4). Apparently, these mRNAs are first repressed and lateron stimulated. Another mRNA encoding the Nuclear Factor Y (NTFY) regulated by mir169, which is critical for symbiotic nodule development [48], was up-regulated during early colonizing periods and declined later (Figure 4).
Furthermore, cell wall metabolism-related genes, such as laccase (LAC) which is targeted by mir397, were up-regulated during the 1 st week, before they returned to the normal level between the 3 rd and 8 th week (Figure 4). Several LRR-kinase genes (LRR), which are targets of mir390, were down-regulated after colonization. F-box genes, targets of mir2950, were also down-regulated between the 1 st and 3 rd week followed by normal levels at the 8 th week (Figure 4). LRR-kinases [68] and F-box proteins [69] derive from large gene families and the proteins are involved in many different biological processes, thus the exact function of these genes during symbiosis needs to be further studied.
Considering that these targets were significantly up-or downregulated and the miRNAs were abundantly detected in the highthroughput sequencing, it is interesting to study the interaction between the miRNAs and their targets. Full understanding of the physiological function of these miRNAs and targets during symbiosis may help to unravel their roles in the symbiosis.

Discussion
Deep Sequencing miRNA in Roots of Oncidium P. indica colonizes the roots of a wide host range, thereby promoting growth and plant performance [7]. This is often associated with alterations in root architecture [13]. In this study, we demonstrate that the fresh weight of Oncidium seedlings was significantly increased after 8 weeks of co-cultivation with the fungus, which includes a higher leaf number, bigger stem diameter, an increase in the root number and root diameter and longer height of the seedlings (Figure 1). Due to the inherently slow growth of orchids, growth promoting effects have practical implications for commercial purposes. P. indica was reported to be isolated from the rhizosphere of the shrubs Zizyphus nummularia and Prosopis juliflora in the Indian desert [3]. It is closely related to the multinuclear fungi Rhizoctonia and Sebacinales, which are naturally associated with diverse orchid plants, and are capable in promoting orchid seed germination [29]. Here, we provide evidence that P. indica promotes growth of both aerial part and roots in orchids.
The purpose of this work is to understand the functions and regulatory mechanisms of miRNAs in Oncidium orchid growth regulation during the symbiosis with P. indica. A high-throughput sequencing and comparative expression analysis were conducted. In total, 17,036,953 unique sequences from 24,570,250 and 24,744,141 clean reads in control and P. indica-colonized root libraries were obtained, and 13,736 unique sequences showed homology to 941 miRNA families distributed in 45 plant species. Above all, there are 46 putative novel miRNAs and 51 corresponding precursors.
Orchids are ancient species which emerged shortly after the mass extinction. They became the most diverse plant family over the world [70]. To date, only a few studies on miRNAs in orchid were performed. They demonstrated that only a very few miRNAs are conserved between species. The majority of known miRNAs are family-or species-specific miRNAs, which derived from reversed duplication events or accumulation of mutations within inverted repeats. The differential expression of miRNAs may be essential for the physiological function and contributes to species evolution [71,72]. To date, miRNA studies have been carried out on the tropic orchid Phalaenopsis [73]. Compared to Oncidium, the most abundant top 20 miRNAs conserved in roots, such as mir528, mir156, mir166, mir167, mir168, mir159, mir529, and mir894, were also abundantly found in Phalaenopsis. The Oncidiumspecific miRNAs are mir397, mir408, mir2916, mir2950, mir5077, and mir5059, while Phalaenopsis-specific miRNAs are mir169, mir172, mir396, mir1318 and mir2911. Furthermore, we also detected a considerable number of new members of conserved miRNA families and novel miRNAs. It is known that the variation in miRNA precursors foldback, processing efficiency and miRNA size may influence the maturation type of miRNAs. The young miRNAs are often expressed in very low levels and in restricted spatial/temporal patterns, and some of the newly identified miRNAs have no obvious function [72]. Most of the newly identified miRNAs in this study had low number of reads and were hardly detectable by qPCR. Further research is required for eludicating their function in Oncidium.

Role of miRNAs and their Targets on Regulating Phytohormone Levels and Root Development during Symbiosis
The miRNAs involved in auxin signalling functions were easily detectable in our experiment, including miR160, mir164, mir167, mir393, miR394 and mir5293. Besides, mir390, another auxinrelated miRNA is related to the Arabidopsis TAS3 RNAs, and predicted to target the mRNA encoding a LRR-kinase, which is responsible for regulating lateral root growth [74]. This miRNA is common in monocotyledonous plants [75]. Phytohormones play key roles in plant growth and development, and are also required for resistance against abiotic and biotic stresses [11,12,76,77]. Previously, phytohormone metabolism genes including auxin, gibberellin (GA) and abscisic acid (ABA) were found to be upregulated after colonization with P. indica in barley roots [11], and auxin was considered to be a key phytohormone in plant growth Table 3. Target genes prediction for miRNAs in P. indica colonized-Oncidium roots. promotion triggered by P. indica [14,64]. In our experiments, auxin-related genes such as TIR1, ARF and CUC, which were regulated by mir393, mir160 and mir164, were up-or downregulated during different phases of the symbiosis, and a putative GA pathway gene, MYB, regulated by mir159, was downregulated during symbiosis. Auxin and GA probably function as a SA suppressor or JA enhancer in symbiotic processes [11,[78][79][80]. Taken together, the high expression level of this group of miRNAs may probably balance growth and defense during symbiosis. Another group of miRNA including mir156, mir166 and mir169, which targets mRNAs for the transcription factor genes SPL, PHV/PHB and NTF, were also abundantly detected in our experiment (Table S5). Also mir162 and mir168 are highly expressed. They function as feed-back regulators to miRNAs and catalyse the AGO1 and DCL mRNAs (Table S5). AGO1 and the transcription factor PHV/B may cooperatively control the development of shoot apical meristem (SAM) and the root apical meristem (RAM) [81]. Plants which fail to express these miRNAtargets have lesions in root development [66,[81][82][83][84][85][86][87]. Furthermore, these miRNAs also accumulated in the nodules of Medicago truncatula after colonization by the endophytic fungus [48,49]. It appears that plants reprogram the development of the SAM and RAM in the symbiotic interaction with microbes by regulating the expression of these miRNA in critical spatial and temporal patterns [48].
MiRNAs such as mir397 or mir408 may be necessary for Oncidium to enhance anti-oxidative capacity and resistance against stress by mediating copper distribution and playing roles in the cell wall metabolism during growth and development processes trigged by P. indica.

Conclusions
We report on the miRNA profiling of Oncidium orchid after colonization by P. indica. The miRNAs and their target genes illustrate that the physiological metabolism of Oncidium is reprogrammed in response to the symbiotic interaction. Genes participating in phytohormone signaling, for cell wall metabolism and regulatory transcription factors are major targets of the P. indica-induced miRNAs in Oncidium roots. Therefore, we propose that P. indica alters the miRNA pattern to establish an intricate network for growth promotion, developmental reprogramming and enhances resistance in the roots of Oncidium. Several novel unique miRNAs were detected, for which a function could not yet be identified. Further investigations on the molecular mechanism of miRNAs in symbiotic interactions are of huge significance.

Plant Materials
Oncidium GR flask seedlings were cultured on 1/2 Murashige-Skoog (MS) medium, and incubated with 16 h photoperiod (100 mmol m 22 s 21 ) at 25uC. Seedlings about 3 cm high were selected and transferred to fresh medium. After 10 days of acclimating cultivation, plants were inoculated with P. indica mycelium as described [52]. Briefly, one agar block of 5 mm in diameter was placed per seedling at a distance of 1 cm from the adventitious root. The seedlings were either exposed to an agar block containing P. indica (+ P. indica) or to a block without the fungus (control; -P. indica). After 8 week of co-cultivation, the biomass of the seedlings was determined. They were then immediately frozen in liquid N 2 and stored at 280uC for RNA extraction and other assays. Anatomic Dissection of Root Tissues Colonized by P. indica Root samples were fixed in 4% formaldehyde in 20 mM PBS, pH 7.4, for 60 minutes at 37uC followed by dehydration with 5% sucrose in 20 mM PBS, pH 7.4 at 4uC overnight. Thin sections were cut by free hand and stained with chitin-specific WGA, Alexa Fluor R 488 conjugate (WGA-AF488, Invitrogen USA) or Alexa Fluor R conjugated to succinylated concanavalin A (ConA-AF633, Invitrogen USA), according to the manufacturers protocol. Sections were photographed by Olympus IX17, WGA-AF488 was excited at 488 nm and detected at 505-540 nm, ConA-AF633 was excited at 546 and detected at 600-660 nm.

Small RNA Library Construction and Sequencing
Total RNA was isolated from control roots and those colonized with P. indica using the pine tree method. The RNA integrity number (RIN) was examined on Agilent 2100; when the RIN was .8.0, the samples were sent to BGI (Shengzheng China) for small RNA Solexa sequencing. In brief: the small molecule RNAs were separated by 15% (w/v) PAGE (18-30 nt), subsequently, the purified small RNAs were ligated to a pair of Solexa adaptors to the 59 and 39 ends, reverse transcribed to cDNA using a RT primer, and finally amplified by PCR and sequenced. The results were deposited in Sequence Read Archive (SRA) on NCBI website. The accession number is SRP031471.

Identification of Conserved and Novel miRNAs
After removal of the adapter sequences, low quality tags and contaminated sequences, rRNA, tRNA, snRNA, etc. were identified by alignment to the Rfam 10.1 and Genbank databases. Conserved miRNA were identified by search against known plant miRNAs in miRBase 19.0 (http://www.mirbase.org). Novel miRNAs were predicted by Mireap (http://sourceforge.net/ projects/mireap/). Secondary structure of all candidate precursors was mapped in RNAstructure Web (http://rna.urmc.rochester. edu/RNAstructureWeb/).

Prediction and Annotation of Potential Targets of miRNAs
The potential targets of miRNAs were predicted with psRNATarget (plantgrn.noble.org/psRNATarget/) by alignment to mRNAs from Oncidium, Oryza stativa and A. thaliana (data were downloaded from the NBCI GenBank), the best hit results were chosen and annotated using Blast 2 GO.

Expression Analysis of miRNAs and their Potential Targets by Quantitative Real Time PCR
For expression analysis of miRNAs and their potential targets, RNA samples were isolated from roots, which were collected from 3 independent replicates of seedlings cultured under the same condition as described above. Small and high molecular weight RNAs were isolated by LiCl and isopropanol methods, respectively [73]. For miRNA quantitative real time PCR (qPCR), the cDNAs were synthesized with the QuantiMir RT kit (RA420AU-3; SBI). Abundantly expressed miRNAs (based on the number of reads) and novel miRNAs were chosen for analysis ( Figure S2).
Using 5.8S rRNA as internal control, expression analyses of these miRNAs were performed with the Applied Biosystems 7500 Fast Real-Time PCR System using the SYBR Green PCR master kit (with ROX). The qPCR reaction was performed in a 30 ml reaction system and programmed as follow: 95uC for 2 min; 40 cycles of 95uC 15 s, 65uC 15 s, 72uC 30 s, following by a dissociation program from 60uC to 95uC. Each reaction was repeated three times.
For potential target qPCR, the cDNAs were synthesized by Revert Aid First Strand cDNA Synthesis kit (#k1622; Fermentas). Using actin as internal control, the qPCR reaction was carried out in a 20 ml reaction system, and the qPCR program is identical to that described above, but annealing was performed at 60uC. Primers information is shown in Table S1.