Deep sequencing and microarray hybridization identify conserved and species-specific microRNAs during somatic embryogenesis in hybrid yellow poplar

Epigenetic modifications and miRNAs determine the transition of somatic cells into somatic embryos

KEY MESSAGE

This review discusses the epigenetic changes during somatic embryo (SE) development, highlights the genes and miRNAs involved in the transition of somatic cells into SEs as a result of epigenetic changes, and draws insights on biotechnological opportunities to study SE development. Somatic embryogenesis from somatic cells occurs in a series of steps. The transition of somatic cells into somatic embryos (SEs) is the most critical step under genetic and epigenetic regulations. Major regulatory genes such as SERK, WUS, BBM, FUS3/FUSA3, AGL15, and PKL, control SE steps and development by turning on and off other regulatory genes. Gene transcription profiles of somatic cells during SE development is the result of epigenetic changes, such as DNA and histone protein modifications, that control and decide the fate of SE formation. Depending on the type of somatic cells and the treatment with plant growth regulators, epigenetic changes take place dynamically. Either hypermethylation or hypomethylation of SE-related genes promotes the transition of somatic cells. For example, the reduced levels of DNA methylation of SERK and WUS promotes SE initiation. Histone modifications also promote SE induction by regulating SE-related genes in somatic cells. In addition, miRNAs contribute to the various stages of SE by regulating the expression of auxin signaling pathway genes (TIR1, AFB2, ARF6, and ARF8), transcription factors (CUC1 and CUC2), and growth-regulating factors (GRFs) involved in SE formation. These epigenetic and miRNA functions are unique and have the potential to regenerate bipolar structures from somatic cells when a pluripotent state is induced. However, an integrated overview of the key regulators involved in SE development and downstream processes is lacking. Therefore, this review discusses epigenetic modifications involved in SE development, SE-related genes and miRNAs associated with epigenetics, and common cis-regulatory elements in the promoters of SE-related genes. Finally, we highlight future biotechnological opportunities to alter epigenetic pathways using the genome editing tool and to study the transition mechanism of somatic cells.

Genome-wide analysis of the GRAS gene family in Liriodendron chinense reveals the putative function in abiotic stress and plant development

Introduction

GRAS genes encode plant-specific transcription factors that play essential roles in plant growth and development. However, the members and the function of the GRAS gene family have not been reported in Liriodendron chinense. L. chinense, a tree species in the Magnolia family that produces excellent timber for daily life and industry. In addition, it is a good relict species for plant evolution research.

Methods

Therefore, we conducted a genome-wide study of the LcGRAS gene family and identified 49 LcGRAS genes in L. chinense.

Results

We found that LcGRAS could be divided into 13 sub-groups, among which there is a unique branch named HAM-t. We carried out RNA sequencing analysis of the somatic embryos from L. chinense and found that LcGRAS genes are mainly expressed after heart-stage embryo development, suggesting that LcGRAS may have a function during somatic embryogenesis. We also investigated whether GRAS genes are responsive to stress by carrying out RNA sequencing (RNA-seq) analysis, and we found that the genes in the PAT subfamily were activated upon stress treatment, suggesting that these genes may help plants survive stressful environments. We found that PIF was downregulated and COR was upregulated after the transient overexpression of PATs, suggesting that PAT may be upstream regulators of cold stress.

Discussion

Collectively, LcGRAS genes are conserved and play essential roles in plant development and adaptation to abiotic stress.

Genome-Wide Identification and Expression Analysis of CCO Gene Family in Liriodendron chinense

Carotenoid cleavage oxygenase (CCO) is an enzyme that can catalyze carotenoids to volatile aromatic substances and participate in the biosynthesis of two important phytohormones, i.e., abscisic acid (ABA) and strigolactone (SL). However, the genome-wide identification and analysis of the CCO gene family in the rare and endangered woody plant Liriodendron chinense has not been reported. Here, we performed a genome-wide analysis of the CCO gene family in the L. chinense genome and examined its expression pattern during different developmental processes and in response to various abiotic stresses. A total of 10 LcCCO genes were identified and divided into 6 subfamilies according to the phylogenetic analysis. Subcellular localization prediction showed that most of the LcCCO proteins were located in the cytoplasm. Gene replication analysis showed that segmental and tandem duplication contributed to the expansion of this gene family in the L. chinense genome. Cis-element prediction showed that cis-elements related to plant hormones, stress and light response were widely distributed in the promoter regions of LcCCO genes. Gene expression profile analysis showed that LcNCED3b was extensively involved in somatic embryogenesis, especially the somatic embryo maturation, as well as in response to heat and cold stress in leaves. Furthermore, qRT-PCR analysis showed that LcNCED3b obviously responded to drought stress in roots and leaves. This study provides a comprehensive overview of the LcCCO gene family and a potential gene target for the optimization of the somatic embryogenesis system and resistance breeding in the valuable forest tree L. chinense.

Uncovering a Phenomenon of Active Hormone Transcriptional Regulation during Early Somatic Embryogenesis in Medicago sativa

Somatic embryogenesis (SE) is a developmental process in which somatic cells undergo dedifferentiation to become plant stem cells, and redifferentiation to become a whole embryo. SE is a prerequisite for molecular breeding and is an excellent platform to study cell development in the majority of plant species. However, the molecular mechanism involved in M. sativa somatic embryonic induction, embryonic and maturation is unclear. This study was designed to examine the differentially expressed genes (DEGs) and miRNA roles during somatic embryonic induction, embryonic and maturation. The cut cotyledon (ICE), non-embryogenic callus (NEC), embryogenic callus (EC) and cotyledon embryo (CE) were selected for transcriptome and small RNA sequencing. The results showed that 17,251 DEGs, and 177 known and 110 novel miRNAs families were involved in embryonic induction (ICE to NEC), embryonic (NEC to EC), and maturation (EC to CE). Expression patterns and functional classification analysis showed several novel genes and miRNAs involved in SE. Moreover, embryonic induction is an active process of molecular regulation, and hormonal signal transduction related to pathways involved in the whole SE. Finally, a miRNA–target interaction network was proposed during M. sativa SE. This study provides novel perspectives to comprehend the molecular mechanisms in M. sativa SE.

miR171 modulates induction of somatic embryogenesis in citrus callus

