The YABBY Genes of Leaf and Leaf-Like Organ Polarity in Leafless Plant Monotropa hypopitys.
Monotropa hypopitys (syn. Hypopitys monotropa) is a member of the flowering seed plant family Ericaceae, which in turn belongs to the order Ericales splitting from the base of the clade Asterids . This mycoheterotrophic, nonphotosynthetic, achlorophyllous plant acquires carbon from the roots of autotrophic trees through monotropoid mycorrhizal symbiosis [2, 3]. The M. hypopitys root system consists of fleshy roots, on which shoot buds develop, and finer mycorrhizal roots . This plant has a typical aboveground structure, although the stem and leaves can be taken for flowering parts--floral axis and sterile bracts [4, 5]. Similar to extant plants, M. hypopitys forms asymmetrical lateral organs on the flanks of a shoot or inflorescence apical meristem, with adaxial and abaxial surfaces adjacent to or distant from, respectively, the meristem. Paleobotanic studies indicate that such structural asymmetry first appeared in a true leaf (euphyll) transformed from a radially symmetric stem as a consequence of the need to absorb more sunlight [6-9]. Studies of asymmetry in plants indicate that during plant development, cell fate is determined mostly by positional signals. The correct maintenance of the apical meristem and abaxial-adaxial differentiation of lateral organs requires reciprocal signal interaction between the meristem and derived structures [10-12].
It has been shown that the polarity of leaves and floral organs is defined by the network of genes encoding the Class III Homeodomain Leucine Zipper (HD-ZIPIII), ASYMMETRIC LEAVES (AS1/AS2), KANADI (KAn), AUXIN RESPONSE FACTOR (ARF3/ARF4), and YABBY families of transcription factors [8, 13-15]. Among them, the HD-ZIPIII REVOLUTA (REV) is expressed in the adaxial domain of lateral organs, whereas the GARP-family transcription factors KAN1-4 are involved in abaxial differentiation. Together, REV and KAN1 antagonistically regulate the expression of a number of genes encoding auxin signaling and transport components [16, 17].
The YABBY genes originating from the lineage leading to seed plants are identified in various green spermatophyte plant species; they are closely associated with the evolutionary emergence of flat-shaped leaves and are presumably diversified during evolution, which resulted in the appearance of family members with specific functions in the leaf, carpel, and ovule [8, 18-20]. YABBY transcription factors are characterized by their nuclear localization and the presence of the [C.sub.2][C.sub.2] zinc-finger and DNA-binding YABBY (High Mobility Group- (HMG-) box-like) domains [21, 22]. In gymnosperms, the YABBY genes are distributed among the A, B, C, and D clades . In extant angiosperms, five YABBY subfamilies, FILAMENTOUS FLOWER (FIL or YAB1)/YABBY3 (YAB3 or AFO), CRABS CLAW (CRC), INNER NO OUTER (INO or YAB4), YAB2, and YAB5  are distinguished by conserved functions in the initiation of lamina outgrowth, polarity maintenance, and establishment of the leaf margin [23-25]. Almeida et al.  and Morioka et al.  have provided evidence for the involvement of the YAB2 and YAB5 genes in the evolutionary diversification of style and filament morphology. The branching of INO and CRC from other YABBY genes has most likely occurred in parallel with the evolution of the carpel and outer integument via modification of reproductive leaf-like sporophylls [20, 28].
There are current theories related to the history of the YABBY genes in angiosperms. Bartholmes et al.  suggested that "vegetative" YABBYs (FIL/YAB3, YAB2 and YAB5) do not form a monophyletic clade and that CRC and FIL evolved from a common ancestor gene, while the INO genes are sisters to that ancestral gene. On the other hand, Finet et al.  clustered INO together with clades YAB5 and YAB2. In addition, two alternative evolutionary scenarios, that is, monophyly or paraphyly of the gymnosperm YABBY family towards angiosperm YABBY genes, suggest that all spermatophyte YABBY genes were derived from one or two, respectively, YABBY genes of the last common ancestor of extant seed plants [20, 23, 29]. The reconstruction of YABBY evolution in spermatophytes based on these theories suggests that at least one YABBY predecessor has already functioned as a polarity regulator and that the diversification of the gene family occurred in both angiosperms and gymnosperms . Although the presence of the YABBY genes is presumably restricted to seed plants [19, 20], genomics studies conducted on marine picoeukaryotes revealed YABBY homologs in Chlorophyta . Phylogenetic analysis of identified sequences suggested, with equal probability that either Chlorophyta YAB genes are evolutionarily related to seed plant YABBYs or emerged independently from ancestral HMG-box sequences .
The expression data available for the angiosperm YABBY genes suggested that the FIL-, YAB2-, and YAB5-like genes retained a more ancestral expression pattern in both vegetative and reproductive tissues, while the expression of the CRC and INO-like genes is more variable . In Eudicots, FIL, YAB3, YAB2, and YAB5 transcripts were detected in the abaxial side of primordia in all aboveground lateral organs (except ovules) determining the abaxial cell fate [31, 32]. The CRC genes are expressed abaxially in the carpel, placenta, and nectaries promoting the development of the gynoecium and abaxial part of the carpel wall, and terminating the floral meristem [33-38]. INO mRNA is detected in the abaxial epidermis of the outer integument [20, 28, 39, 40]. The YABBY expression pattern differences between cereal monocots and other angiosperms indicate the modification of genetic pathways involving YABBYs during the process of angiosperm diversification [21, 27, 41-46]. It is assumed that FIL, together with REVOLUTA (REV), APETALA1 (AP1), and LEAFY (LFY), corrects the spatial activity of the AGAMOUS (AG), AP3, PISTILLATA (PI), and SUPERMAN (SUP) genes, and, thus, is involved in the initiation of floral organ primordia at the correct position and numbers, defining the fate of appropriate cells [31, 44, 47].
Thus, the bifunctional YABBY transcription factors have an important role in driving the evolution of the leaf and gynoecium, as well as in the initiation, growth, and structural organization of almost all aboveground lateral organs, and in the control of shoot apical meristem organization and activity.
In the present study, we identified and phylogenetically classified seven YABBY members from M. hypopitys and characterized their expression profiles in various tissues during flowering. The structural features and composition of conserved motifs belonging to the predicted MhyYABBY proteins were also analyzed. Our data should further the understanding of possible links between polarity determination and the physiology of achlorophyllous mycoheterotrophic plants.
2. Materials and Methods
2.1. Plants and Transcriptomes. The previous study divided M. hypopitys specimens into a North American cluster and two Eurasian (excluding Russia) sister lineages (Swedish and pan-Eurasian) . Analysis of M. hypopitys from the European part of Russia revealed two types, A and B, which showed 99 and 100% homology with specimens from Japan, Finland, and Great Britain, and with Swedish specimens, respectively [49-51]. The study of H. monotropa specimens from Northern Ireland showed that they occur in small, highly fragmented populations, and exhibit a relatively high level of within-population genetic diversity and a low level of clonality .
In the present study, two M. hypopitys plants of type B from one clone of the same genet were used. Flowering plants were collected in a coniferous forest, Kaluga region, Russia, in August, 2015. The individual plant was a 15 cm reproductive axis with bracts and raceme of 10-12 flowers (each of 4 sepals, 4 petals, 8 stamens, and 4 fused carpels), and root system comprised mycorrhizal and fleshy roots with adventitious buds. The annual floral axes arise from adventitious buds on the perennial roots and carry the laminar appendages (there are no flowers in their axils, but they are above the soil level), which are termed sterile bracts (in our study, bracts, for simplicity) .
Plants were dissected into flowers, bracts, fleshy roots with adventitious buds, and predominantly haustoria-enriched roots, immediately frozen and homogenized in liquid nitrogen, and stored at -80[degrees]C. Total RNA was isolated from tissue of each M. hypopitys bracts (two individual plants), flowers (two individual plants), roots containing buds (individual plant), and haustoria-enriched roots (individual plant) and used for mRNA library preparation, which was sequenced on the Illumina HiSeq2500 platform (Illumina Inc., San Diego, CA, USA). The M. hypopitys RNA-seq data for each of six transcriptomes were assembled into the 98,350 unigenes with a length of 201-12,993 bp [51,53]. Individual reads were mapped on contigs using Bowtie 2 , and protein-coding genes in contigs were identified using TransDecoder (https://transdecoder.github.io/).
2.2. Identification and Bioinformatics Characterization of M. hypopitys YABBY-Coding Sequences. To identify the M. hypopitys genes homologous to the known organ polarity genes, we searched unique transcripts revealed by the RNA-seq against the NCBI database (http://blast.ncbi.nlm.nih. gov/). To predict YABBY transcripts in M. hypopitys, we additionally searched the assembled transcriptomes with the known YABBY-related sequences coding for conserved zinc-finger and HMG-like domains extracted from the NCBI database. Selected YABBY candidates were examined for open reading frames (ORFs), translated using Clone Manager v.7.11 (http://clone-manager-professional.software.informer .com/), and the conserved domains of putative MhyYABBY proteins were identified using the NCBI-CDD analyzer (http://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi) and specified according to .
To evaluate the overlap between transcriptomes, Venn diagrams were generated using the online program Venny . To illustrate the transcriptome-based gene expression pattern, the data were clustered with the Average linkage method and Spearman rank correlation (distance measurement method), and visualized as a heat map (http://www2. heatmapper.ca/) .
Conserved MhyYABBY amino acid motifs were identified using the MEME (Multiple Expectation Maximization for Motif Elicitation) 4.11.2 online analysis (http://meme-suite. org/tools/meme)  and used to construct a schematic diagram. To search for motifs, the "Normal" motif discovery mode, the default width range of 6-50 amino acids (aa), and the motif site distribution "zero or one per sequence" were used. The identified motifs were manually compared with previously suggested specific motifs. Initial search was performed on a set of 34 complete sequences, including identified MhyYABBYs, independently of YABBY clade affiliation. In addition, variable regions between conserved domains within the same group of proteins were searched. Since most of the identified motifs were clade-specific, we further analyzed individual YABBY clades.
To investigate the evolutionary relationship of the MhyYABBY genes, the MhyYABBY proteins and YABBY homologs from other species available in NCBI were aligned using ClustalX . For analysis, full-size amino acid sequences, as well as conserved regions consisting of the zinc-finger and HMG-like domains were used. Evolutionary divergence between the genes and proteins was estimated using the maximum composite likelihood and equal input models, respectively, in MEGA7 [59-62]. The phylogenetic tree topology was estimated using the maximal likelihood method based on the JTT matrix-based model in MEGA7 .
2.3. Analysis of Tissue-Specific Gene Expression. The MhyYABBY gene expression was calculated in each transcriptome. Transcript quantification based on RNA-seq data was performed without a reference genome using the RSEM  and Bowtie 2  programs, including normalization of transcripts per kilobase of exon per million fragments mapped (FPKM) values and between samples.
To perform quantitative real-time PCR (qRT-PCR), the first strand cDNA was synthesized from 1 [micro]g of each mixture of two-root, two-bract, and two-flower RNA preparations using the Reverse Transcription System (Promega, Madison, WI, USA) and an oligo-dT primer, and quantified using the Qubit[R] Fluorometer.
Based on the identified YABBY-like transcripts and corresponding draft genomic sequences (our unpublished data), gene-specific primers separated by at least one big intron were designed to amplify parts of gene-coding sequences (Supplementary Table 1). The qRT-PCR was performed in three technical replicates using 2.5 ng of cDNA and an SYBR Green and ROX RT-PCR mixture (Syntol, Moscow, Russia) at the following cycling conditions: initial denaturation at 95[degrees]C for 5 min, 40 cycles of denaturation at 95[degrees]C for 15 sec, and annealing/synthesis at 60[degrees]C for 40 sec. Obtained PCR-fragments were additionally purified and sequenced to confirm certain gene specificities. Gene expression levels were normalized to those of the reference pinesap Actin5, Actin3, and SAND genes (primers are provided in Supplementary Table 1), which transcripts were evenly represented in six transcriptomes . Normalized expression data were statistically evaluated using GraphPad Prism version 7.02 (San Diego, CA, USA; https://www.graphpad.com/scientificsoftware/prism/). Three values (3 technical replicates) were used for SD calculation. The error bars were generated based on mean with SD calculation. Significance of the qRT-PCR data within the same tissue between species was estimated by unequal variance Welch's t-test and additionally treated with Bonferroni's correction: if any of the t-tests in the list had p [less than or equal to] 0.05/number of t-tests in the list, then the null hypothesis was rejected; that is, the difference between samples was recognized as significant.
3.1. M. hypopitys Organ-Polarity Genes. The obtained six transcriptomes  showed significant overlap between the whole set of reads of paired libraries (74-77%), and in three libraries for each of the two plants (>80%) (Figure 1), and reflected a number of the known plant organ polarity genes that are being expressed in different tissues of flowering pinesap (Supplementary Table 2). Among them, the genes of the KAN and REV transcription factors, responsible for the abaxial and adaxial cell identity in lateral organs, respectively, and their common targets (ARF, SAUR, Aux/IAA, PID, PIN, NPY, GH3, etc.) involved in auxin biology  were identified. The genes encoding SEUSS (SEU) and LEUNIG (LUG) involved in petal polarity determination along the adaxial/ abaxial axis , AP2-like transcription factor AINTEGUMENTA (ANT) which contributes to organ polarity , the ADP ribosylation factor guanine nucleotide exchange factor GNOM essential for basal polarity establishment in A. thaliana , ULTRAPETALA1 (ULT1) which acts antagonistically with KAN1 to pattern the adaxial-abaxial polarity axis but jointly to pattern the apical-basal axis restricting the expression domain of the SPATULA gene , and HDZIPIII transcription factors HAT and ATHB positively regulated by REV  were also found (Supplementary Table 2). Heat map-based clustering of transcriptomic reads associated with organ polarity genes revealed similarities between pair libraries (according to the column dendrogram in Figure 2(a)). For most genes, the expression levels in flowers were higher than those in bracts (Figure 2(b)).
All six transcriptomes of bracts, flowers and roots contained seven unique YABBY-like transcripts, which were considered as putative MhyYABBYs (MhyYAB1-MhyYAB7). The size of putative MhyYABBY ORFs varied from 502 to 673 bp, and the length of predicted proteins was from 166 to 224 aa. One transcript, MhyYAB4, had a partial 3-truncated coding sequence (CDS), while the remaining six mRNAs contained complete CDSs. Putative MhyYABBY proteins included both a conserved N-terminal 37-aa [C.sub.2][C.sub.2] zinc finger-like domain and a C-terminal 48-aa helix-loop-helix domain resembling a part of an HMG box (Figure 3) [21, 33]. A cluster of amino acids at the beginning of the HMG-like domain could potentially serve as a nuclear localization signal [33, 69].
Among MhyYABBYs, we distinguished two groups (comprising MhyYAB2, MhyYAB5, and MhyYAB6; and MhyYAB3 and MhyYAB4) based on high identity outside the conserved domains (Supplementary Figures 1a, 1b). The estimated low evolutionary pairwise divergence in the nucleotide and amino acid sequences between the members of each group compared to that between the members of different groups suggested the presence of two sets of paralogs (Supplementary Tables 4 and 5). MhyYAB1 showed high pairwise divergence with other MhyYABBYs, which indicates the affiliation of MhyYAB1 to the separate clade.
3.2. Phylogeny of the M. hypopitys YABBY Family. Previous phylogenetic studies considered only conserved YABBY domains; however, it has also been shown that variable protein regions contain clade-specific conserved motifs of potential functional importance . In transcription factors, variable regions are often essential for their activity and/or formation of multimeric protein complexes, as it has been shown for MADS-box transcription factors . We aligned complete amino acid sequences and only conserved domains of MhyYABBYs. For the group comprising MhyYAB2, MhyYAB5, and MhyYAB6 paralogs, two different results within the same clade were obtained. The full-size proteins were orthologous to FIL, while the conserved domains were orthologous to YAB3, another member of the FIL/YAB3 clade. We decided to use complete sequences to increase the sensitivity of phylogenetic analysis.
The generated tree, rooted with the Micromonas commode (Chlorophyta) YABBY-like protein, classified MhyYABBY1-7 by comparing them with YABBY-like proteins of A. thaliana and other angiosperm species belonging to Asterids and Rosids. As a result, MhyYABBY transcripts were distributed into the FIL, INO, CRC, and YAB5 clades and renamed accordingly as MhyCRC (MhyYAB1), MhyINO1 and MhyINO2 (MhyYAB3 and MhyYAB4, resp.), MhyYAB5 (MhyYAB7), and MhyFIL1, MhyFIL2, and MhyFIL3 (MhyYAB5, MhyYAB6, and MhyYAB2, resp.). It should be noted that no YABBY2-like transcripts were found in all six transcriptomes. The sequences have been deposited in GenBank (KX12839-KX12841, KX12843-KX12846; Supplementary Table 3). The maximum likelihood reconstruction of the MhyYABBY family with bootstrap values at tree nodes is shown in Figure 4. Within the individual clades, MhyYABBY sequences are sisters to other YABBY-like sequences from Asterids, and the closest homologs are YABBYs from Ericales.
3.3. MhyYABBY-Specific Motifs Identified outside of the YABBY Domains. To further investigate the structural divergence of pinesap YABBY proteins, we searched for the motifs conserved within individual clades or the whole MhyYABBY group. Comparison of complete sequences of YABBY orthologs (Supplementary Table 6) revealed two major domains, zinc-finger and YABBY, specific to all YABBY proteins. The number of motifs in variable regions was 4 to 9 with the length from 5 to 26 residues. The obtained data were compared with previously defined motifs specific to individual YABBY clades .
In MhyINO1/2, the known INO-A motif was found immediately after the zinc-finger domain. Comparison of MhyINO1/2 with other INO-like sequences in the NCBI database revealed novel putative INO-specific motifs. At the N-terminus and between the zinc-finger and YABBY domains, a highly conserved 11-aa INO-B motif and an Eudicot-specific 15-aa INO-C motif were identified. In addition, we found a C-terminal 26-aa INO-D motif specific to MhyINO1/2 and INO-like proteins in Solanaceae. All previously predicted CRC-specific motifs, CRC-A, CRC-B, CRC-C, and CRC-D, were found in MhyCRC. In addition, we identified a putative conserved Eudicot-specific C-terminal 15-aa CRC-E motif. In MhyYAB5, the presence of the known YAB5-specific motifs YAB5-A, GY/YAB2/5-A, YAB5-B, and GY/YAB2/5-B was confirmed. The YAB5-D motif was not detected, while a modified YAB5-C sequence was identified before the YABBY domain as a 12-aa YAB5-Cm motif. In MhyFIL1/2/3, we found all previously predicted motifs, FIL-A, FIL-B, FIL-C, FIL-D, FIL-E, FIL-F, and FIL-G. In addition, we suggested two motifs, 12-aa FIL-I (monocot and eudicot-specific) and 11-aa FIL-H (eudicot-specific) located between the zinc-finger and YABBY domains, as candidates for conserved FIL-specific motifs.
Thus, we characterized MhyYABBY proteins by sequence conservation on the clade level. The motifs identified in MhyYABBYs were shown in Figure 5 and sequences of the predicted novel motifs are presented in Supplementary Figure 2).
3.4. Expression Pattern of "Vegetative" and Flower-Specific M. hypopitys YABBY Genes. The transcriptome-based data on MhyYABBY expression in the roots and buds, haustoria-enriched roots, bracts, and flowers (Figure 6, Supplementary Table 3) showed that the MhyYABBY mRNAs (except MhyFIL2) were present in the flowers. Except for INO- and CRC-like MhyYABBYs, all other five transcripts were detected in the bracts. The highest bract-specific expression was observed for MhyYAB5, while the remaining four genes were transcribed at similarly low levels. Finally, in roots and buds, only MhyYAB5 and MhyFIL2 mRNAs were expressed at very low levels, while in the haustoria-enriched roots none of genes were expressed. The MhyFIL1 and MhyFIL3 transcripts were increased from the bract to the flower, maintaining the same expression profile. In contrast with this, the number of MhyYAB5 transcripts was decreased from the bract to the flower.
Quantitative (q) RT-PCR data on MhyFIL3, MhyYAB5, MhyINO1, MhyINO2, and MhyCRC expression are represented at Figure 6(d) and Supplementary Table 7. The relative expression of MhyYABBY genes was estimated in the flowers, bracts, and roots and buds. All analyzed genes were expressed in the flowers with the highest MhyYAB5 level, and the lowest MhyINO2 and MhyFIL3 levels. In the bracts, the MhyYAB5 gene was also highly expressed, but only traces of the MhyINO2 and MhyFIL3 mRNAs were observed. In the roots and buds, the only MhyYAB5 mRNA was detected at low level. The difference in the MhyYAB5 gene expression between the pinesap tissues was statistically significant. The flower-specific expression of MhyCRC and MhyINO1 was significantly different from that in bracts and roots. The expression of MhyINO2 and MhyFIL3 was similar between tissues (Supplementary Table 7). All the analyzed gene expression modes were the same as it was shown in their transcriptome-based patterns, except for MhyINO2. In the flower, the measured qRT-PCR expression of this gene was equally low. Given transcriptomic data, MhyINO2 was absent in flower 2, which may be due to the quality of the libraries or their sequencing. Also, considering the low evolutionary pairwise divergence in the INO1 and INO2 sequences, a higher level of INO1 expression compared to the level of INO2 may indicate that INO1 may be more required than INO2 during plant development.
The emergence of photosynthesis has become the most significant event in the evolution of plants. The majority of extant plants are autotrophic, except for about 1% of flowering heterotrophic plants. Among the latter, obligate mycoheterotrophs are the results of replicated deevolutionary events of the photosynthetic ability loss, triggering the degradation of both cytoplasmic and nuclear genomes . Full mycoheterotrophs demonstrate a wide range of deevolutionary outcomes such as abrupt morphophysiological changes , genome rearrangements, and massive gene loss [71, 72].
In the large and diverse eudicot family Ericaceae with a nearly worldwide distribution, two of nine subfamilies, Pyroloideae and Monotropoideae, contain partial and full mycoheterotrophs, respectively . In Monotropoideae, M. hypopitys represents a unique obligate mycoheterotroph. Recent studies on the M. hypopitys plastid genome and its comparison with that of photosynthetic relative Pyrola rotundifolia indicated that this plant is at the final stages of plastome degradation, which is expressed in highly reduced size and content, dramatic structural rearrangements, and acceleration of nucleotide substitutions in all protein-coding genes [74-76]. Furthermore, the coordinated loss of photosynthesis-related functions in both plastome and nuclear genomes of M. hypopitys is a sign of ongoing changes in the nuclear genome of this mycoheterotrophic plant .
It is generally accepted that mycoheterotrophic plants have evolved from photosynthetic mycorrhizal lineages, as mycoheterotrophy helps to succeed in the low-light conditions of the forest . It has been established that dark-induced leaf senescence leads to a significant chlorophyll loss and photosynthesis inactivation [78, 79]. During evolution, a M. hypopitys ancestor (already with megaphylls) growing in shaded habitats lost the genetic ability to photosynthesize due to symbiosis with fungi, which provided a sufficient amount of carbon to pinesap from the roots of autotrophic trees. It is shown that the loss of photosynthetic ability and full heterotrophy are linked to the degradation and/or modification of vegetative structures . It is believed that in M. hypopitys, an elongated raceme emerges instead of a true stem, developing directly from the adventitious bud on the roots .
The photosynthetic ability is closely related to the origin of asymmetrical leaves (providing the absorption of sufficient light energy by seed plants [80, 81]), in particular, due to the YABBY genes' evolutionary duplication and diversification [6-9, 25, 82, 83]. Although the role of the YABBY genes in plant evolutionary adaptation to light perception is established and they have been systematically studied in model and nonmodel species [27, 84-87], up to now, no YABBY genes have been described in mycoheterotrophic plants. It was interesting to figure out, if these genes and, therefore, the conserved mechanism of leaf polarity determination, were exposed to the adaptive deevolution in leafless mycoheterotroph M. hypopitys. Therefore, in this study we focused on the diversity and expression profile of the YABBY genes in a M. hypopitys that may further the understanding of the development and evolution of this plant group.
Peripheral cells of the shoot apical meristem give rise to the leaves that develop along three axes and acquire the adaxial-abaxial and proximal-distal asymmetry and the mediolateral symmetry . The elongated M. hypopitys raceme carries bracts below flowers and leaf-like sterile bracts . Bracts are thin, 8-15 mm long, 3-15 mm broad, ovate, and expanding to the top, with irregularly toothed margins , rudimentary midvein, and parallel veins of similar thickness (Figures 7(a) and 7(b)). Interestingly, sterile bracts are not exactly leaves as they have the genetic signatures of reproductive organs, which are manifested in the expression of the floral organ identity MADS-box genes  (Figure 2).
The lack of leaves in M. hypopitys may correlate with possible changes in the conserved genetic network of lateral organ polarity. In pinesap transcriptomes, we found the number of genes associated with this network, including seven YABBY-like sequences encoding proteins, which contain conserved domains and nuclear localization signals characteristic for YABBY transcription factors.
Initially, "adaxial" and "abaxial" genes are expressed throughout the leaf primordium, and, as the leaf develops, their expression becomes restricted to their respective domains due to the mutually exclusive actions of their protein products . YABBYs are "abaxial" genes involved in stimulation of the cellular division during lamina outgrowth in all aboveground lateral organs, vegetative or reproductive (e.g., [25, 91-93]). M. hypopitys flowers do not show any visible abnormalities compared to those of other eudicots. Thence, MhyYABBYs may play common roles in the proper development of the floral meristem into a mature flower. However, the question arises how did the lack of leaves affect the function of the "vegetative" MhyFIL1-3 and MhyYAB5 genes.
In A. thaliana leaves and sepals, FIL and YAB3 genes are upregulated by KAN1 and ARF4, and, in turn, FIL and YAB3 stimulate the expression of ARF4, KAN1, and AS1 [92, 94, 95], and besides, in complex with LUG and SEU promote not only organ polarity, but embryonic shoot apical meristem initiation and maintenance . "Abaxial" KAN represses the "adaxial" HD-ZIPIII genes [90, 96, 97]. At the boundary between adaxial and abaxial tissues, the FIL/YAB3 and KAN, respectively, up- and downregulate the WOX1 and WOX3 genes that specify redundantly lateral lamina outgrowth and leaf margin cell fate . The YAB5, KAN2, ARF3, and ARF4 genes are repressed by the "adaxial" AS1 and AS2 implicated in the proper leaf formation along all three axes [99-101].
During flowering, YABBYs are required to establish a correctly developed flower primordium through the interaction with REV, KAN4, SEU, LUG, ANT, SUP, LFY, and the floral homeotic MADS-box genes [28, 47, 64, 65, 91, 93]. The SEU and LUG are needed to promote and maintain the FIL/YAB3 and HD-ZIPIII expression . In turn, FIL in combination with ANT acts to upregulate the "adaxial" gene PHB and MADS-box gene AP3 , and together with REV, AP1, and LFY, spatially regulates the transcription of the SUP and the MADS-box genes AG, AP3, and PI [31, 44, 47]. To maintain the polar development of the ovule outer integument, the INO interacts with SEU and LUG, but its expression is restricted by KAN4, REV, and SUP . CRC, upregulated by AP3/PI/SEP, is involved in the control of radial and longitudinal gynoecium growth, carpel fusion, and nectary location, and participates in the floral meristem termination through the WUS repression [33, 35, 36, 102].
The finding of almost all the above-described polarity genes in the analyzed transcriptomes (Supplementary Table 2; Figure 2) suggests that the polarity of M. hypopitys bracts and floral organs is under the control of conserved mechanisms (Figure 7), with the exception of the absence of transcripts PHB and PHV (HD-ZIPIII), WOX1 and WOX3 (homeodomain protein), and AS2 (LBD domain transcription factor). The PHB and PHV function may be replaced by another member of the HD-ZIPIII family, REV, since these three genes can function redundantly . Similarly, the MhyWOX genes may perform the functions of WOX1 and WOX3 . Although 12 AS2-like (ASL) genes were found in M. hypopitys transcriptomes, the AS2 cannot be functionally replaced by other family members . The lack of AS2 transcripts may contribute to a high level of MhyYAB5 expression compared to other "vegetative" MhyYABBYs, since the AS1/AS2 complex suppresses YAB5 . Moreover, the lack of AS2 activity may be related to the characteristics of M. hypopitys bracts having a plump lamina base, the midvein indistinguishable from parallel veins, and the absence of petioles. This conclusion may be supported by the phenotype of the mutant as2, which has a significantly rudimentary leaf midvein (and several parallel veins of very similar thickness), shortened petioles and leaf blades, and a plump and swelled leaf lamina base .
Most of the genes that define organ polarity existed before the emergence of a flat leaf in seed plants. Based on the analyses of the families of organ polarity genes, such as HD-ZIPIII , ARFs , and ASLs , it is assumed that after the ferns' divergence, multiple paralogs arose in the seed-plant common ancestor . Unlike other polarity genes, YABBYs originated in the lineage leading to seed plants, and it is proposed that they are implicated in the transition of an ancestral shoot-specific network into a leaf-specific one [19, 25]. At least four gene duplication events in the YABBY family led to the emergence of at least five YABBY genes with both novel and redundant functions in the last common ancestor of extant flowering plants [23, 29].
The structural and phylogenetic analysis based on comparison with the YABBY orthologs revealed that each MhyYABBY belonged to one of the four highly conserved clades in the angiosperm YABBY family . In the dendrogram, it is possible to single out a cluster consisting of the FIL/YAB3, YAB5, and YAB2 clades (Figure 3). CRC- and INO-orthologs have formed separate clusters, which corresponds to the previously proposed origin of the CRC and INO genes from different ancestors [20, 23]. The tree composition was not completely congruent with the data of other studies with observations of the two clusters CRC/FIL and YAB5/YAB2/INO , or the two clusters FIL/CRC/ INO and YAB2/YAB5 [20, 29], probably due to the inherent instability of the tree topology depending on the composition of taxa and the mode of analysis.
The obtained tree was consistent with the established phylogenetic relationships among higher plants. The presence of the MhyYABBY paralogs, which are coorthologous to the FIL and INO clades, indicates that the MhyFIL1/ MhyFIL2/MhyFIL3 and MhyINO1/MhyINO2 groups could represent allelic variants (for FIL group), alternative splicing variants (for INO group), or may have originated as a result of a recent gene duplication event unique to the Ericales order.
The simplest explanation of the absence of the YABBY2-like transcripts in all analyzed pinesap transcriptomes may be the low abundance and insufficient transcriptome size. It is also possible that YAB2 homologous genes are expressed at the earlier developmental stages during the formation of lateral organ primordia, which were not analyzed. One of three possible evolutionary scenarios explaining the YABBYs' diversification suggests that YAB2 is the result of the earliest duplication of the YABBY ancestor gene , and therefore, can be associated with the evolution of a leaf stronger than other YABBYs. The YAB2 ortholog is present in photosynthetic Ericaceae relative species Vaccinum corymbosum. Hence, the absence of the YAB2 gene in M. hypopitys may be due to the loss of the gene during the adaptive evolution (deevolution of the genome) of the autotrophic ancestor of M. hypopitys accompanied by the loss of the leaf. MhyYAB2 functions could be partially complemented through neofunctionalization of MhyYAB5 or MhyFIL paralogs. Studies in Oryza sativa and other plants suggest that within the multi-component regulatory network composed of homo- and heterodimers formed by "vegetative" FIL/YAB3, YAB2, and YAB5 orthologs, protein substitutions and replacements are possible [29, 108]. It was shown that in O. sativa the loss of the YAB5 genes was complemented by the FIL and YAB2 paralogs . Accordingly, M. hypopitys MhyYAB2 could have been replaced by the MhyYAB5 or the three MhyFILs. It is also known that YAB2 and YAB5 are important for laminar style and filament morphology evolution in angiosperms, when YAB2 expression over a certain threshold disturbs the balance in the regulatory network, leading to radialization of the laminar structure [26, 27]. Therefore, given the M. hypopitys style and filament radial structure and the expression of the MhyYAB5 gene in flower tissue, the absence of MhyYAB2 transcripts may also indicate a possible substitution of YAB2 by YAB5 in M. hypopitys.
Bioinformatics analysis of the MhyYABBY structural organization revealed the presence of 18 previously known conserved motifs  and 7 putative novel candidate sequences for clade-specific motifs (Figure 5), which may be used as markers to identify appropriate genes. Given the shared evolutionary history of YABBYs, the novel motifs could be biologically relevant and involved in subfamily-specific functions, which need further investigation.
The MhyYABBY gene orthology data are supported by the MhyYABBY expression patterns (Figure 6). In Arabidopsis, "vegetative" genes FIL, YAB3, and YAB5 are expressed in leaves and leaf-like cotyledons, sepals, petals, stamens, and carpels, whereas expression of CRC and INO is restricted to specific floral organs that are evolutionarily derived from leaves . MhyYABBY transcripts were also found in aboveground tissues (bracts and flowers). In the case of the MhyYAB5 gene, its atypical expression in roots indicates that it may have some roles in the development of the M. hypopitys root system. On the other hand, the perennial plant M. hypopitys commonly develops underground adventitious buds on the roots, which presumably contain an embryonic inflorescence , and, thus, the MhyYAB5 gene can be expressed in the buds. Interestingly, in bracts, the expression level of MhyYAB5 is much higher than that of MhyFILs. Given the possible synergy of their action, it can be assumed that the low expression of MhyFILs was compensated by an increase in the expression of the MhyYAB5 gene not only in bracts, but also in roots, or rather in adventitious buds.
The MhyFIL expression profiles indicate the possible subfunctionalization of the paralogs. The MhyFIL3 gene, which according to phylogenetic analysis is at the base of the M. hypopitys FIL clade, is expressed approximately at the same level as MhyFIL1, while the extremely low number of transcripts of the third paralog MhyFIL2 is present only in two of the six transcriptomes (Supplementary Table 3). It is possible that MhyFIL3 and MhyFIL1 may have redundant functions, and MhyFIL2 may be a pseudogene. Given trace amounts of MhyFIL2 transcripts in the root and bud library, similar to MhyYAB5, it is likely that MhyFIL2 may be involved in the development of inflorescence at the early stages after bud dormancy release .
The MhyCRC expression was detected only in flower tissue, confirming its potentially conserved roles in carpel fusion, style/stigma and nectary development, and in the floral meristem termination as it was shown for A. thaliana CRC [33, 38, 102], as well as in vascular development, as indicated by a recent report on the functional role of Pisum sativum cRc . Similar to A. thaliana INO [40, 112], MhyINO1 and MhyINO2 may redundantly define and promote the outer ovule integument growth in M. hypopitys, while MhyFIL1/2/3 and MhyYAB5 may influence the abaxial cell fate in all aboveground lateral organs like their corresponding YAB1/3 and YAB5 orthologs .
It has recently been shown that "vegetative" YABBYs act as transcriptional activators of jasmonate-triggered responses. Jasmonate-induced degradation releases YABBYs from complexes with JAZ3 to mediate anthocyanin accumulation and chlorophyll breakdown . The analysis of M. hypopitys transcriptome data did not reveal JAZ3-like transcripts, which may be consistent with complete chlorophyll loss in M. hypopitys. Thus, such mechanism may become evolutionarily obsolete in mycoheterotrophic plants.
The current study is the first to identify the YABBY genes in a mycoheterotrophic plant devoid of vegetative leaf-like organs. Seven MhyYABBY members were detected and classified in M. hypopitys, and putative protein structure, conserved motifs, and phylogenetic relationship were systematically analyzed. MhyYABBY transcription profiling in different plant tissues indicated the involvement of MhyYABBY proteins in the regulatory network controlling bract and flower formation. Our findings should further the investigation of YABBY functional roles in the regulation of developmental and physiological processes in achlorophyllous plant species and help to reveal possible differences in generally conserved molecular mechanisms underlying plant development and evolution.
Abbreviations CRC: CRABS CLAW FIL: FILAMENTOUS FLOWER HMG: High Mobility Group INO: INNER NO OUTER MhyYABBY: M. hypopitys YABBY YAB: YABBY qRT-PCR: Quantitative RT-PCR.
Conflicts of Interest
The authors declare no conflict of interest.
Konstantin G. Skryabin and Nikolay V. Ravin conceived and designed the research. Elena Z. Kochieva and Mikhail A. Filyushin provided the plant material. Andrey A. Mardanov, Elena Z. Kochieva, Marya A. Slugina, Mikhail A. Filyushin, Olga A. Shulga, and Anna V. Shchennikova performed the transcriptome sequencing and qRT-PCR experiment. Alexey V. Beletsky and Anna V. Shchennikova analyzed the data. Anna V. Shchennikova wrote the paper. All authors read and approved the manuscript.
This work was supported by the Russian Science Foundation (Grant no. 14-24-00175) and was performed using the experimental climate control facility of the Institute of Bioengineering, Research Center of Biotechnology, RAS. The authors would like to thank Dr. Marina Chuenkova for providing language help.
Supplementary 1. Tables 1-7: primers used for qRT-PCR analysis of pinesap gene expression (1); assembled transcripts homologous to the known organ polarity genes in six pinesap transcriptomes (2); MhyYABBY characteristics (3); estimates of evolutionary divergence between nucleotide sequences of the MhyYABBY genes (4); estimates of evolutionary divergence between amino acid sequences of the MhyYABBY proteins (5); set of YABBY orthologues from different plant species used for MEME-mediated identification of clade-specific conserved motifs (6); the significance (p-value) of qRT-PCR results for one gene expression between pinesap tissues (7) (in Microsoft Excel worksheet file (.xls)).
Supplementary 2. Figures 1-2: structural analysis of the MhyYABBY genes and encoded proteins. (1) Alignment of M. hypopitys YABBYgenes (a) and encoded putative proteins (b); (2) novel amino acid motifs predicted in the sequence of the MhyYABBY proteins. (a)-(g) Logos created from aligned sequences and copied directly from MEME graphically represent amino acid conservation: FIL-H (a); FIL-I (b); YAB5-Cm (c); CRC-E (d); INO-B (e); INO-C (f); and INOD (g). The height of the letters in each stack indicates the relative frequency of individual residues at the position (in portable document format (.pdf)).
 A. A. Anderberg, C. Rydin, and M. Kallersjo, "Phylogenetic relationships in the order Ericales s.l.: analyses of molecular data from five genes from the plastid and mitochondrial genomes," American Journal of Botany, vol. 89, no. 4, pp. 677-687, 2002.
 J. R. Leake, "The biology of myco-heterotrophic ("saprophytic") plants," New Phytologist, vol. 127, no. 2, pp. 171-216, 1994.
 E. A. Devers, J. Teply, A. Reinert, N. Gaude, and F. Krajinski, "An endogenous artificial microRNA system for unraveling the function of root endosymbioses related genes in Medicago truncatula," BMC Plant Biology, vol. 13, no. 1, p. 82, 2013.
 G. D. Wallace, "Studies of the Monotropoidiae (Ericaceae): taxonomy and distribution," The Wassman Journal of Biology, vol. 33, pp. 1-88, 1975.
 V. S. F. T. Merckx, J. V. Freudenstein, J. Kissling et al., "Taxonomy and classification," in Mycoheterotrophy, V. S. F. T. Merckx, Ed., pp. 19-101, The Biology of Plants Living on Fungi. Springer Science+Buisness Media, New York, 2013.
 W. N. Stewart and G. W. Rothwell, Paleobotany and the Evolution of Plants, Cambridge University Press, New York, 1993.
 Q. C. B. Cronk, "Plant evolution and development in a post-genomic context," Nature Reviews. Genetics, vol. 2, no. 8, pp. 607-619, 2001.
 J. L. Bowman, Y. Eshed, and S. F. Baum, "Establishment of polarity in angiosperm lateral organs," Trends in Genetics, vol. 18, no. 3, pp. 134-141, 2002.
 D. Beerling and A. Fleming, "Zimmermann's telome theory of megaphyll leaf evolution: a molecular and cellular critique," Current Opinion in Plant Biology, vol. 10, no. 1, pp. 4-12, 2007.
 M. Snow and R. Snow, "The dorsiventrality of leaf primordia," New Phytologist, vol. 58, no. 2, pp. 188-207, 1959.
 J. R. McConnell and M. K. Barton, "Leaf polarity and meristem formation in Arabidopsis," Development, vol. 125, no. 15, pp. 2935-2942, 1998.
 K. Lynn, A. Fernandez, M. Aida et al., "The PINHEAD/ ZWILLE gene acts pleiotropically in Arabidopsis development and has overlapping functions with the ARGONAUTE1 gene," Development, vol. 126, no. 3, pp. 469-481, 1999.
 A. Y. Husbands, D. H. Chitwood, Y. Plavskin, and M. C. P. Timmermans, "Signals and prepatterns: new insights into organ polarity in plants," Genes & Development, vol. 23, no. 17, pp. 1986-1997, 2009.
 T. Yamaguchi, A. Nukazuka, and H. Tsukaya, "Leaf adaxial-abaxial polarity specification and lamina outgrowth: evolution and development," Plant & Cell Physiology, vol. 53, no. 7, pp. 1180-1194, 2012.
 X. Sun, Z. Feng, and L. Meng, "Ectopic expression of the Arabidopsis ASYMMETRIC LEAVES2-LIKE5 (ASL5) gene in cockscomb (Celosia cristata) generates vascular-pattern modifications in lateral organs," Plant Cell, Tissue and Organ Culture (PCTOC), vol. 110, no. 1, pp. 163-169, 2012.
 T. Huang, Y. Harrar, C. Lin et al., "Arabidopsis KANADI1 acts as a transcriptional repressor by interacting with a specific cis-element and regulates auxin biosynthesis, transport, and signaling in opposition to HD-ZIPIII factors," The Plant Cell, vol. 26, no. 1, pp. 246-262, 2014.
 P. Merelo, E. B. Paredes, M. G. Heisler, and S. Wenkel, "The shady side of leaf development: the role of the REVOLUTA/ KANADI1 module in leaf patterning and auxin-mediated growth promotion," Current Opinion in Plant Biology, vol. 35, pp. 111-116, 2017.
 E. M. Meyerowitz, "Genetic control of cell division patterns in developing plants," Cell, vol. 88, no. 3, pp. 299-308, 1997.
 S. K. Floyd and J. L. Bowman, "Gene expression patterns in seed plant shoot meristems and leaves: homoplasy or homology?," Journal of Plant Research, vol. 123, no. 1, pp. 43-55, 2010.
 T. Yamada, S.'y. Yokota, Y. Hirayama, R. Imaichi, M. Kato, and C. S. Gasser, "Ancestral expression patterns and evolutionary diversification of YABBY genes in angiosperms," The Plant Journal, vol. 67, no. 1, pp. 26-36, 2011.
 J. L. Bowman, "The YABBY gene family and abaxial cell fate," Current Opinion in Plant Biology, vol. 3, no. 1, pp. 17-22, 2000.
 E. Kanaya, N. Nakajima, and K. Okada, "Non-sequence-specific DNA binding by the FILAMENTOUS FLOWER protein from Arabidopsis thaliana is reduced by EDTA," The Journal of Biological Chemistry, vol. 277, no. 14, pp. 11957-11964, 2002.
 C. Finet, S. K. Floyd, S. J. Conway, B. Zhong, C. P. Scutt, and J. L. Bowman, "Evolution of the YABBY gene family in seed plants," Evolution & Development, vol. 18, no. 2, pp. 116126, 2016.
 H. L. Liu, Y. Y. Xu, Z. H. Xu, and K. Chong, "A rice YABBY gene, OsYABBY4, preferentially expresses in developing vascular tissue," Development Genes and Evolution, vol. 217, no. 9, pp. 629-637, 2007.
 R. Sarojam, P. G. Sappl, A. Goldshmidt et al., "Differentiating Arabidopsis shoots from leaves by combined YABBY activities," Plant Cell, vol. 22, no. 7, pp. 2113-2130, 2010.
 A. M. R. de Almeida, R. Yockteng, J. Schnable, E. R. Alvarez-Buylla, M. Freeling, and C. D. Specht, "Co-option of the polarity gene network shapes filament morphology in angiosperms," Scientific Reports, vol. 4, no. 1, 2014.
 K. Morioka, R. Yockteng, A. M. R. Almeida, and C. D. Specht, "Loss of YABBY2-like gene expression may underlie the evolution of the laminar style in Canna and contribute to floral morphological diversity in the Zingiberales," Frontiers in Plant Science, vol. 6, 2015.
 D. R. Kelley, D. J. Skinner, and C. S. Gasser, "Roles of polarity determinants in ovule development," The Plant Journal, vol. 57, no. 6, pp. 1054-1064, 2009.
 C. Bartholmes, O. Hidalgo, and S. Gleissberg, "Evolution of the YABBY gene family with emphasis on the basal eudicot Eschscholzia californica (Papaveraceae)," Plant Biology, vol. 14, no. 1, pp. 11-23, 2012.
 A. Z. Worden, J.-H. Lee, T. Mock et al., "Green evolution and dynamic adaptations revealed by genomes of the marine picoeukaryotes Micromonas," Science, vol. 324, no. 5924, pp. 268-272, 2009.
 S. Sawa, K. Watanabe, K. Goto et al., "FILAMENTOUS FLOWER, a meristem and organ identity gene of Arabidopsis, encodes a protein with a zinc finger and HMG-related domains," Genes & Development, vol. 13, no. 9, pp. 1079-1088, 1999.
 J. F. Golz, M. Roccaro, R. Kuzoff, and A. Hudson, "GRAMINIFOLIA promotes growth and polarity of Antirrhinum leaves," Development, vol. 131, no. 15, pp. 36613670, 2004.
 J. L. Bowman and D. R. Smyth, "CRABS CLAW, a gene that regulates carpel and nectary development in Arabidopsis, encodes a novel protein with zinc finger and helix-loophelix domains," Development, vol. 126, no. 11, pp. 2387-2396, 1999.
 Y. Eshed, S. F. Baum, and J. L. Bowman, "Distinct mechanisms promote polarity establishment in carpels of Arabidopsis," Cell, vol. 99, no. 2, pp. 199-209, 1999.
 J. Y. Lee, S. F. Baum, S. H. Oh, C. Z. Jiang, J. C. Chen, and J. L. Bowman, "Recruitment of CRABS CLAW to promote nectary development within the eudicot clade," Development, vol. 132, no. 22, pp. 5021-5032, 2005.
 S. Orashakova, M. Lange, S. Lange, S. Wege, and A. Becker, "The CRABS CLAW ortholog from California poppy (Eschscholzia californica, Papaveraceae), EcCRC, is involved in floral meristem termination, gynoecium differentiation and ovule initiation," The Plant Journal, vol. 58, no. 4, pp. 682-693, 2009.
 H. Li, W. Liang, Y. Hu et al., "Rice MADS6 interacts with the floral homeotic genes SUPERWOMAN1, MADS3, MADS58, MADS13, and DROOPING LEAF in specifying floral organ identities and meristem fate," Plant Cell, vol. 23, no. 7, pp. 2536-2552, 2011.
 W. Sun, W. Huang, Z. Li, H. Lv, H. Huang, and Y. Wang, "Characterization of a Crabs Claw gene in basal eudicot species Epimedium sagittatum (Berberidaceae)," International Journal of Molecular Sciences, vol. 14, no. 1, pp. 1119-1131, 2013.
 S. C. Baker, K. Robinson-Beers, J. M. Villanueva, J. C. Gaiser, and C. S. Gasser, "Interactions among genes regulating ovule development in Arabidopsis thaliana," Genetics, vol. 145, no. 4, pp. 1109-1124, 1997.
 J. M. Villanueva, J. Broadhvest, B. A. Hauser, R. J. Meister, K. Schneitz, and C. S. Gasser, "INNER NO OUTER regulates abaxial-adaxial patterning in Arabidopsis ovules," Genes & Development, vol. 13, no. 23, pp. 3160-3169, 1999.
 M. T. Juarez, R. W. Twigg, and M. C. P. Timmermans, "Specification of adaxial cell fate during maize leaf development," Development, vol. 131, no. 18, pp. 4533-4544, 2004.
 W. Zhao, H. Y. Su, J. Song, X. Y. Zhao, and X. S. Zhang, "Ectopic expression of TaYAB1, a member of YABBY gene family in wheat, causes the partial abaxialization of the adaxial epidermises of leaves and arrests the development of shoot apical meristem in Arabidopsis," Plant Science, vol. 170, no. 2, pp. 364-371, 2006.
 M. Dai, Y. Hu, Y. Zhao, and D. X. Zhou, "Regulatory networks involving YABBY genes in rice shoot development," Plant Signaling & Behavior, vol. 2, no. 5, pp. 399-400, 2007.
 Q. Chen, A. Atkinson, D. Otsuga, T. Christensen, L. Reynolds, and G. N. Drews, "The Arabidopsis FILAMENTOUS FLOWER gene is required for flower formation," Development, vol. 126, no. 12, pp. 2715-2726, 1999.
 B. Cong, L. S. Barrero, and S. D. Tanksley, "Regulatory change in YABBY-like transcription factor led to evolution of extreme fruit size during tomato domestication," Nature Genetics, vol. 40, no. 6, pp. 800-804, 2008.
 H. Nakayama, T. Yamaguchi, and H. Tsukaya, "Expression patterns of AaDL, a CRABS CLAW ortholog in Asparagus asparagoides (Asparagaceae), demonstrate a stepwise evolution of CRC/DL subfamily of YABBY genes," American Journal of Botany, vol. 97, no. 4, pp. 591-600, 2010.
 S. Sawa, T. Ito, Y. Shimura, and K. Okada, "FILAMENTOUS FLOWER controls the formation and development of Arabidopsis inflorescences and floral meristems," Plant Cell, vol. 11, no. 1, pp. 69-86, 1999.
 M. I. Bidartondo and T. D. Bruns, "Extreme specificity in epiparasitic Monotropoideae (Ericaceae): widespread phylogenetic and geographical structure," Molecular Ecology, vol. 10, no. 9, pp. 2285-2295, 2001.
 M. A. Filyushin, N. M. Reshetnikova, E. Z. Kochieva, and K. G. Skryabin, "Polymorphism of sequences and the secondary structure of b/c intron of mitochondrial gene nad1 in Monotropa hypopitys and related Ericaceae species," Biology Bulletin, vol. 43, no. 3, pp. 271-275, 2016.
 M. A. Filyushin, E. Z. Kochieva, and K. G. Skryabin, "Polymorphism of the chloroplast gene rps2 in parasitic plant Monotropa hypopitys L. from the European Russian populations," Russian Journal of Genetics, vol. 53, no. 3, pp. 400405, 2017.
 N. V. Ravin, E. V. Gruzdev, A. V. Beletsky et al., "The loss of photosynthetic pathways in the plastid and nuclear genomes of the non-photosynthetic mycoheterotrophic eudicot Monotropa hypopitys," BMC Plant Biology, vol. 16, Suppl 3, p. 238, 2016.
 G. E. Beatty and J. Provan, "High clonal diversity in threatened peripheral populations of the yellow bird's nest (Hypopitys monotropa; syn. Monotropa hypopitys)," Annals of Botany, vol. 107, no. 4, pp. 663-670, 2011.
 A. V. Beletsky, M. A. Filyushin, E. V. Gruzdev et al., "De novo transcriptome assembly of the mycoheterotrophic plant Monotropa hypopitys," Genomics Data, vol. 11, pp. 60-61, 2017.
 B. Langmead and S. L. Salzberg, "Fast gapped-read alignment with Bowtie 2," Nature Methods, vol. 9, no. 4, pp. 357-359, 2012.
 J. C. Oliveros and Venny, "An interactive tool for comparing lists with Venn's diagrams (2007-2015)," http://bioinfogp .cnb.csic.es/tools/venny/index.html.
 S. Babicki, D. Arndt, A. Marcu et al., "Heatmapper: web-enabled heat mapping for all," Nucleic Acids Research, vol. 44, no. W1, pp. W147-W153, 2016.
 T. L. Bailey and C. Elkan, "Fitting a mixture model by expectation maximization to discover motifs in biopolymers," Proceedings, International Conference on Intelligent Systems for Molecular Biology, vol. 2, pp. 28-36, 1994.
 M. A. Larkin, G. Blackshields, N. P. Brown et al., "Clustal W and Clustal X version 2.0," Bioinformatics, vol. 23, no. 21, pp. 2947-2948, 2007.
 F. Tajima and M. Nei, "Estimation of evolutionary distance between nucleotide sequences," Molecular Biology and Evolution, vol. 1, no. 3, pp. 269-285, 1984.
 K. Tamura and S. Kumar, "Evolutionary distance estimation under heterogeneous substitution pattern among lineages," Molecular Biology and Evolution, vol. 19, no. 10, pp. 1727-1736, 2002.
 K. Tamura, M. Nei, and S. Kumar, "Prospects for inferring very large phylogenies by using the neighbor-joining method," Proceedings of the National Academy of Sciences, vol. 101, no. 30, pp. 11030-11035, 2004.
 S. Kumar, G. Stecher, and K. Tamura, "MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets," Molecular Biology and Evolution, vol. 33, no. 7, pp. 1870-1874, 2016.
 B. Li and C. N. Dewey, "RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome," BMC Bioinformatics, vol. 12, no. 1, p. 323, 2011.
 R. G. Franks, Z. Liu, and R. L. Fischer, "SEUSS and LEUNIG regulate cell proliferation, vascular development and organ polarity in Arabidopsis petals," Planta, vol. 224, no. 4, pp. 801-811, 2006.
 S. Nole-Wilson and B. A. Krizek, "AINTEGUMENTA contributes to organ polarity and regulates growth of lateral organs in combination with YABBY genes," Plant Physiology, vol. 141, no. 3, pp. 977-987, 2006.
 S. M. Doyle, A. Haeger, T. Vain et al., "An early secretory pathway mediated by GNOM-LIKE 1 and GNOM is essential for basal polarity establishment in Arabidopsis thaliana," Proceedings of the National Academy of Sciences of the United States of America, vol. 112, no. 7, pp. E806-E815, 2015.
 H. R. Pires, M. M. Monfared, E. A. Shemyakina, and J. C. Fletcher, "ULTRAPETALA trxG genes interact with KANADI transcription factor genes to regulate Arabidopsis gynoecium patterning," Plant Cell, vol. 26, no. 11, pp. 4345-4361, 2014.
 J. Bou-Torrent, M. Salla-Martret, R. Brandt et al., "ATHB4 and HAT3, two class II HD-ZIP transcription factors, control leaf development in Arabidopsis," Plant Signaling & Behavior, vol. 7, no. 11, pp. 1382-1387, 2012.
 N. Raikhel, "Nuclear targeting in plants," Plant Physiology, vol. 100, no. 4, pp. 1627-1632, 1992.
 R. G. H. Immink, I. A. N. Tonaco, S. de Folter et al., "SEPALLATA3: the "glue" for MADS box transcription factor complex formation," Genome Biology, vol. 10, no. 2, p. R24, 2009.
 S. Wicke, K. F. Muller, C. W. dePamphilis, D. Quandt, S. Bellot, and G. M. Schneeweiss, "Mechanistic model of evolutionary rate variation en route to a nonphotosynthetic lifestyle in plants," Proceedings of the National Academy of Sciences of the United States of America, vol. 113, no. 32, pp. 9045-9050, 2016.
 S. W. Graham, V. K. Y. Lam, and V. S. F. T. Merckx, "Plastomes on the edge: the evolutionary breakdown of mycoheterotroph plastid genomes," New Phytologist, vol. 214, no. 1, pp. 48-55, 2017.
 K. A. Kron, W. S. Judd, P. F. Stevens et al., "Phylogenetic classification of Ericaceae: molecular and morphological evidence," Botanical Review, vol. 68, no. 3, pp. 335-423, 2002.
 T. W. A. Braukmann, M. B. Broe, S. Stefanovic, and J. V. Freudenstein, "On the brink: the highly reduced plastomes of nonphotosynthetic Ericaceae," New Phytologist, vol. 216, no. 1, pp. 254-266, 2017.
 E. V. Gruzdev, A. V. Mardanov, A. V. Beletsky, E. Z. Kochieva, N. V. Ravin, and K. G. Skryabin, "The complete chloroplast genome of parasitic flowering plant Monotropa hypopitys: extensive gene losses and size reduction," Mitochondrial DNA part B Resources, vol. 1, no. 1, pp. 212-213, 2016.
 M. D. Logacheva, M. I. Schelkunov, V. Y. Shtratnikova, M. V. Matveeva, and A. A. Penin, "Comparative analysis of plastid genomes of non-photosynthetic Ericaceae and their photosynthetic relatives," Scientific Reports, vol. 6, no. 1, article 30042, 2016.
 M. I. Bidartondo, "The evolutionary ecology of myco-heterotrophy," New Phytologist, vol. 167, no. 2, pp. 335-352, 2005.
 V. Buchanan-Wollaston, T. Page, E. Harrison et al., "Comparative transcriptome analysis reveals significant differences in gene expression and signalling pathways between developmental and dark/starvation-induced senescence in Arabidopsis," The Plant Journal, vol. 42, no. 4, pp. 567-585, 2005.
 H. Zhang and C. Zhou, "Signal transduction in leaf senescence," Plant Molecular Biology, vol. 82, no. 6, pp. 539-545, 2013.
 H. Tsukaya, "Leaf shape: genetic controls and environmental factors," The International Journal of Developmental Biology, vol. 49, no. 5-6, pp. 547-555, 2005.
 H. Tsukaya, "Mechanism of leaf-shape determination," Annual Review of Plant Biology, vol. 57, no. 1, pp. 477-496, 2006.
 N. A. Eckardt, "YABBY genes and the development and origin of seed plant leaves," Plant Cell, vol. 22, no. 7, p. 2103, 2010.
 S. Mathews and E. M. Kramer, "The evolution of reproductive structures in seed plants: a re-examination based on insights from developmental genetics," The New Phytologist, vol. 194, no. 4, pp. 910-923, 2012.
 K. R. Siegfried, Y. Eshed, S. F. Baum, D. Otsuga, G. N. Drews, and J. L. Bowman, "Members of the YABBY gene family specify abaxial cell fate in Arabidopsis," Development, vol. 126, no. 18, pp. 4117-4128, 1999.
 Y. Wang, H. Huang, Y. Ma, J. Fu, L. Wang, and S. Dai, "Construction and de novo characterization of a transcriptome of Chrysanthemum lavandulifolium: analysis of gene expression patterns in floral bud emergence," Plant Cell, Tissue and Organ Culture (PCTOC), vol. 116, no. 3, pp. 297-309, 2014.
 H. Q. Han, Y. Liu, M. M. Jiang, H. Y. Ge, and H. Y. Chen, "Identification and expression analysis of YABBY family genes associated with fruit shape in tomato (Solanum lycopersicum L.)," Genetics and Molecular Research, vol. 14, no. 2, pp. 7079-7091, 2015.
 J. Kim, J. Yang, R. Yang, R. C. Sicher, C. Chang, and M. L. Tucker, "Transcriptome analysis of soybean leaf abscission identifies transcriptional regulators of organ polarity and cell fate," Frontiers in Plant Science, vol. 7, 2016.
 C. Machida, A. Nakagawa, S. Kojima, H. Takahashi, and Y. Machida, "The complex of ASYMMETRIC LEAVES (AS) proteins plays a central role in antagonistic interactions of genes for leaf polarity specification in Arabidopsis," Wiley Interdisciplinary Reviews: Developmental Biology, vol. 4, no. 6, pp. 655-671, 2015.
 O. A. Shulga, A. V. Shchennikova, A. V. Beletsky et al., "Transcriptome-wide characterization of the MADS-box family in pinesap Monotropa hypopitys reveals flowering conservation in non-photosynthetic myco-heterotrophs," Journal of Plant Growth Regulation, pp. 1-16, 2017.
 J. F. Emery, S. K. Floyd, J. Alvarez et al., "Radial patterning of Arabidopsis shoots by class IIIHD-ZIP and KANADI genes," Current Biology, vol. 13, no. 20, pp. 1768-1774, 2003.
 M. K. Simon, D. J. Skinner, T. L. Gallagher, and C. S. Gasser, "Integument development in Arabidopsis depends on interaction of YABBY protein INNER NO OUTER with co-activators and co-repressors," Genetics, vol. 207, no. 4, pp. 1489-1500, 2017.
 Y. Eshed, A. Izhaki, S. F. Baum, S. K. Floyd, and J. L. Bowman, "Asymmetric leaf development and blade expansion in Arabidopsis are mediated by KANADI and YABBY activities," Development, vol. 131, no. 12, pp. 2997-3006, 2004.
 M. I. Stahle, J. Kuehlich, L. Staron, A. G. von Arnim, and J. F. Golz, "YABBYs and the transcriptional corepressors LEUNIG and LEUNIG_HOMOLOG maintain leaf polarity and meristem activity in Arabidopsis," Plant Cell, vol. 21, no. 10, pp. 3105-3118, 2009.
 O. Bonaccorso, J. E. Lee, L. Puah, C. P. Scutt, and J. F. Golz, "FILAMENTOUS FLOWER controls lateral organ development by acting as both an activator and a repressor," BMC Plant Biology, vol. 12, no. 1, p. 176, 2012.
 C. La Rota, J. Chopard, P. Das et al., "A data-driven integrative model of sepal primordium polarity in Arabidopsis," Plant Cell, vol. 23, no. 12, pp. 4318-4333, 2011.
 I. Pekker, J. P. Alvarez, and Y. Eshed, "Auxin response factors mediate Arabidopsis organ asymmetry via modulation of KANADI activity," Plant Cell, vol. 17, no. 11, pp. 2899-2910, 2005.
 B. J. Reinhart, T. Liu, N. R. Newell et al., "Establishing a framework for the ad/abaxial regulatory network of Arabidopsis: ascertaining targets of class III homeodomain leucine zipper and KANADI regulation," Plant Cell, vol. 25, no. 9, pp. 3228-3249, 2013.
 M. Nakata, N. Matsumoto, R. Tsugeki, E. Rikirsch, T. Laux, and K. Okada, "Roles of the middle domain-specific WUSCHEL-RELATED HOMEOBOX genes in early development of leaves in Arabidopsis," Plant Cell, vol. 24, no. 2, pp. 519-535, 2012.
 H. Iwakawa, M. Iwasaki, S. Kojima et al., "Expression of the ASYMMETRIC LEAVES2 gene in the adaxial domain of Arabidopsis leaves represses cell proliferation in this domain and is critical for the development of properly expanded leaves," The Plant Journal, vol. 51, no. 2, pp. 173-184, 2007.
 B. Xu, Z. Li, Y. Zhu et al., "Arabidopsis genes AS1, AS2, and JAG negatively regulate boundary specifying genes to promote sepal and petal development," Plant Physiology, vol. 146, no. 2, pp. 566-575, 2008.
 M. Iwasaki, H. Takahashi, H. Iwakawa et al., "Dual regulation of ETTIN (ARF3) gene expression by AS1-AS2, which maintains the DNA methylation level, is involved in stabilization of leaf adaxial-abaxial partitioning in Arabidopsis," Development, vol. 140, no. 9, pp. 1958-1969, 2013.
 N. Prunet, P. Morel, A. M. Thierry et al., "REBELOTE, SQUINT, and ULTRAPETALA1 function redundantly in the temporal regulation of floral meristem termination in Arabidopsis thaliana," Plant Cell, vol. 20, no. 4, pp. 901-919, 2008.
 A. V. Shchennikova, O. A. Shulga, E. Z. Kochieva et al., "Homeobox genes encoding WOX transcription factors in the flowering parasitic plant Monotropa hypopitys," Russian Journal of Genetics: Applied Research, vol. 7, no. 7, pp. 781-788, 2017.
 Y. Matsumura, H. Iwakawa, Y. Machida, and C. Machida, "Characterization of genes in the ASYMMETRIC LEAVES2/LATERAL ORGAN BOUNDARIES (AS2/LOB) family in Arabidopsis thaliana, and functional and molecular comparisons between AS2 and other family members," The Plant Journal, vol. 58, no. 3, pp. 525-537, 2009.
 E. Semiarti, Y. Ueno, H. Tsukaya, H. Iwakawa, C. Machida, and Y. Machida, "The ASYMMETRICLEAVES2 gene of Arabidopsis thaliana regulates formation of a symmetric lamina, establishment of venation and repression of meristem-related homeobox genes in leaves," Development, vol. 128, no. 10, pp. 1771-1783, 2001.
 M. J. Prigge and S. E. Clark, "Evolution of the class III HD-Zip gene family in land plants," Evolution & Development, vol. 8, no. 4, pp. 350-361, 2006.
 A. S. Chanderbali, F. He, P. S. Soltis, and D. E. Soltis, "Out of the water: origin and diversification of the LBD gene family," Molecular Biology and Evolution, vol. 32, no. 8, pp. 1996-2000, 2015.
 T. Toriba, K. Harada, A. Takamura et al., "Molecular characterization the YABBY gene family in Oryza sativa and expression analysis of OsYABBY1," Molecular Genetics and Genomics, vol. 277, no. 5, pp. 457-468, 2007.
 F. Vandelook and J. A. Van Assche, "Temperature conditions control embryo growth and seed germination of Corydalis solida (L.) Clairv., a temperate forest spring geophyte," Plant Biology (Stuttgart, Germany), vol. 11, no. 6, pp. 899-906, 2009.
 J. R. Leake, S. L. McKendrick, M. Bidartondo, and D. J. Read, "Symbiotic germination and development of the myco-heterotroph Monotropa hypopitys in nature and its requirement for locally distributed Tricholoma spp," New Phytologist, vol. 163, no. 2, pp. 405-423, 2004.
 C. Fourquin, A. Primo, I. Martinez-Fernandez, E. Huet-Trujillo, and C. Ferrandiz, "The CRC orthologue from Pisum sativum shows conserved functions in carpel morphogenesis and vascular development," Annals of Botany, vol. 114, no. 7, pp. 1535-1544, 2014.
 M. K. Simon, L. A. Williams, K. Brady-Passerini, R. H. Brown, and C. S. Gasser, "Positive- and negative-acting regulatory elements contribute to the tissue-specific expression of INNER NO OUTER, a YABBY-type transcription factor gene in Arabidopsis," BMC Plant Biology, vol. 12, no. 1, p. 214, 2012.
 M. Boter, J. F. Golz, S. Gimenez-Ibanez, G. Fernandez-Barbero, J. M. Franco-Zorrilla, and R. Solano, "FILAMENTOUS FLOWER is a direct target of JAZ3 and modulates responses to Jasmonate," Plant Cell, vol. 27, no. 11, pp. 3160-3174, 2015.
Anna V. Shchennikova (iD), (1) Marya A. Slugina, (1,2) Alexey V. Beletsky, (1) Mikhail A. Filyushin, (1) Andrey A. Mardanov, (1) Olga A. Shulga, (1) Elena Z. Kochieva, (1,2) Nikolay V. Ravin (iD), (1) and Konstantin G. Skryabin (1,2)
(1) Federal State Institution "Federal Research Centre "Fundamentals of Biotechnology" of the Russian Academy of Sciences", Moscow 119071, Russia
(2) Lomonosov Moscow State University, Moscow 119991, Russia
Correspondence should be addressed to Anna V. Shchennikova; firstname.lastname@example.org
Received 16 January 2018; Revised 2 March 2018; Accepted 18 March 2018; Published 24 April 2018
Academic Editor: Ferenc Olasz
Caption: Figure 1: Venn diagrams of the overlap between the whole set of reads in paired transcriptomes (a) and in three transcriptomes for each of the two plants analyzed (b).
Caption: Figure 2: Heat maps of the organ polarity gene expression in six transcriptomes (a) and the mean expression of these genes in paired transcriptomes (b).
Caption: Figure 3: Alignment of conserved zinc-finger and HMG-like domains from putative M. hypopitys YABBY proteins 1-7, known as A. thaliana YABBYs (AthINO, AF195047; AthCRC, AF132606; AthFIL, AF136538; AthYAB3, AF136540; AthYAB5, NP_850081; AthYAB2, AF136539), and Pinus taeda PtYAB (DR100835; gymnosperms).
Caption: Figure 4: A phylogenetic tree of MhyYABBY proteins was generated with 83 full-size YABBY-like amino acid sequences from M. hypopitys, A. thaliana, and other angiosperm species. The M. commode (Chlorophyta) YABBY-like protein was used as outgroup. The Amborella trichopoda, Ericales, Rosids, Lamiids (euAsterids), and Campanulids (euAsterids) sequences are colored in dark blue, blue, burgundy, green, and red, respectively. The numbers next to the nodes represent bootstrap values from 1000 replicates.
Caption: Figure 5: Amino acid clade-specific motifs predicted in M. hypopitys YABBY proteins. MhyFIL1/2/3 (a); MhyYAB5 (b); MhyCRC (c); MhyINO1/2 (d). Previously suggested clade-specific motifs  are shown as boxes. MEME-predicted novel motifs are marked by arrowheads and shown as letter sequences.
Caption: Figure 6: Expression profiles of M. hypopitys YABBYgenes in the bract, flower, root and bud, and haustoria-enriched root tissues. Expression was estimated as the transcript number per million equal to the sum of total transcripts in the tissue. The M. hypopitys adult plant; br--bract, fl--flower (a); transcriptome-based expression pattern of the MhyFIL1, MhyFIL2, MhyFIL3, MhyYAB5, MhyCRC, MhyINO1, and MhyINO2 genes in pinesap tissues (b); the M. hypopitys roots with adventitious buds; ab--adventitious bud, rs--growing reproductive stem; scale bar = 1 cm (c); relative expression (qRT-PCR) of the pattern of the MhyCRC, MhyINO1, MhyINO2, MhyFIL3, and MhyYAB5 genes in pinesap tissues, F--flower, B--bract, R--root (d).
Caption: Figure 7: A hypothetical scenario for the functioning of the MhyYABBYs during the development of M. hypopitys bracts and flowers. (a) M. hypopitys raceme. (b) M. hypopitys sterile bract. (c) M. hypopitys flower with a partially removed perianth. (d) Scheme of possible relations between "abaxial" and "adaxial" factors in the M. hypopitys bracts. (e) Scheme of possible interactions of MhyYABBYs in the M. hypopitys flowers. Br--bract, StBr--sterile bract, Ad--adaxial side, Ab--abaxial side, LM--leaf margin, Se--sepal, Pe--petal, St--stamen, Pi--pistil, Int--integument. Scale bar = 0.5 cm.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Shchennikova, Anna V.; Slugina, Marya A.; Beletsky, Alexey V.; Filyushin, Mikhail A.; Mardanov, Andr|
|Publication:||International Journal of Genomics|
|Date:||Jan 1, 2018|
|Previous Article:||Expanding the miRNA Transcriptome of Human Kidney and Renal Cell Carcinoma.|
|Next Article:||Gene Coexpression Network Comparison via Persistent Homology.|