Overexpression of miR171 restored SE competence in the recalcitrant citrus callus, and inhibition of miR171 function weakened SE competence in the strongly embryogenic citrus callus. Somatic embryogenesis (SE) is an important way of in vitro regeneration for plants. For perennial woody crops such as citrus, embryogenic callus is usually induced from unfertilized aborted ovules and widely used in biotechnology aided breeding. However, SE capacity always declines in callus during subculture, which makes regeneration difficult and hinders the application of biotechnology. We previously found that miR171 may be a regulator of SE in citrus, based on the abundant expression of csi-miR171c in the embryogenic callus and during SE of citrus. Here, we report that miR171 promotes SE and is required for SE in citrus. Overexpression of miR171 restored SE competence in the recalcitrant callus of 'Guoqing No.1' Satsuma mandarin (G1), whereas inhibition of miR171 function by Short Tandem Target Mimic (STTM) weakened SE competence in the strongly embryogenic callus of 'Valencia' sweet orange (V). The comparative transcriptomic analysis in miR171 overexpressed callus line (OE) and the wild type callus (WT) indicated that overexpression of miR171 decreased the expression level of its SCARECROW-LIKE (CsSCL) targets, and activated stress response related biological processes and metabolic processes that are required for cell differentiation. However, CsSCLs were up-regulated in the OE callus during SE induction process, which activated the cell division and developmental processes that are required for embryogenesis progress. Our results validate the function of miR171 in regulation of SE and reveal the biological responses provoked by miR171 in citrus that may promote SE.

Agrobacterium-Mediated Genetic Transformation of Embryogenic Callus in a Liriodendron Hybrid (L. Chinense × L. Tulipifera)

A highly efficient genetic transformation system of Liriodendron hybrid embryogenic calli through Agrobacterium-mediated genetic transformation was established and optimized. The Agrobacterium tumefaciens strain EHA105, harboring the plasmid pBI121, which contained the ß-glucuronidase (GUS) gene and neomycin phosphotransferase II (npt II) gene under the control of the CaMV35S promoter, was used for transformation. Embryogenic calli were used as the starting explant to study several factors affecting the Agrobacterium-mediated genetic transformation of the Liriodendron hybrid, including the effects of various media, selection by different Geneticin (G418) concentrations, pre-culture period, Agrobacterium optical density, infection duration, co-cultivation period, and delayed selection. Transformed embryogenic calli were obtained through selection on medium containing 90 mg L-1 G418. Plant regeneration was achieved and selected via somatic embryogenesis on medium containing 15 mg L-1 G418. The optimal conditions included a pre-culture time of 2 days, a co-culture time of 3 days, an optimal infection time of 10 min, and a delayed selection time of 7 days. These conditions, combined with an OD600 value of 0.6, remarkably enhanced the transformation rate. The results of GUS chemical tissue staining, polymerase chain reaction (PCR), and southern blot analysis demonstrated that the GUS gene was successfully expressed and integrated into the Liriodendron hybrid genome. A transformation efficiency of 60.7% was achieved for the regenerated callus clumps. Transgenic plantlets were obtained in 5 months, and the PCR analysis showed that 97.5% of plants from the tested G418-resistant lines were PCR positive. The study of the Liriodendron hybrid reported here will facilitate the insertion of functional genes into the Liriodendron hybrid via Agrobacterium-mediated transformation.

miR172 Regulates WUS during Somatic Embryogenesis in Arabidopsis via AP2

In plants, the embryogenic transition of somatic cells requires the reprogramming of the cell transcriptome, which is under the control of genetic and epigenetic factors. Correspondingly, the extensive modulation of genes encoding transcription factors and miRNAs has been indicated as controlling the induction of somatic embryogenesis in Arabidopsis and other plants. Among the MIRNAs that have a differential expression during somatic embryogenesis, members of the MIRNA172 gene family have been identified, which implies a role of miR172 in controlling the embryogenic transition in Arabidopsis. In the present study, we found a disturbed expression of both MIRNA172 and candidate miR172-target genes, including AP2, TOE1, TOE2, TOE3, SMZ and SNZ, that negatively affected the embryogenic response of transgenic explants. Next, we examined the role of AP2 in the miR172-mediated mechanism that controls the embryogenic response. We found some evidence that by controlling AP2, miR172 might repress the WUS that has an important function in embryogenic induction. We showed that the mechanism of the miR172-AP2-controlled repression of WUS involves histone acetylation. We observed the upregulation of the WUS transcripts in an embryogenic culture that was overexpressing AP2 and treated with trichostatin A (TSA), which is an inhibitor of HDAC histone deacetylases. The increased expression of the WUS gene in the embryogenic culture of the hdac mutants further confirmed the role of histone acetylation in WUS control during somatic embryogenesis. A chromatin-immunoprecipitation analysis provided evidence about the contribution of HDA6/19-mediated histone deacetylation to AP2-controlled WUS repression during embryogenic induction. The upstream regulatory elements of the miR172-AP2-WUS pathway might involve the miR156-controlled SPL9/SPL10, which control the level of mature miR172 in an embryogenic culture.

Transcriptome and proteome analysis suggest enhanced photosynthesis in tetraploid Liriodendron sino-americanum

Polyploidy generally provides an advantage in phenotypic variation and growth vigor. However, the underlying mechanisms remain poorly understood. The tetraploid Liriodendron sino-americanum (Liriodendron × sinoamericanum P.C Yieh ex C.B. Shang & Zhang R.Wang) exhibits altered morphology compared with its diploid counterpart, including larger, thicker and deeper green leaves, bigger stomata, thicker stems and increased tree height. Such characteristics can be useful in ornamental and industrial applications. To elucidate the molecular mechanisms behind this variation, we performed a comparative transcriptome and proteome analysis. Our transcriptome data indicated that some photosynthesis genes and pathways were differentially altered and enriched in tetraploid L. sino-americanum, mainly related to F-type ATPase, the cytochrome b6/f complex, photosynthetic electron transport, the light harvesting chlorophyll protein complexes, and photosystem I and II. Most of the differentially expressed proteins we could identify are also involved in photosynthesis. Our physiological results showed that tetraploids have an enhanced photosynthetic capacity, concomitant with great levels of sugar and starch in leaves. This suggests that tetraploid L. sino-americanum might experience comprehensive transcriptome reprogramming of genes related to photosynthesis. This study has especially emphasized molecular changes involved in photosynthesis that accompany polyploidy, and provides a possible explanation for the altered phenotype of polyploidy plants in comparison with their diploid form.

Selection of reference genes for gene expression analysis in Liriodendron hybrids' somatic embryogenesis and germinative tissues

The differential expression of genes is crucial for plant somatic embryogenesis (SE), and the accurate quantification of gene expression levels relies on choosing appropriate reference genes. To select the most suitable reference genes for SE studies, 10 commonly used reference genes were examined in synchronized somatic embryogenic and subsequent germinative cultures of Liriodendron hybrids by using quantitative real-time reverse transcription PCR. Four popular normalization algorithms: geNorm, NormFinder, Bestkeeper and Delta-Ct were used to select and validate the suitable reference genes. The results showed that elongation factor 1-gamma, histone H1 linker protein, glyceraldehyde-3-phosphate dehydrogenase and α-tubulin were suitable for SE tissues, while elongation factor 1-gamma and actin were best for the germinative organ tissues. Our work will benefit future studies of gene expression and functional analyses of SE in Liriodendron hybrids. It is also serves as a guide of reference gene selection in early embryonic gene expression analyses for other woody plant species.

miRNA expression profiling and zeatin dynamic changes in a new model system of in vivo indirect regeneration of tomato

Callus formation and adventitious shoot differentiation could be observed on the cut surface of completely decapitated tomato plants. We propose that this process can be used as a model system to investigate the mechanisms that regulate indirect regeneration of higher plants without the addition of exogenous hormones. This study analyzed the patterns of trans-zeatin and miRNA expression during in vivo regeneration of tomato. Analysis of trans-zeatin revealed that the hormone cytokinin played an important role in in vivo regeneration of tomato. Among 183 miRNAs and 1168 predicted target genes sequences identified, 93 miRNAs and 505 potential targets were selected based on differential expression levels for further characterization. Expression patterns of six miRNAs, including sly-miR166, sly-miR167, sly-miR396, sly-miR397, novel 156, and novel 128, were further validated by qRT-PCR. We speculate that sly-miR156, sly-miR160, sly-miR166, and sly-miR397 play major roles in callus formation of tomato during in vivo regeneration by regulating cytokinin, IAA, and laccase levels. Overall, our microRNA sequence and target analyses of callus formation during in vivo regeneration of tomato provide novel insights into the regulation of regeneration in higher plants.

miRNAs as key regulators via targeting the phytohormone signaling pathways during somatic embryogenesis of plants

Somatic embryogenesis is the regeneration of embryos from the somatic cell via dedifferentiation and redifferentiation without the occurrence of fertilization. A complex network of genes regulates the somatic embryogenesis process. Especially, microRNAs (miRNAs) have emerged as key regulators by affecting phytohormone biosynthesis, transport and signal transduction pathways. miRNAs are small, non-coding small RNA regulatory molecules involved in various developmental processes including somatic embryogenesis. Several types of miRNAs such as miR156, miR157, miR 159, miR 160, miR165, miR166, miR167, miR390, miR393 and miR396 have been reported to intricate in regulating somatic embryogenesis via targeting the phytohormone signaling pathways. Here we review current research progress on the miRNA-mediated regulation involved in somatic embryogenesis via regulating auxin, ethylene, abscisic acid and cytokinin signaling pathways. Further, we also discussed the possible role of other phytohormone signaling pathways such as gibberellins, jasmonates, nitric oxide, polyamines and brassinosteroids. Finally, we conclude by discussing the expression of miRNAs and their targets involved in somatic embryogenesis and possible regulatory mechanisms cross talk with phytohormones during somatic embryogenesis.

AGL15 Controls the Embryogenic Reprogramming of Somatic Cells in Arabidopsis through the Histone Acetylation-Mediated Repression of the miRNA Biogenesis Genes

The embryogenic transition of somatic cells requires an extensive reprogramming of the cell transcriptome. Relevantly, the extensive modulation of the genes that have a regulatory function, in particular the genes encoding the transcription factors (TFs) and miRNAs, have been indicated as controlling somatic embryogenesis (SE) that is induced in vitro in the somatic cells of plants. Identifying the regulatory relationships between the TFs and miRNAs during SE induction is of central importance for understanding the complex regulatory interplay that fine-tunes a cell transcriptome during the embryogenic transition. Hence, here, we analysed the regulatory relationships between AGL15 (AGAMOUS-LIKE 15) TF and miR156 in an embryogenic culture of Arabidopsis. Both AGL15 and miR156 control SE induction and AGL15 has been reported to target the MIR156 genes in planta. The results showed that AGL15 contributes to the regulation of miR156 in an embryogenic culture at two levels that involve the activation of the MIR156 transcription and the containment of the abundance of mature miR156 by repressing the miRNA biogenesis genes DCL1 (DICER-LIKE1), SERRATE and HEN1 (HUA-ENHANCER1). To repress the miRNA biogenesis genes AGL15 seems to co-operate with the TOPLESS co-repressors (TPL and TPR1-4), which are components of the SIN3/HDAC silencing complex. The impact of TSA (trichostatin A), an inhibitor of the HDAC histone deacetylases, on the expression of the miRNA biogenesis genes together with the ChIP results implies that histone deacetylation is involved in the AGL15-mediated repression of miRNA processing. The results indicate that HDAC6 and HDAC19 histone deacetylases might co-operate with AGL15 in silencing the complex that controls the abundance of miR156 during embryogenic induction. This study provides new evidence about the histone acetylation-mediated control of the miRNA pathways during the embryogenic reprogramming of plant somatic cells and the essential role of AGL15 in this regulatory mechanism.

Epigenetic Regulation of Auxin-Induced Somatic Embryogenesis in Plants

Somatic embryogenesis (SE) that is induced in plant explants in response to auxin treatment is closely associated with an extensive genetic reprogramming of the cell transcriptome. The significant modulation of the gene transcription profiles during SE induction results from the epigenetic factors that fine-tune the gene expression towards embryogenic development. Among these factors, microRNA molecules (miRNAs) contribute to the post-transcriptional regulation of gene expression. In the past few years, several miRNAs that regulate the SE-involved transcription factors (TFs) have been identified, and most of them were involved in the auxin-related processes, including auxin metabolism and signaling. In addition to miRNAs, chemical modifications of DNA and chromatin, in particular the methylation of DNA and histones and histone acetylation, have been shown to shape the SE transcriptomes. In response to auxin, these epigenetic modifications regulate the chromatin structure, and hence essentially contribute to the control of gene expression during SE induction. In this paper, we describe the current state of knowledge with regard to the SE epigenome. The complex interactions within and between the epigenetic factors, the key SE TFs that have been revealed, and the relationships between the SE epigenome and auxin-related processes such as auxin perception, metabolism, and signaling are highlighted.

Genome-Wide microRNA Profiling Using Oligonucleotide Microarray Reveals Regulatory Networks of microRNAs in Nicotiana benthamiana During Beet Necrotic Yellow Vein Virus Infection

Beet necrotic yellow vein virus (BNYVV) infections induce stunting and leaf curling, as well as root and floral developmental defects and leaf senescence in Nicotiana benthamiana. A microarray analysis with probes capable of detecting 1596 candidate microRNAs (miRNAs) was conducted to investigate differentially expressed miRNAs and their targets upon BNYVV infection of N. benthamiana plants. Eight species-specific miRNAs of N. benthamiana were identified. Comprehensive characterization of the N. benthamiana microRNA profile in response to the BNYVV infection revealed that 129 miRNAs were altered, including four species-specific miRNAs. The targets of the differentially expressed miRNAs were predicted accordingly. The expressions of miR164, 160, and 393 were up-regulated by BNYVV infection, and those of their target genes, NAC21/22, ARF17/18, and TIR, were down-regulated. GRF1, which is a target of miR396, was also down-regulated. Further genetic analysis of GRF1, by Tobacco rattle virus-induced gene silencing, assay confirmed the involvement of GRF1 in the symptom development during BNYVV infection. BNYVV infection also induced the up-regulation of miR168 and miR398. The miR398 was predicted to target umecyanin, and silencing of umecyanin could enhance plant resistance against viruses, suggesting the activation of primary defense response to BNYVV infection in N. benthamiana. These results provide a global profile of miRNA changes induced by BNYVV infection and enhance our understanding of the mechanisms underlying BNYVV pathogenesis.

The miR396-GRF Regulatory Module Controls the Embryogenic Response in Arabidopsis via an Auxin-Related Pathway

In plants, microRNAs have been indicated to control various developmental processes, including somatic embryogenesis (SE), which is triggered in the in vitro cultured somatic cells of plants. Although a transcriptomic analysis has indicated that numerous MIRNAs are differentially expressed in the SE of different plants, the role of specific miRNAs in the embryogenic reprogramming of the somatic cell transcriptome is still poorly understood. In this study, we focused on performing a functional analysis of miR396 in SE given that the transcripts of MIR396 genes and the mature molecules of miR396 were found to be increased during an SE culture of Arabidopsis. In terms of miR396 in embryogenic induction, we observed the SE-associated expression pattern of MIR396b in explants of the β-glucuronidase (GUS) reporter line. In order to gain insight into the miR396-controlled mechanism that is involved in SE induction, the embryogenic response of mir396 mutants and the 35S:MIR396b overexpressor line to media with different 2,4-Dichlorophenoxyacetic acid (2,4-D) concentrations was evaluated. The results suggested that miR396 might contribute to SE induction by controlling the sensitivity of tissues to auxin treatment. Within the targets of miR396 that are associated with SE induction, we identified genes encoding the GROWTH-REGULATING FACTOR (GRF) transcription factors, including GRF1, GRF4, GRF7, GRF8, and GRF9. Moreover, the study suggested a regulatory relationship between miR396, GRF, and the PLETHORA (PLT1 and PLT2) genes during SE induction. A complex regulatory relationship within the miR396–GRF1/4/8/9–PLT1/2 module that involves the negative and positive control of GRFs and PLT (respectively) by miR396 might be assumed.

Exosomal miRNAs in central nervous system diseases: biomarkers, pathological mediators, protective factors and therapeutic agents

Exosomes are small bilipid layer-enclosed extracellular vesicles that can be found in tissues and biological fluids. As a key cell-to-cell and distant communication mediator, exosomes are involved in various central nervous system (CNS) diseases, potentially through transferring their contents such as proteins, lipids and nucleic acids to the target cells. Exosomal miRNAs, which are small non-coding RNAs in the exosomes, are known to be more stable than free miRNAs and therefore have lasting effects on disease-related gene expressions. There are distinct profiles of exosomal miRNAs in different types of CNS diseases even before the onset of irreversible neurological damages, indicating that exosomal miRNAs within tissues and biological fluids could serve as promising biomarkers. Emerging evidence has also demonstrated the pathological effects of several exosomal miRNAs in CNS diseases via specific modulation of disease-related factors. Moreover, exosomes carry therapeutically beneficial miRNAs across the blood-brain-barrier, which can be exploited as a powerful drug delivery tool to help alleviating multiple CNS diseases. In this review, we summarize the recent progress made in understanding the biological roles of exosomal miRNAs as potential diagnostic biomarkers, pathological regulators, and therapeutic targets/drugs for CNS diseases. A comprehensive discussion of the main concerns and challenges for the applications of exosomal miRNAs in the clinical setting is also provided.

New clues into the mechanisms of rice domestication

Domestication of rice involved incorporation of specific yield-related changes in wild species of rice. This agricultural process has been of significant interest for plant biologists. The recent advance in genomics has provided new tools to investigate the genetic basis and consequences of domestication. Several genes involved in domestication and diversification process have been characterized, and as expected, this list is over-represented by transcription factors and their cofactors. Most often the modification orchestrated expression levels of genes such as those coding for transcription factors. It has been proposed that transcriptional regulators and their regulation is likely a major theme controlling morphological differences between crops and their progenitors. However, recent data indicate that single amino acid changes in genes coding for key proteins as well as epigenetic and small RNA-mediated pathways also contributed towards domesticationassociated phenotypes.

The role of miRNA in somatic embryogenesis

Somatic embryogenesis (SEG) is one of the best techniques for mass production of economically important plants. It is also used for the study of morphology, anatomy, physiology, genetics and molecular mechanism of embryo development. Somatic Embryos (SE) are bipolar structures that develop from a cell other than a gamete or zygote. SEG reflects the unique developmental potential of plant somatic cells, resulting in the transition of the differentiated somatic cells to embryogenic cells to follow the zygotic embryo stages. There are several biochemical and physiological processes that transformed a single somatic cell to a whole plant. SE studies provide insight into cell mechanisms governing the totipotency process in plants. Previously, in vitro studies have suggested the role of various regulatory genes in embryogenic transition that are triggered by plant hormones in response to stress. The omic studies identify the specific genes, transcripts, and proteins required for somatic embryogenesis development. MicroRNAs (miRNAs) are small, 19-24 nucleotides (nt), non-coding small RNA regulatory molecules controlling a large number of biological processes. In addition to their role in SEG, miRNAs play vital role in plant development, secondary metabolite synthesis and metabolism of macromolecules, hormone signal transduction, and tolerance of plants to biotic and abiotic stresses. During last decade several types of miRNAs involved in SEG have been reported. Among these miRNAs, miR156, miR162, miR166a, miR167, miR168, miR171a/b, miR171c, miR393, miR397 and miR398 played very active role during various stages of SEG. In this review, we highlighted the role of these as well as other miRNAs in some economically important plants.

Major Domestication-Related Phenotypes in Indica Rice Are Due to Loss of miRNA-Mediated Laccase Silencing

Domestication of rice (Oryza sativa) included conversion of perennial wild species with few seeds to short plants that produced abundant seeds. Most domestication-associated changes were due to variations in transcription factors and other key proteins such as enzymes. Here, we show that multiple yield-related traits associated with indica rice domestication are linked to micro (mi) RNA-mediated regulation. Analysis of small (s) RNA data sets from cultivated indica rice lines, a few landraces, and two wild relatives of rice revealed the presence of abundant 22-nucleotide (nt) reads in wild relatives that mapped to miR397 precursors. miR397 was expressed at very high levels in wild relatives and at negligible levels in high-yielding cultivated lines. In its genera-specific form of 22-nt, miR397 targeted mRNAs encoding laccases that decayed and induced robust secondary cascade silencing in wild species that required RNA-dependent RNA polymerase 6. In wild species of rice, reduced expression of laccases resulted in low lignification. As expected, overexpression of miR397 induced de-domestication phenotypes. At least 26 uncharacterized QTLs previously implicated in rice yield overlapped with laccases and miR397 genes. These results suggest that miRNAs contribute to rice domestication-associated phenotypes.

Whole Mount in situ Localization of miRNAs and mRNAs During Somatic Embryogenesis in Arabidopsis

Somatic embryogenesis (SE) results from the transition of differentiated plant somatic cells into embryogenic cells that requires the extensive reprogramming of the somatic cell transcriptome. Commonly, the SE-involved genes are identified by analyzing the heterogeneous population of explant cells and thus, it is necessary to validate the expression of the candidate genes in the cells that are competent for embryogenic transition. Here, we optimized and implemented the whole mount in situ hybridization (WISH) method (Bleckmann and Dresselhaus, 2016; Dastidar et al., 2016) in order to analyze the spatiotemporal localization of miRNAs (miR156, miR166, miR390, miR167) and mRNAs such as WOX5 and PHABULOSA-target of miR165/166 during the SE that is induced in Arabidopsis explants. This study presents a detailed step-by-step description of the WISH procedure in which DIG-labeled LNA and RNA probes were used to detect miRNAs and mRNAs, respectively. The usefulness of the WISH in the functional analysis of the SE-involved regulatory pathways is demonstrated and the advantages of this method are highlighted: (i) the ability to analyze intact non-sectioned plant tissue; (ii) the specificity of transcript detection; (iii) the detection of miRNA; and (iv) a semi-quantitative assessment of the RNA abundance.

miR160 and miR166/165 Contribute to the LEC2-Mediated Auxin Response Involved in the Somatic Embryogenesis Induction in Arabidopsis

MicroRNAs are non-coding small RNA molecules that are involved in the post-transcriptional regulation of the genes that control various developmental processes in plants, including zygotic embryogenesis (ZE). miRNAs are also believed to regulate somatic embryogenesis (SE), a counterpart of the ZE that is induced in vitro in plant somatic cells. However, the roles of specific miRNAs in the regulation of the genes involved in SE, in particular those encoding transcription factors (TFs) with an essential function during SE including LEAFY COTYLEDON2 (LEC2), remain mostly unknown. The aim of the study was to reveal the function of miR165/166 and miR160 in the LEC2-controlled pathway of SE that is induced in in vitro cultured Arabidopsis explants.In ZE, miR165/166 controls the PHABULOSA/PHAVOLUTA (PHB/PHV) genes, which are the positive regulators of LEC2, while miR160 targets the AUXIN RESPONSE FACTORS (ARF10, ARF16, ARF17) that control the auxin signaling pathway, which plays key role in LEC2-mediated SE. We found that a deregulated expression/function of miR165/166 and miR160 resulted in a significant accumulation of auxin in the cultured explants and the spontaneous formation of somatic embryos. Our results show that miR165/166 might contribute to SE induction via targeting PHB, a positive regulator of LEC2 that controls embryogenic induction via activation of auxin biosynthesis pathway (Wójcikowska et al., 2013). Similar to miR165/166, miR160 was indicated to control SE induction through auxin-related pathways and the negative impact of miR160 on ARF10/ARF16/ARF17 was shown in an embryogenic culture. Altogether, the results suggest that the miR165/166- and miR160-node contribute to the LEC2-mediated auxin-related pathway of embryogenic transition that is induced in the somatic cells of Arabidopsis. A model summarizing the suggested regulatory interactions between the miR165/166-PHB and miR160-ARF10/ARF16/ARF17 nodes that control SE induction in Arabidopsis was proposed.

Endogenous small interfering RNAs associated with maize embryonic callus formation

The induction efficiency of maize embryonic callus is highly dependent on the genotype, and only a few lines possess a high capacity for callus formation. Although certain genes and pathways have been reported to contribute to the regulation of callus induction, to the best of our knowledge, the functions of the small interfering RNAs (siRNAs) involved in this process remain unknown. In this study, we identified 861 differentially expressed siRNAs and 576 target genes in the callus induction process. These target genes were classified into 3 clusters, and their functions involve controlling metalloexopeptidase activity, catalase activity, transcription regulation, and O-methyltransferase activity. In addition, certain genes related to auxin transport and stem cell or meristem development (e.g., PLT5-like, ARF15, SAUR-like, FAS1-like, Fea3, SCL5, and Zmwox2A) were regulated by the differentially expressed siRNAs. Moreover, zma-siR004119-2 directly cleaves the 5' UTR of Homeobox-transcription factor 25, which further leads to the down-regulation of its expression. Twelve 24-nt-siRNAs led to the hyper-methylation of GRMZM2G013465, which further decreases its expression. These results suggest that differentially expressed siRNAs regulate callus formation by controlling the expression of their target genes.

Identification of novel miRNAs and miRNA expression profiling in embryogenic tissues of Picea balfouriana treated by 6-benzylaminopurine

Here, we compared miRNA expression profiles in embryonic cell cultures of the conifer Picea balfouriana following application of the synthetic cytokinin 6-benzylaminopurine (6-BAP). We used next-generation sequencing to analyze three libraries of small RNAs from the treated embryogenic cell cultures and generated 24,000,000 raw reads from each of the libraries. Over 70 differentially regulated micro RNA (miRNA) families (≥2 fold change in expression) were identified between pairs of treatments. A quantitative analysis showed that miR3633 and miR1026 were upregulated in tissues with the highest embryogenic ability. These two miRNAs were predicted to target genes encoding receptor-like protein kinase and GAMYB transcription factors, respectively. In one library, miR1160, miR5638, miR1315, and miR5225 were downregulated. These four miRNAs were predicted to target genes encoding APETALA2, calmodulin-binding protein, and calcium-dependent protein kinase transcription factors. The expression patterns of the miRNAs and their targets were negatively correlated. Approximately 181 potentially novel P. balfouriana miRNAs were predicted from the three libraries, and seven were validated during the quantitative analysis. This study is the first report of differential miRNA regulation in tissues treated with 6-BAP during somatic embryogenesis. The differentially expressed miRNAs will be of value for investigating the mechanisms of embryogenic processes that are responsive to 6-BAP in P. balfouriana.

Small RNA and Transcriptome Sequencing Reveal a Potential miRNA-Mediated Interaction Network That Functions during Somatic Embryogenesis in Lilium pumilum DC. Fisch

Plant somatic embryos are widely used in the fields of germplasm conservation, breeding for genetic engineering and artificial seed production. MicroRNAs (miRNAs) play pivotal roles in somatic embryogenesis (SE) regulation. However, their regulatory roles during various stages of SE remain unclear. In this study, six types of embryogenic samples of Lilium pumilum DC. Fisch., including organogenic callus, embryogenic callus induced for 4 weeks, embryogenic callus induced for 6 weeks, globular embryos, torpedo embryos and cotyledon embryos, were prepared for small RNA sequencing. The results revealed a total of 2,378,760 small RNA reads, among which the most common size was 24 nt. Four hundred and fifty-two known miRNAs, belonging to more than 86 families, 57 novel miRNAs and 40 miRNA*s were identified. The 86 known miRNA families were sorted according to an alignment with their homologs across 24 land plants into the following four categories: 23 highly conserved, 4 moderately conserved, 15 less conserved and 44 species-specific miRNAs. Differentially expressed known miRNAs were identified during various stages of SE. Subsequently, the expression levels of 12 differentially expressed miRNAs and 4 targets were validated using qRT-PCR. In addition, six samples were mixed in equal amounts for transcript sequencing, and the sequencing data were used as transcripts for miRNA target prediction. A total of 66,422 unigenes with an average length of 800 bp were assembled from 56,258,974 raw reads. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment indicated that 38,004 and 15,497 unigenes were successfully assigned to GO terms and KEGG pathways, respectively. Among the unigenes, 2,182 transcripts were predicted to be targets for 396 known miRNAs. The potential targets of the identified miRNAs were mostly classified into the following GO terms: cell, binding and metabolic process. Enriched KEGG analysis demonstrated that carbohydrate metabolism was the predominant pathway in Lilium SE. Thus, we performed systemic characterization, homology comparisons and profiling of miRNA expression, and we constructed an miRNA-target network during Lilium SE for the first time. Our findings establish a foundation for the further exploration of critical genes and elucidation of SE in Lilium.

Small RNA profiling and degradome analysis reveal regulation of microRNA in peanut embryogenesis and early pod development

Background

As a typical geocarpic plant, peanut embryogenesis and pod development are complex processes involving many gene regulatory pathways and controlled by appropriate hormone level. MicroRNAs (miRNAs) are small non-coding RNAs that play indispensable roles in post-transcriptional gene regulation. Recently, identification and characterization of peanut miRNAs has been described. However, whether miRNAs participate in the regulation of peanut embryogenesis and pod development has yet to be explored.

Results

In this study, small RNA and degradome libraries from peanut early pod of different developmental stages were constructed and sequenced. A total of 70 known and 24 novel miRNA families were discovered. Among them, 16 miRNA families were legume-specific and 12 families were peanut-specific. 30 known and 10 novel miRNA families were differentially expressed during pod development. In addition, 115 target genes were identified for 47 miRNA families by degradome sequencing. Several new targets that might be specific to peanut were found and further validated by RNA ligase-mediated rapid amplification of 5′ cDNA ends (RLM 5′-RACE). Furthermore, we performed profiling analysis of intact and total transcripts of several target genes, demonstrating that SPL (miR156/157), NAC (miR164), PPRP (miR167 and miR1088), AP2 (miR172) and GRF (miR396) are actively modulated during early pod development, respectively.

Conclusions

Large numbers of miRNAs and their related target genes were identified through deep sequencing. These findings provided new information on miRNA-mediated regulatory pathways in peanut pod, which will contribute to the comprehensive understanding of the molecular mechanisms that governing peanut embryo and early pod development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3587-8) contains supplementary material, which is available to authorized users.

miR393 contributes to the embryogenic transition induced in vitro in Arabidopsis via the modification of the tissue sensitivity to auxin treatment

miR393 was found to control embryogenic transition in somatic cells in Arabidopsis via control of the TIR1 and AFB2 auxin receptors genes of the F-box family. miR393 molecules are believed to regulate the expression of the auxin receptors of the TAAR clade. Considering the central role of auxin in the induction of somatic embryogenesis (SE) in plant explants cultured in vitro, the involvement of miR393 in the embryogenic transition of somatic cells has been hypothesised. To verify this assumption, the reporter, overexpressor and mutant lines in genes encoded MIR393 and TIR1/AFB proteins of the F-box family were analysed during SE in Arabidopsis. Expression profiling of MIR393a and MIR393b, mature miR393 and the target genes (TIR1, AFB1, AFB2, AFB3) were investigated in explants undergoing SE. In addition, the embryogenic potential of various genotypes with a modified activity of the MIR393 and TIR1/AFB targets was evaluated. The distinct increase in the accumulation of miR393 that was coupled with a notable down-regulation of TIR1 and AFB2 targets was observed at the early phase of SE induction. Relevant to this observation, the GUS/GFP monitored expression of MIR393, TIR1 and AFB2 transcripts was localised in explant tissue undergoing SE induction. The results suggest the miR393-mediated regulation of TIR1 and AFB2 during embryogenic transition induced in Arabidopsis and a modification of the explant sensitivity to auxin treatment is proposed as underlying this regulatory pathway.

Identification of microRNAs and their corresponding targets involved in the susceptibility interaction of wheat response to Puccinia striiformis f. sp. tritici

MicroRNAs (miRNAs) play very important roles in plant defense responses. However, little is known about their roles in the susceptibility interaction between wheat and Puccinia striiformis f. sp. tritici (Pst). In this study, two miRNA libraries were constructed from the leaves of the cultivar Xingzi 9104 inoculated with the virulent Pst race CYR32 and sterile water, respectively. A total of 1316 miRNA candidates, including 173 known miRNAs that were generated from 98 pre-miRNAs, were obtained. The remaining 1143 miRNA candidates included 145 conserved and 998 wheat-specific miRNAs that were generated from 87 and 1088 pre-miRNAs, respectively. The 173 known and 145 conserved miRNAs were sub-classified into 63 miRNA families. The target genes of wheat miRNAs were also confirmed using degradome sequencing technology. Most of the annotated target genes were related to signal transduction or energy metabolism. Additionally, we found that miRNAs and their target genes form complicated regulation networks. The expression profiles of miRNAs and their corresponding target genes were further analyzed by quantitative real-time polymerase chain reaction (qRT-PCR), and the results indicate that some miRNAs are involved in the compatible wheat-Pst susceptibility interaction. Importantly, tae-miR1432 was highly expressed when wheat was challenged with CYR32, and the corresponding target gene, predicted to be a calcium ion-binding protein, also exhibited upregulated expression but a divergent expression trend. PC-3P-7484, a specific wheat miRNA, was highly expressed in the wheat response to Pst infection, while the expression of the corresponding target gene ubiquillin was dramatically downregulated. These data provide the foundation for evaluating the important regulatory roles of miRNAs in wheat-Pst susceptibility interaction.

De novo comparative transcriptome analysis provides new insights into sucrose induced somatic embryogenesis in camphor tree (Cinnamomum camphora L.)

BACKGROUND

Somatic embryogenesis is a notable illustration of cell totipotency, by which somatic cells undergo dedifferentiation and then differentiate into somatic embryos. Our previous work demonstrated that pretreatment of immature zygotic embryos with 0.5 M sucrose solution for 72 h efficiently induced somatic embryo initiation in camphor tree. To better understand the molecular basis of somatic embryogenesis induced by osmotic stress, de novo transcriptome sequencing of three tissues of camphor tree (immature zygotic embryos, sucrose-pretreated immature zygotic embryos, and somatic embryos induced from sucrose-pretreated zygotic embryos) were conducted using Illumina Hiseq 2000 platform.

RESULTS

A total of 30.70 G high quality clean reads were obtained from cDNA libraries of the three samples. The overall de novo assembly of cDNA sequence data generated 205592 transcripts, with an average length of 998 bp. 114229 unigenes (55.56 % of all transcripts) with an average length of 680 bp were annotated with gene descriptions, gene ontology terms or metabolic pathways based on Blastx search against Nr, Nt, Swissprot, GO, COG/KOG, and KEGG databases. CEGMA software identified 237 out of 248 ultra-conserved core proteins as 'complete' in the transcriptome assembly, showing a completeness of 95.6 %. A total of 897 genes previously annotated to be potentially involved in somatic embryogenesis were identified. Comparative transcriptome analysis showed that a total of 3335 genes were differentially expressed in the three samples. The differentially expressed genes were divided into six groups based on K-means clustering. Expression level analysis of 52 somatic embryogenesis-related genes indicated a high correlation between RNA-seq and qRT-PCR data. Gene enrichment analysis showed significantly differential expression of genes responding to stress and stimulus.

CONCLUSIONS

The present work reported a de novo transcriptome assembly and global analysis focused on gene expression changes during initiation and formation of somatic embryos in camphor tree. Differential expression of somatic embryogenesis-related genes indicates that sucrose induced somatic embryogenesis may share or partly share the mechanisms of somatic embryogenesis induced by plant hormones. This study provides comprehensive transcript information and gene expression data for camphor tree. It could also serve as an important platform resource for further functional studies in plant embryogenesis.

MicroRNA functions in plant embryos

Plant miRNAs are short non-coding RNAs that mediate the repression of hundreds of genes. The basic plant body plan is established during early embryogenesis, and recent results have demonstrated that miRNAs play pivotal roles during both embryonic pattern formation and developmental timing. Multiple miRNAs appear to specifically repress transcription factor families during early embryogenesis. Therefore miRNAs probably have a large influence on the gene regulatory networks that contribute to the earliest cellular differentiation events in plants.

Regulation of RNA metabolism is important for in vitro dedifferentiation of plant cells

The characteristic high regenerative ability of plants has been exploited to develop in vitro plant regeneration techniques, which are usually initiated by an in vitro dedifferentiation step induced by artificial phytohormone treatment. Recent advances in plant molecular biological and genetic technologies have revealed the importance of the regulation of RNA metabolism, including the control of rRNA biosynthesis, pre-mRNA splicing, and miRNA-based RNA decay, in successful in vitro dedifferentiation. This review provides a brief overview of current knowledge of the roles of RNA metabolism in the dedifferentiation of plant cells in vitro. In addition, the possibility that plant-specific aspects of RNA metabolism regulation are linked closely to their high regenerative ability is discussed.

Genomewide analysis of small RNAs in nonembryogenic and embryogenic tissues of citrus: microRNA- and siRNA-mediated transcript cleavage involved in somatic embryogenesis

Somatic embryogenesis (SE) is a process of somatic cells becoming dedifferentiated and generating embryos. SE has been widely used in biotechnology as a powerful way of regeneration and a model system for studying plant embryogenesis, but the controlling mechanisms of SE are far from clear. Here, we show the genomewide profiles of miRNAs/siRNAs and their target genes in nonembryogenic and embryogenic tissues of 'Valencia' sweet orange. By high-throughput sequencing (HTS) of small RNAs and RNA degradome tags, we identified 50 known and 45 novel miRNAs, 130 miniature inverted-repeat transposable elements (MITEs) derived, 94 other and 235 phased small interfering RNAs (siRNAs), as well as 203 target genes. The majority of the abundantly expressed miRNAs/siRNAs exhibit lower expression levels in embryogenic callus (EC) or during SE process than in nonembryogenic callus (NEC), which is supposed to derepress the target genes that are involved in development and stress response, thus to activate the biological processes required for cell differentiation. However, the conserved csi-miR156a/b, miR164b and 171c directed suppression of specific transcription factors (TFs) are supposed to inactivate the postembryonic growth thus to maintain normal SE. In this study, miRNA- and siRNA-mediated silencing of target genes was found under sophisticated regulation in citrus SE system; the enhancement effect of specific conserved miRNAs on SE was discussed, providing new clues for future investigation of mechanisms that control SE.

Maize miRNA and target regulation in response to hormone depletion and light exposure during somatic embryogenesis

Maize somatic embryogenesis (SE) is induced from the immature zygotic embryo in darkness and under the appropriate hormones' levels. Small RNA expression is reprogrammed and certain miRNAs become particularly enriched during induction while others, characteristic to the zygotic embryo, decrease. To explore the impact of different environmental cues on miRNA regulation in maize SE, we tested specific miRNA abundance and their target gene expression in response to photoperiod and hormone depletion for two different maize cultivars (VS-535 and H-565). The expression levels of miR156, miR159, miR164, miR168, miR397, miR398, miR408, miR528, and some predicted targets (SBP23, GA-MYB, CUC2, AGO1c, LAC2, SOD9, GR1, SOD1A, PLC) were examined upon staged hormone depletion in the presence of light photoperiod or darkness. Almost all examined miRNA, except miR159, increased upon hormone depletion, regardless photoperiod absence/presence. miR528, miR408, and miR398 changed the most. On the other hand, expression of miRNA target genes was strongly regulated by the photoperiod exposure. Stress-related miRNA targets showed greater differences between cultivars than development-related targets. miRNA/target inverse relationship was more frequently observed in darkness than light. Interestingly, miR528, but not miR159, miR168 or miR398, was located on polyribosome fractions suggesting a role for this miRNA at the level of translation. Overall our results demonstrate that hormone depletion exerts a great influence on specific miRNA expression during plant regeneration independently of light. However, their targets are additionally influenced by the presence of photoperiod. The reproducibility or differences observed for particular miRNA-target regulation between two different highly embryogenic genotypes provide clues for conserved miRNA roles within the SE process.

Differential miRNA expression profiles in the longissimus dorsi muscle between intact and castrated male pigs

MicroRNAs (miRNAs) are important modulators of skeletal muscle development in multiple mammalian species, but their role in skeletal muscle growth in castrated male pigs has not been well studied. The aim of the present study was to determine the role of miRNAs in longissimus dorsi muscle under castration. Our results showed that castration caused a significant decrease in serum testosterone levels as well as carcass lean mass, but led to an increase in carcass fat mass. Moreover, miRNA expression profiles in skeletal muscle were significantly altered by castration, and seven differentially expressed miRNAs were discovered. More importantly, functional analysis suggested that these differentially expressed miRNAs and their targets are involved in the regulation of skeletal muscle contractile function and fat metabolism. Taken together, these results demonstrate altered miRNA expression in skeletal muscle of castrated male pigs, and suggest a potential mechanism underlying the effects of castration on porcine skeletal muscle growth.

Analysis of miRNAs and their targets during adventitious shoot organogenesis of Acacia crassicarpa

Organogenesis is an important process for plant regeneration by tissue or cell mass differentiation to regenerate a complete plant. MicroRNAs (miRNAs) play an essential role in regulating plant development by mediating target genes at transcriptional and post-transcriptional levels, but the diversity of miRNAs and their potential roles in organogenesis of Acacia crassicarpa have rarely been investigated. In this study, approximately 10 million sequence reads were obtained from a small RNA library, from which 189 conserved miRNAs from 57 miRNA families, and 7 novel miRNAs from 5 families, were identified from A. crassicarpa organogenetic tissues. Target prediction for these miRNAs yielded 237 potentially unique genes, of which 207 received target Gene Ontology annotations. On the basis of a bioinformatic analysis, one novel and 13 conserved miRNAs were selected to investigate their possible roles in A. crassicarpa organogenesis by qRT-PCR. The stage-specific expression patterns of the miRNAs provided information on their possible regulatory functions, including shoot bud formation, modulated function after transfer of the culture to light, and regulatory roles during induction of organogenesis. This study is the first to investigate miRNAs associated with A. crassicarpa organogenesis. The results provide a foundation for further characterization of miRNA expression profiles and roles in the regulation of diverse physiological pathways during adventitious shoot organogenesis of A. crassicarpa.

Conservation and divergence of transcriptomic and epigenomic variation in maize hybrids

Background

Recent genome-wide studies suggested that in addition to genetic variations, epigenetic variations may also be associated with differential gene expression and growth vigor in plant hybrids. Maize is an ideal model system for the study of epigenetic variations in hybrids given the significant heterotic performance, the well-known complexity of the genome, and the rich history in epigenetic studies. However, integrated comparative transcriptomic and epigenomic analyses in different organs of maize hybrids remain largely unexplored.

Results

Here, we generated integrated maps of transcriptomes and epigenomes of shoots and roots of two maize inbred lines and their reciprocal hybrids, and globally surveyed the epigenetic variations and their relationships with transcriptional divergence between different organs and genotypes. We observed that whereas histone modifications vary both between organs and between genotypes, DNA methylation patterns are more distinguishable between genotypes than between organs. Histone modifications were associated with transcriptomic divergence between organs and between hybrids and parents. Further, we show that genes up-regulated in both shoots and roots of hybrids were significantly enriched in the nucleosome assembly pathway. Interestingly, 22- and 24-nt siRNAs were shown to be derived from distinct transposable elements, and for different transposable elements in both shoots and roots, the differences in siRNA activity between hybrids and patents were primarily driven by different siRNA species.

Conclusions

These results suggest that despite variations in specific genes or genomic loci, similar mechanisms may account for the genome-wide epigenetic regulation of gene activity and transposon stability in different organs of maize hybrids.