Genomic organization, differential expression, and interaction of SQUAMOSA promoter-binding-like transcription factors and microRNA156 in rice

SPL6 represses signalling outputs of ER stress in control of panicle cell death in rice

Inositol-requiring enzyme 1 (IRE1) is the most conserved transducer of the unfolded protein response that produces either adaptive or death signals depending on the amplitude and duration of its activation. Here, we report that SQUAMOSA PROMOTER-BINDING PROTEIN-LIKE 6 (SPL6)-deficient plants displayed hyperactivation of the endoplasmic reticulum (ER) stress sensor IRE1, leading to cell death in rice panicles, indicating that SPL6 is an essential survival factor for the suppression of persistent or intense ER stress conditions. Importantly, knockdown of the hyperactivated mRNA level of IRE1 rescues panicle apical abortion in the spl6-1 transgenic plants harbouring the IRE1-RNAi constructs, establishing the genetic linkage between the hyperactivation of IRE1 and cell death in spl6-1. Our findings reveal a novel cell survival machinery in which SPL6 represses the transcriptional activation of the ER stress sensor IRE1 in control of ER stress signalling outputs that hinge on a balance between adaptive and death signals for determining cell fates during ER stress.

Analysis of microRNA reveals cleistogamous and chasmogamous floret divergence in dimorphic plant

Cleistogenes songorica, a grass species that exhibits two spatially different type of inflorescence, chastogamy (CH), flowers localized at the top, and cleistogamy (CL) flowers embedded in leaf sheath. This study aimed at dissecting reasons underlying these distinct floral development patterns at morphological and microRNA level. Phenotyping for CH and CL was conducted and four small RNA libraries were constructed from the CH and CL flowers for high-throughput sequencing to identify the differentiated miRNAs. As results, spikelet, stigma, anther, lemma and lodicule length of CH flowers were found larger than that of CL, and so was seed setting. Also, 17 flower-related differential expression miRNAs were identified which were associated with floral organ development and morphogenesis, and the flower development. Further results showed that miR159a.1-CL3996.Contig2 pair was related to anther development, miR156a-5p-CL1954.Contig2 was linked to response to high light intensity, miR408-3p/miR408d-Unigene429 was related to pollination and Unigene429 positively regulated flower development. To our knowledge, this is the first study on differential miRNA accumulation between CH and CL flowers and our study serves as a foundation to the future elucidation of regulatory mechanisms of miRNAs in the divergent development of CL and CH flowers in a single plant.

From model to crop: functional characterization of SPL8 in M. truncatula led to genetic improvement of biomass yield and abiotic stress tolerance in alfalfa

Biomass yield, salt tolerance and drought tolerance are important targets for alfalfa (Medicago sativa L.) improvement. Medicago truncatula has been developed into a model plant for alfalfa and other legumes. By screening a Tnt1 retrotransposon-tagged M. truncatula mutant population, we identified three mutants with enhanced branching. Branch development determines shoot architecture which affects important plant functions such as light acquisition, resource use and ultimately impacts biomass production. Molecular analyses revealed that the mutations were caused by Tnt1 insertions in the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 8 (SPL8) gene. The M. truncatula spl8 mutants had increased biomass yield, while overexpression of SPL8 in M. truncatula suppressed branching and reduced biomass yield. Scanning electron microscopy (SEM) analysis showed that SPL8 inhibited branching by directly suppressing axillary bud formation. Based on the M. truncatula SPL8 sequence, alfalfa SPL8 (MsSPL8) was cloned and transgenic alfalfa plants were produced. MsSPL8 down-regulated or up-regulated alfalfa plants exhibited similar phenotypes to the M. truncatula mutants or overexpression lines, respectively. Specifically, the MsSPL8 down-regulated alfalfa plants showed up to 43% increase in biomass yield in the first harvest. The impact was even more prominent in the second harvest, with up to 86% increase in biomass production compared to the control. Furthermore, down-regulation of MsSPL8 led to enhanced salt and drought tolerance in transgenic alfalfa. Results from this research offer a valuable approach to simultaneously improve biomass production and abiotic stress tolerance in legumes.

Overexpression of a SBP-Box Gene (VpSBP16) from Chinese Wild Vitis Species in Arabidopsis Improves Salinity and Drought Stress Tolerance

Salinity and drought are two major abiotic stresses that limit grape productivity. Responses to stress in grape are known to be regulated by several families of transcription factors. However, little is known about the role of grape Squamosa promoter binding protein (SBP)-box transcription factor genes in response to abiotic stress. To better understand the functions of the grape SBP-box genes in abiotic stress tolerance, a full-length complementary DNA (cDNA) sequence of the putative SBP-box transcription factor gene, VpSBP16 was amplified from Chinese wild grapevine Vitis pseudoreticulata clone "Baihe-35-1". We observed that the VpSBP16 protein fused to the green fluorescent protein (GFP) reporter accumulated in the nucleus when transiently expressed in onion epidermal cells. Moreover, VpSBP16 was shown to have transcriptional activation activity using a yeast trans-activation assay. We performed a VpSBP16 functional analysis through the characterization of transgenic Arabidopsis thaliana plants constitutively over-expressing VpSBP16. The transgenic lines had longer roots and the seeds had a higher germination rate than the wild type (WT) under osmotic stress. In addition, the accumulation of reactive oxygen species (ROS) of transgenic seedlings was significantly lower than WT in the transgenic lines, as was electrolyte leakage. VpSBP16 overexpression also elevated expression levels of stress-response genes involved in the salt overly sensitive (SOS) pathway. These results indicate that overexpression VpSBP16 in A. thaliana enhances tolerance of salt and drought stress during seed germination, as well in seedlings and mature plants, by regulating SOS and ROS signaling cascades.

Genome-wide identification and characterization of the SBP-box gene family in Petunia

Background

SQUAMOSA PROMOTER BINDING PROTEIN (SBP)-box genes encode a family of plant-specific transcription factors (TFs) that play important roles in many growth and development processes including phase transition, leaf initiation, shoot and inflorescence branching, fruit development and ripening etc. The SBP-box gene family has been identified and characterized in many species, but has not been well studied in Petunia, an important ornamental genus.

Results

We identified 21 putative SPL genes of Petunia axillaris and P. inflata from the reference genome of P. axillaris N and P. inflata S6, respectively, which were supported by the transcriptome data. For further confirmation, all the 21 genes were also cloned from P. hybrida line W115 (Mitchel diploid). Phylogenetic analysis based on the highly conserved SBP domains arranged PhSPLs in eight groups, analogous to those from Arabidopsis and tomato. Furthermore, the Petunia SPL genes had similar exon-intron structure and the deduced proteins contained very similar conserved motifs within the same subgroup. Out of 21 PhSPL genes, fourteen were predicted to be potential targets of PhmiR156/157, and the putative miR156/157 response elements (MREs) were located in the coding region of group IV, V, VII and VIII genes, but in the 3’-UTR regions of group VI genes. SPL genes were also identified from another two wild Petunia species, P. integrifolia and P. exserta, based on their transcriptome databases to investigate the origin of PhSPLs. Phylogenetic analysis and multiple alignments of the coding sequences of PhSPLs and their orthologs from wild species indicated that PhSPLs were originated mainly from P. axillaris. qRT-PCR analysis demonstrated differential spatiotemperal expression patterns of PhSPL genes in petunia and many were expressed predominantly in the axillary buds and/or inflorescences. In addition, overexpression of PhSPL9a and PhSPL9b in Arabidopsis suggested that these genes play a conserved role in promoting the vegetative-to-reproductive phase transition.

Conclusion

Petunia genome contains at least 21 SPL genes, and most of the genes are expressed in different tissues. The PhSPL genes may play conserved and diverse roles in plant growth and development, including flowering regulation, leaf initiation, axillary bud and inflorescence development. This work provides a comprehensive understanding of the SBP-box gene family in Petunia and lays a significant foundation for future studies on the function and evolution of SPL genes in petunia.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4537-9) contains supplementary material, which is available to authorized users.

Genome-Wide Analysis Characterization and Evolution of SBP Genes in Fragaria vesca, Pyrus bretschneideri, Prunus persica and Prunus mume

The SQUAMOSA promoter binding protein (SBP)-box proteins are plant-specific transcriptional factors in plants. SBP TFs are known to play important functions in a diverse development process and also related in the process of evolutionary novelties. SBP gene family has been characterized in several plant species, but little is known about molecular evolution, functional divergence and comprehensive study of SBP gene family in Rosacea. We carried out genome-wide investigations and identified 14, 32, 17, and 17 SBP genes from four Rosacea species (Fragaria vesca, Pyrus bretschneideri, Prunus persica and Prunus mume, respectively). According to phylogenetic analysis arranged the SBP protein sequences in seven groups. Localization of SBP genes presented an uneven distribution on corresponding chromosomes of Rosacea species. Our analyses designated that the SBP genes duplication events (segmental and tandem) and divergence. In addition, due to highly conserved structure pattern of SBP genes, recommended that highly conserved region of microsyneteny in the Rosacea species. Type I and II functional divergence was detected among various amino acids in SBP proteins, while there was no positive selection according to substitutional model analysis using PMAL software. These results recommended that the purifying selection might be leading force during the evolution process and dominate conservation of SBP genes in Rosacea species according to environmental selection pressure analysis. Our results will provide basic understanding and foundation for future research insights on the evolution of the SBP genes in Rosacea.

Small but powerful: function of microRNAs in plant development

MicroRNAs (miRNAs) are a group of endogenous noncoding small RNAs frequently 21 nucleotides long. miRNAs act as negative regulators of their target genes through sequence-specific mRNA cleavage, translational repression, or chromatin modifications. Alterations of the expression of a miRNA or its targets often result in a variety of morphological and physiological abnormalities, suggesting the strong impact of miRNAs on plant development. Here, we review the recent advances on the functional studies of plant miRNAs. We will summarize the regulatory networks of miRNAs in a series of developmental processes, including meristem development, establishment of lateral organ polarity and boundaries, vegetative and reproductive organ growth, etc. We will also conclude the conserved and species-specific roles of plant miRNAs in evolution and discuss the strategies for further elucidating the functional mechanisms of miRNAs during plant development.

Regulation of Long Noncoding RNAs Responsive to Phytoplasma Infection in Paulownia tomentosa

Paulownia witches' broom caused by phytoplasma infection affects the production of Paulownia trees worldwide. Emerging evidence showed that long noncoding RNAs (lncRNA) play a protagonist role in regulating the expression of genes in plants. So far, the identification of lncRNAs has been limited to a few model plant species, and their roles in mediating responses to Paulownia tomentosa that free of phytoplasma infection are yet to be characterized. Here, whole-genome identification of lncRNAs, based on strand-specific RNA sequencing, from four Paulownia tomentosa samples, was performed and identified 3689 lncRNAs. These lncRNAs showed low conservation among plant species and some of them were miRNA precursors. Further analysis revealed that the 112 identified lncRNAs were related to phytoplasma infection. We predicted the target genes of these phytoplasma-responsive lncRNAs, and our analysis showed that 51 of the predicted target genes were alternatively spliced. Moreover, we found the expression of the lncRNAs plays vital roles in regulating the genes involved in the reactive oxygen species induced hypersensitive response and effector-triggered immunity in phytoplasma-infected Paulownia. This study indicated that diverse sets of lncRNAs were responsive to Paulownia witches' broom, and the results will provide a starting point to understand the functions and regulatory mechanisms of Paulownia lncRNAs in the future.

Genome-wide identification and characterization of SPL transcription factor family and their evolution and expression profiling analysis in cotton

Plant specific transcription factors, SQUAMOSA promoter-binding protein-like (SPL), are involved in many biological processes. However, no systematical study has been reported in cotton. In this study, a total of 177 SPL genes were identified, including 29, 30, 59 and 59 SPLs in Gossypium arboreum, G. raimondii, G. barbadense, and G. hirsutum, respectively. These SPL genes were classified into eight phylogenetical groups. The gene structure, conserved motif, and clustering were highly conserved within each orthologs. Two zinc finger-like structures (Cys3His and Cys2HisCys) and NLS segments were existed in all GrSPLs. Segmental duplications play important roles in SPL family expansion, with 20 genes involved in segmental duplications and 2 in tandem duplications, and ten ortholog pairs in syntenic regions between G. raimondii and A. thaliana. Several putative cis-elements, involved in light, stresses and phytohormones response, were found in the promoter regions of GhSPLs, suggesting that plant responses to those environmental changes may be induced through targeting SPL transcription factors. RNA-seq analysis shows that SPL genes were differentially expressed in cotton; some were highly expressed during fiber initiation and early development. Comparing with other plants, SPL genes show subfunctionalization, lost and/or gain functions in cotton during long-term domestication and evolution.

Molecular Characterization of SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) Gene Family in Betula luminifera

As a major family of plant-specific transcription factors, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes play vital regulatory roles in plant growth, development and stress responses. In this study, 18 SPL genes were identified and cloned from Betula luminifera. Two zinc finger-like structures and a nuclear location signal (NLS) segments were existed in the SBP domains of all BlSPLs. Phylogenetic analysis showed that these genes were clustered into nine groups (group I-IX). The intron/exon structure and motif composition were highly conserved within the same group. 12 of the 18 BlSPLs were experimentally verified as the targets of miR156, and two cleavage sites were detected in these miR156-targeted BlSPL genes. Many putative cis-elements, associated with light, stresses and phytohormones response, were identified in the promoter regions of BlSPLs, suggesting that BlSPL genes are probably involved in important physiological processes and developmental events. Tissue-specific expression analysis showed that miR156-targeted BlSPLs exhibited a more differential expression pattern, while most miR156-nontargeted BlSPLs tended to be constitutively expressed, suggesting the distinct roles of miR156-targeted and nontargeted BlSPLs in development and growth of B. luminifera. Further expression analysis revealed that miR156-targeted BlSPLs were dramatically up-regulated with age, whereas mature BlmiR156 level was apparently declined with age, indicating that miR156/SPL module plays important roles in vegetative phase change of B. luminifera. Moreover, yeast two-hybrid assay indicated that several miR156-targeted and nontargeted BlSPLs could interact with two DELLA proteins (BlRGA and BlRGL), which suggests that certain BlSPLs take part in the GA regulated processes through protein interaction with DELLA proteins. All these results provide an important basis for further exploring the biological functions of BlSPLs in B. luminifera.

Under a New Light: Regulation of Light-Dependent Pathways by Non-coding RNAs

The biological relevance of non-protein coding RNAs in the regulation of critical plant processes has been firmly established in recent years. This has been mostly achieved with the discovery and functional characterization of small non-coding RNAs, such as small interfering RNAs and microRNAs (miRNAs). However, recent next-generation sequencing techniques have widened our view of the non-coding RNA world, which now includes long non-coding RNAs (lncRNAs). Small and lncRNAs seem to diverge in their biogenesis and mode of action, but growing evidence highlights their relevance in developmental processes and in responses to particular environmental conditions. Light can affect MIRNA gene transcription, miRNA biogenesis, and RNA-induced silencing complex (RISC) activity, thus controlling not only miRNA accumulation but also their biological function. In addition, miRNAs can mediate several light-regulated processes. In the lncRNA world, few reports are available, but they already indicate a role in the regulation of photomorphogenesis, cotyledon greening, and photoperiod-regulated flowering. In this review, we will discuss how light controls MIRNA gene expression and the accumulation of their mature forms, with a particular emphasis on those miRNAs that respond to different light qualities and are conserved among species. We will also address the role of small non-coding RNAs, particularly miRNAs, and lncRNAs in the regulation of light-dependent pathways. We will mainly focus on the recent progress done in understanding the interconnection between these non-coding RNAs and photomorphogenesis, circadian clock function, and photoperiod-dependent flowering.

Epigenetic Regulation of Juvenile-to-Adult Transition in Plants

Epigenetic regulation is referred to as changes in gene function that do not involve changes in the DNA sequence, it is usually accomplished by DNA methylation, histone modifications (repressive marks such as H3K9me, H3K27me, H2Aub, or active marks such as H3K4me, H3K36me, H3Ac), and chromatin remodeling (nucleosome composition, occupancy, and location). In plants, the shoot apex produces different lateral organs during development to give rise to distinguishable phases of a juvenile, an adult and a reproductive phase after embryogenesis. The juvenile-to-adult transition is a key developmental event in plant life cycle, and it is regulated by a decrease in the expression of a conserved microRNA-miR156/157, and a corresponding increase in the expression of its target genes encoding a set of plant specific SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) proteins. Recent work has revealed that the miR156/157-SPL pathway is the master regulator of juvenile-to-adult transition in plants, and genes in this pathway are subjected to epigenetic regulation, such as DNA methylation, histone modifications, and chromatin remodeling. In this review, we summarized the recent progress in understanding the epigenetic regulation of the miR156/157-SPL pathway during juvenile-to-adult transition and bring forward some perspectives of future research in this field.

SPL13 regulates shoot branching and flowering time in Medicago sativa

Our results show SPL13 plays a crucial role in regulating vegetative and reproductive development in Medicago sativa L. (alfalfa), and that MYB112 is targeted and downregulated by SPL13 in alfalfa. We previously showed that transgenic Medicago sativa (alfalfa) plants overexpressing microRNA156 (miR156) show a bushy phenotype, reduced internodal length, delayed flowering time, and enhanced biomass yield. In alfalfa, transcripts of seven SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors, including SPL13, are targeted for cleavage by miR156. Thus, association of each target SPL gene to a trait or set of traits is essential for developing molecular markers for alfalfa breeding. In this study, we investigated SPL13 function using SPL13 overexpression and silenced alfalfa plants. Severe growth retardation, distorted branches and up-curled leaves were observed in miR156-impervious 35S::SPL13m over-expression plants. In contrast, more lateral branches and delayed flowering time were observed in SPL13 silenced plants. SPL13 transcripts were predominantly present in the plant meristems, indicating that SPL13 is involved in regulating shoot branch development. Accordingly, the shoot branching-related CAROTENOID CLEAVAGE DIOXYGENASE 8 gene was found to be significantly downregulated in SPL13 RNAi silencing plants. A R2R3-MYB gene MYB112 was also identified as being directly silenced by SPL13 based on Next Generation Sequencing-mediated transcriptome analysis and chromatin immunoprecipitation assays, suggesting that MYB112 may be involved in regulating alfalfa vegetative growth.

Hybrid incompatibilities in interspecific crosses between tetraploid wheat and its wild diploid relative Aegilops umbellulata

Hybrid abnormalities, severe growth abortion and grass-clump dwarfism, were found in the tetraploid wheat/Aegilops umbellulata hybrids, and the gene expression changes were conserved in the hybrids with those in other wheat synthetic hexaploids. Aegilops umbellulata Zhuk., a diploid goatgrass species with a UU genome, has been utilized as a genetic resource for wheat breeding. Here, we examine the reproductive barriers between tetraploid wheat cultivar Langdon (Ldn) and various Ae. umbellulata accessions by conducting interspecific crossings. Through systematic cross experiments, three types of hybrid incompatibilities were found: seed production failure in crosses, hybrid growth abnormalities and sterility in the ABU hybrids. Hybrid incompatibilities were widely distributed over the entire range of the natural species, and in about 50% of the cross combinations between tetraploid Ldn and Ae. umbellulata accessions, ABU F₁ hybrids showed one of two abnormal growth phenotypes: severe growth abortion (SGA) or grass-clump dwarfism. Expression of the shoot meristem maintenance-related and cell cycle-related genes was markedly repressed in crown tissues of hybrids showing SGA, suggesting dysfunction of mitotic cell division in the shoot apices. The grass-clump dwarf phenotype may be explained by down-regulation of wheat APETALA1-like MADS box genes, which act as flowering promoters, and altered expression in crown tissues of the miR156/SPLs module, which controls tiller number and branching. These gene expression changes in growth abnormalities were well conserved between the Ldn/Ae. umbellulata plants and interspecific hybrids from crosses of Ldn and wheat D-genome progenitor Ae. tauschii.

Small RNAs: Big Impact on Plant Development

While the role of proteins in determining cell identity has been extensively studied, the contribution of small noncoding RNA molecules such as miRNAs and siRNAs has been also recognised. miRNAs bind to complementary sites in target mRNA molecules to trigger the degradation or translational inhibition of those targets. Recent studies have revealed that miRNAs play pivotal roles in key developmental processes such as patterning of the embryo, meristem, leaf, and flower. Furthermore, these miRNAs have been recruited throughout plant evolution into pathways that create diverse plant organ forms and shapes. This review focuses on the roles of miRNAs in establishing plant cell identity during key plant development processes and creating morphological diversity during plant evolution.

Identification of miRNA-mediated drought responsive multi-tiered regulatory network in drought tolerant rice, Nagina 22

Comparative characterization of microRNA-mediated stress regulatory networks in contrasting rice cultivars is critical to decipher plant stress response. Consequently, a multi-level comparative analysis, using sRNA sequencing, degradome analysis, enzymatic and metabolite assays and metal ion analysis, in drought tolerant and sensitive rice cultivars was conducted. The study identified a group of miRNAs "Cultivar-specific drought responsive" (CSDR)-miRNAs (osa-miR159f, osa-miR1871, osa-miR398b, osa-miR408-3p, osa-miR2878-5p, osa-miR528-5p and osa-miR397a) that were up-regulated in the flag-leaves of tolerant cultivar, Nagina 22 (N22) and Vandana, but down-regulated in the sensitive cultivar, Pusa Basmati 1 (PB1) and IR64, during drought. Interestingly, CSDR-miRNAs target several copper-protein coding transcripts like plantacyanins, laccases and Copper/Zinc superoxide dismutases (Cu/Zn SODs) and are themselves found to be similarly induced under simulated copper-starvation in both N22 and PB1. Transcription factor OsSPL9, implicated in Cu-homeostasis also interacted with osa-miR408-3p and osa-miR528-5p promoters. Further, N22 flag leaves showed lower SOD activity, accumulated ROS and had a higher stomata closure. Interestingly, compared to PB1, internal Cu levels significantly decreased in the N22 flag-leaves, during drought. Thus, the study identifies the unique drought mediated dynamism and interplay of Cu and ROS homeostasis, in the flag leaves of drought tolerant rice, wherein CSDR-miRNAs play a pivotal role.

Phenanthrene-responsive microRNAs and their targets in wheat roots

MicroRNAs (miRNAs) play key roles in plant growth, development and responses to abiotic stress. Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental pollutants. However, it is yet unknown how miRNAs work during PAH uptake by plant roots. Thus, in this study we ascertain phenanthrene (a model PAH)-responsive miRNAs using small RNA high-throughput deep sequencing and their target genes in wheat roots. We identified 108 conserved and non-conserved miRNA members belonging to 82 miRNA families and found 11 differentially expressed miRNAs, among which four miRNAs (miR156, miR164, miR171a and miR9678-3p) were up-regulated and the other seven miRNAs (miR398, miR531, miR1121, miR5048-5p, miR9653b, miR9773 and miR9778) were down-regulated. ABC-transporter-related Gene CA704421 and CA697226 did not respond to phenanthrene exposure. miR156 and miR164 might regulate directly the growth and development of wheat roots by targeting SPL and NAC, respectively. miR398 and miR1121 could regulate oxidative reactions to respond to phenanthrene stress. Additionally, miR9773 might involve phenanthrene metabolism through acting on CYP450. Therefore, it is concluded that phenanthrene triggers variation in miRNA expression, which is associated with uptake of and response to phenanthrene. These findings are of significance for further understanding miRNA regulation mechanisms on PAH uptake, and providing guidance for screening of resistant cultivars in crop production and phytoremediation of PAH-contaminated soils or water at genetic level.

The miR156-SPL4 module predominantly regulates aerial axillary bud formation and controls shoot architecture

Grasses possess basal and aerial axillary buds. Previous studies have largely focused on basal bud (tiller) formation but scarcely touched on aerial buds, which may lead to aerial branch development. Genotypes with and without aerial buds were identified in switchgrass (Panicum virgatum), a dedicated bioenergy crop. Bud development was characterized using scanning electron microscopy. Microarray, RNA-seq and quantitative reverse transcription polymerase chain reaction (RT-qPCR) were used to identify regulators of bud formation. Gene function was characterized by down-regulation and overexpression. Overexpression of miR156 induced aerial bud formation in switchgrass. Various analyses revealed that SQUAMOSA PROMOTER BINDING PROTEIN LIKE4 (SPL4), one of the miR156 targets, directly regulated aerial axillary bud initiation. Down-regulation of SPL4 promoted aerial bud formation and increased basal buds, while overexpression of SPL4 seriously suppressed bud formation and tillering. RNA-seq and RT-qPCR identified potential downstream genes of SPL4. Unlike all previously reported genes acting as activators of basal bud initiation, SPL4 acts as a suppressor for the formation of both aerial and basal buds. The miR156-SPL4 module predominantly regulates aerial bud initiation and partially controls basal bud formation. Genetic manipulation of SPL4 led to altered plant architecture with increased branching, enhanced regrowth after cutting and improved biomass yield.

Comparative analysis of molecular and physiological traits between perennial Arabis alpina Pajares and annual Arabidopsis thaliana Sy-0

Annual plants complete life cycle in a year while perennial plants maintain growth for several years. Arabis alpina, a polycarpic perennial, is a close relative of monocarpic annual Arabidopsis. Pajares is an accession of A. alpina requiring vernalization, a long-term cold for flowering. Arabidopsis shows holistic flowering whereas Pajares shows idiographic flowering, producing axillary branches under variable developmental phases from juvenile, adult vegetative to reproductive phases. To understand the molecular mechanism behind diverse phases of axillary branches, we analyzed the levels of primary miR156 expressions because miR156-SPL module is a key regulator for developmental phase transition. We found that in Pajares, miR156 levels were highly variable among the axillary branches, which causes differential sensitivity to vernalization. Thus, the axillary branches expressing high levels of miR156 remain in juvenile phase even after vernalization, whereas the axillary branches expressing low levels of miR156 produce flowers after vernalization. In contrast, every axillary branches of Arabidopsis winter annual Sy-0 expressed similar levels of miR156 and synchronously responded to vernalization, which causes holistic flowering. Therefore, we suggest that variable miR156 expression levels and the resulting differential response to vernalization among axillary branches are distinctive features determining polycarpic perenniality of A. alpina Pajares.

Characterization of microRNAs of Beta macrocarpa and their responses to Beet necrotic yellow vein virus infection

Plant microRNAs (miRNAs) are a class of non-coding RNAs that play important roles in plant development, defense, and symptom development. Here, 547 known miRNAs representing 129 miRNA families, and 282 potential novel miRNAs were identified in Beta macrocarpa using small RNA deep sequencing. A phylogenetic analysis was performed, and 8 Beta lineage-specific miRNAs were identified. Through a differential expression analysis, miRNAs associated with Beet necrotic yellow vein virus (BNYVV) infection were identified and confirmed using a microarray analysis and stem-loop RT-qPCR. In total, 103 known miRNAs representing 38 miRNA families, and 45 potential novel miRNAs were differentially regulated, with at least a two-fold change, in BNYVV-infected plants compared with that of the mock-inoculated control. Targets of these differentially expressed miRNAs were also predicted by degradome sequencing. These differentially expressed miRNAs were involved in hormone biosynthesis and signal transduction pathways, and enhanced axillary bud development and plant defenses. This work is the first to describe miRNAs of the plant genus Beta and may offer a reference for miRNA research in other species in the genus. It provides valuable information on the pathogenicity mechanisms of BNYVV.

Control of chrysanthemum flowering through integration with an aging pathway

Age, as a threshold of floral competence acquisition, prevents precocious flowering when there is insufficient biomass, and ensures flowering independent of environmental conditions; however, the underlying regulatory mechanisms are largely unknown. In this study, silencing the expression of a nuclear factor gene, CmNF-YB8, from the short day plant chrysanthemum (Chrysanthemum morifolium), results in precocious transition from juvenile to adult, as well as early flowering, regardless of day length conditions. The expression of SQUAMOSA PROMOTER BINDING-LIKE (SPL) family members, SPL3, SPL5, and SPL9, is upregulated in CmNF-YB8-RNAi plants, while expression of the microRNA, cmo-MIR156, is downregulated. In addition, CmNF-YB8 is shown to bind to the promoter of the cmo-MIR156 gene. Ectopic expression of cmo-miR156, using a virus-based microRNA expression system, restores the early flowering phenotype caused by CmNF-YB8 silencing. These results show that CmNF-YB8 influences flowering time through directly regulating the expression of cmo-MIR156 in the aging pathway.

The mechanisms by which plant age regulates flowering remain incompletely understood. Here the authors show that age dependent regulation of SPL transcription factors by miR156 influence flowering via control of NF-YB8 expression in Chrysanthemum.

Genomic Survey, Characterization, and Expression Profile Analysis of the SBP Genes in Pineapple (Ananas comosus L.)

Gene expression is regulated by transcription factors, which play many significant developmental processes. SQUAMOSA promoter-binding proteins (SBP) perform a variety of regulatory functions in leaf, flower, and fruit development, plant architecture, and sporogenesis. 16 SBP genes were identified in pineapple and were divided into four groups on basis of phylogenetic analysis. Five paralogs in pineapple for SBP genes were identified with Ka/Ks ratio varied from 0.20 for AcSBP14 and AcSBP15 to 0.36 for AcSBP6 and AcSBP16, respectively. 16 SBP genes were located on 12 chromosomes out of 25 pineapple chromosomes with highly conserved protein sequence structures. The isoionic points of SBP ranged from 6.05 to 9.57, while molecular weight varied from 22.7 to 121.9 kD. Expression profiles of SBP genes revealed that AcSBP7 and AcSBP15 (leaf), AcSBP13, AcSBP12, AcSBP8, AcSBP16, AcSBP9, and AcSBP11 (sepal), AcSBP6, AcSBP4, and AcSBP10 (stamen), AcSBP14, AcSBP1, and AcSBP5 (fruit) while the rest of genes showed low expression in studied tissues. Four genes, that is, AcSBP11, AcSBP6, AcSBP4, and AcSBP12, were highly expressed at 4°C, while AcSBP16 were upregulated at 45°C. RNA-Seq was validated through qRT-PCR for some genes. Salt stress-induced expression of two genes, that is, AcSBP7 and AcSBP14, while in drought stress, AcSBP12 and AcSBP15 were highly expressed. Our study lays a foundation for further gene function and expression studies of SBP genes in pineapple.

Genome-wide identification of leaf abscission associated microRNAs in sugarcane (Saccharum officinarum L.)

Background

Sugarcane (Saccharum officinarum L.) is an economically important crop, mainly due to the production of sugar and biofuel (Azevedo RA, Carvalho RF, Cia MC, & Gratão PL, Trop Plant Biol 4:42-51, 2011). Grown mainly in tropical and subtropical countries, sugarcane is a highly polyploid plant with up to ten copies of each chromosome, which increases the difficulties of genome assembly and genetic, physiologic and biochemical analyses. The increasing demands of sugar and the increasing cost of sugarcane harvest require sugarcane varieties which can shed their leaves during the maturity time, so it is important to study the mechanism of leaf abscission in sugarcane.

Results

To improve the understanding of miRNA roles in sugarcane leaf abscission, we reported the genome-wide characterization of miRNAs and their putative targets in sugarcane using deep sequencing for six small RNA libraries. In total, 93 conserved miRNAs and 454 novel miRNAs were identified in sugarcane using previously reported transcriptome as reference. Among them, 25 up-regulated and 13 down-regulated miRNAs were identified in leaf abscission sugarcane plants (LASP) compared to leaf packaging sugarcane plants (LPSP). Target prediction revealed several miRNA-mRNA modules including miR156-SPL, miR319-TPR2, miR396-GRF and miR408-LAC3 might be involved in the sugarcane leaf abscission. KEGG pathway enrichment analysis showed differentially expressed miRNAs may regulate pathways like “plant hormone signal transduction” and “plant-pathogen interaction”, which is consistent with previous transcriptome study. In addition, we identified 96 variant miRNAs with 135 single nucleotide polymorphisms (SNPs). The expression of sugarcane miRNAs and variant miRNAs were confirmed by qRT-PCR. We identified a possible poaceae specific miRNA called miR5384 for the first time in sugarcane.

Conclusions

We not only reported miR5384, a possible poaceae specific miRNA, for the first time in sugarcane but also presented some miRNA-mRNA modules including miR156-SPL, miR319-TPR2, miR396-GRF and miR408-LAC in sugarcane. These modules might be involved in the regulation of sugarcane leaf abscission during the maturity time. All of these findings may lay ground work for future application of sugarcane breeding program and benefit research studies of sugarcane miRNAs.

Electronic supplementary material

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

The use of high-throughput small RNA sequencing reveals differentially expressed microRNAs in response to aster yellows phytoplasma-infection in Vitis vinifera cv. 'Chardonnay'

Phytoplasmas are cell wall-less plant pathogenic bacteria responsible for major crop losses throughout the world. In grapevine they cause grapevine yellows, a detrimental disease associated with a variety of symptoms. The high economic impact of this disease has sparked considerable interest among researchers to understand molecular mechanisms related to pathogenesis. Increasing evidence exist that a class of small non-coding endogenous RNAs, known as microRNAs (miRNAs), play an important role in post-transcriptional gene regulation during plant development and responses to biotic and abiotic stresses. Thus, we aimed to dissect complex high-throughput small RNA sequencing data for the genome-wide identification of known and novel differentially expressed miRNAs, using read libraries constructed from healthy and phytoplasma-infected Chardonnay leaf material. Furthermore, we utilised computational resources to predict putative miRNA targets to explore the involvement of possible pathogen response pathways. We identified multiple known miRNA sequence variants (isomiRs), likely generated through post-transcriptional modifications. Sequences of 13 known, canonical miRNAs were shown to be differentially expressed. A total of 175 novel miRNA precursor sequences, each derived from a unique genomic location, were predicted, of which 23 were differentially expressed. A homology search revealed that some of these novel miRNAs shared high sequence similarity with conserved miRNAs from other plant species, as well as known grapevine miRNAs. The relative expression of randomly selected known and novel miRNAs was determined with real-time RT-qPCR analysis, thereby validating the trend of expression seen in the normalised small RNA sequencing read count data. Among the putative miRNA targets, we identified genes involved in plant morphology, hormone signalling, nutrient homeostasis, as well as plant stress. Our results may assist in understanding the role that miRNA pathways play during plant pathogenesis, and may be crucial in understanding disease symptom development in aster yellows phytoplasma-infected grapevines.

Genome-Wide Organization and Expression Profiling of the SBP-Box Gene Family in Chinese Jujube (Ziziphus jujuba Mill.)

Transcription factors play vital roles in the developmental processes of plants. The SQUAMOSA promoter binding protein (SBP) genes encode a family of plant-specific transcription factors and plays many crucial roles in plant development. In this study, 16 SBP-box gene family members were identified in Ziziphus jujuba Mill. Dongzao (Dongzao), which were distributed over 8 chromosomes. They were classified into seven groups according to their phylogenetic relationships with other SBP-box gene families. Within each group, genes shared similar exon-intron structures and motif locations. The number of exons varied among the groups. We identified 12 homologous gene pairs between Dongzao and Arabidopsis. Expression profiling revealed that ZjSBP02 and ZjSBP14 expressed highly in mature fruits, ZjSBP01 expressed higher in mature leaves than other tissues and the expression level of ZjSBP12 was much higher in the flowers. The transcriptome analysis indicated that ZjSBPs had different expression patterns in various tissues. This study represents the first systematic analysis of the SBP-box gene family in Z. jujuba. The data presented here provides a foundation for understanding the crucial roles of ZjSBP genes in plant growth and development.

IPA1 functions as a downstream transcription factor repressed by D53 in strigolactone signaling in rice

Strigolactones (SLs), a group of carotenoid derived terpenoid lactones, are root-to-shoot phytohormones suppressing shoot branching by inhibiting the outgrowth of axillary buds. DWARF 53 (D53), the key repressor of the SL signaling pathway, is speculated to regulate the downstream transcriptional network of the SL response. However, no downstream transcription factor targeted by D53 has yet been reported. Here we report that Ideal Plant Architecture 1 (IPA1), a key regulator of the plant architecture in rice, functions as a direct downstream component of D53 in regulating tiller number and SL-induced gene expression. We showed that D53 interacts with IPA1 in vivo and in vitro and suppresses the transcriptional activation activity of IPA1. We further showed that IPA1 could directly bind to the D53 promoter and plays a critical role in the feedback regulation of SL-induced D53 expression. These findings reveal that IPA1 is likely one of the long-speculated transcription factors that act with D53 to mediate the SL-regulated tiller development in rice.

Genome-wide dissection of heterosis for yield traits in two-line hybrid rice populations

Heterosis has been widely utilized in agriculture and is important for world food safety. Many genetic models have been proposed as mechanisms underlying heterosis during the past century, yet more evidence is needed to support such models. To investigate heterosis in two-line hybrid rice, we generated a partial diallel crossing scheme, which consisted of approximately 500 F1 hybrids derived from 14 male sterile lines and 39 restorer lines. In this population, increased panicle number played the most important role in yield heterosis of hybrid rice. Genome-wide association studies identified many QTLs related to the yield traits of F1 hybrids, better paternal heterosis and special combining ability. Relevant genes, including Hd3a, qGL3, OsmiR156h, and LAX2, were identified as candidates within these QTLs. Nearly forty percent of the QTLs had only two genotypes in the F1 hybrids, mainly because the maternal lines were under intense selective pressure. Further analysis found male sterile lines and restorer lines made different superior allele contributions to F1 hybrids, and their contributions varied among different traits. These results extend our understanding of the molecular basis of heterosis in two-line hybrid rice.

Response of miR156-SPL Module during the Red Peel Coloration of Bagging-Treated Chinese Sand Pear (Pyrus pyrifolia Nakai)

MicroRNA156 is an evolutionarily highly conserved plant micro-RNA (miRNA) that controls an age-dependent flowering pathway. miR156 and its target SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes regulate anthocyanin accumulation in plants, but it is unknown whether this process is affected by light. Red Chinese sand pear (Pyrus pyrifolia) fruits exhibit a unique coloration pattern in response to bagging treatments, which makes them appropriate for studying the molecular mechanism underlying light-induced anthocyanin accumulation in fruit. Based on high-throughput miRNA and degradome sequencing data, we determined that miR156 was expressed in pear fruit peels, and targeted four SPL genes. Light-responsive elements were detected in the promoter regions of the miR156a and miR156ba precursors. We identified 19 SPL genes using the “Suli” pear (Pyrus pyrifolia Chinese White Pear Group) genome database, of which seven members were putative miR156 targets. The upregulated expression of anthocyanin biosynthetic and regulatory genes and downregulated expression of PpSPL2, PpSPL5, PpSPL7, PpSPL9, PpSPL10, PpSPL13, PpSPL16, PpSPL17, and PpSPL18 were observed in pear fruits after bags were removed from plants during the anthocyanin accumulation period. Additionally, miR156a/ba/g/s/sa abundance increased after bags were removed. Yeast two-hybrid results suggested that PpMYB10, PpbHLH, and PpWD40 could form a protein complex, probably involved in anthocyanin biosynthesis. Additionally, PpSPL10 and PpSPL13 interacted with PpMYB10. The results obtained in this study are helpful in understanding the possible role of miR156 and its target PpSPL genes in regulating light-induced red peel coloration and anthocyanin accumulation in pear.

The SBP-Box Gene VpSBP11 from Chinese Wild Vitis Is Involved in Floral Transition and Affects Leaf Development

Flowering occurs in angiosperms during a major developmental transition from vegetative growth to the reproductive phase. Squamosa promoter binding protein (SBP)-box genes have been found to play critical roles in regulating flower and fruit development, but their roles in grapevine have remained unclear. To better understand the functions of the grape SBP-box genes in both vegetative and reproductive growth phases, a full-length complementary DNA (cDNA) sequence of the putative SBP-box transcription factor gene, VpSBP11, was obtained from Chinese wild grapevine Vitis pseudoreticulata Wen Tsai Wang (W. T. Wang) clone ‘Baihe-35-1’. VpSBP11 encoded a putative polypeptide of 170 amino acids with a highly conserved SBP-domain with two zinc-binding sites of the Cx2C-x3-H-x11-C-x6-H (C2HCH) type and a nuclear localization signal. We confirmed that the VpSBP11 protein was targeted to the nucleus and possessed transcriptional activation activity by subcellular localization and trans-activation assay. Over-expression of VpSBP11 in Arabidopsis thaliana was shown to activate the FUL gene, and subsequently the AP1 and LFY genes, all of which were floral meristem identity genes, and to cause earlier flowering than in wild type (WT) plants. The pattern of vegetative growth was also different between the transgenic and WT plants. For example, in the VpSBP11 over-expressing transgenic plants, the number of rosette leaves was less than that of WT; the petiole was significantly elongated; and the rosette and cauline leaves curled upwards or downwards. These results were consistent with VpSBP11 acting as a transcription factor during the transition from the vegetative stage to the reproductive stage.

Crop-model assisted phenomics and genome-wide association study for climate adaptation of indica rice. 2. Thermal stress and spikelet sterility

Low night and high day temperatures during sensitive reproductive stages cause spikelet sterility in rice. Phenotyping of tolerance traits in the field is difficult because of temporal interactions with phenology and organ temperature differing from ambient. Physiological models can be used to separate these effects. A 203-accession indica rice diversity panel was phenotyped for sterility in ten environments in Senegal and Madagascar and climate data were recorded. Here we report on sterility responses while a companion study reported on phenology. The objectives were to improve the RIDEV model of rice thermal sterility, to estimate response traits by fitting model parameters, and to link the response traits to genomic regions through genome-wide association studies (GWAS). RIDEV captured 64% of variation of sterility when cold acclimation during vegetative stage was simulated, but only 38% when it was not. The RIDEV parameters gave more and stronger quantitative trait loci (QTLs) than index variables derived more directly from observation. The 15 QTLs identified at P<1 × 10-5 (33 at P<1 × 10-4) were related to sterility effects of heat, cold, cold acclimation, or unexplained causes (baseline sterility). Nine annotated genes were found on average within the 50% linkage disequilibrium (LD) region. Among them, one to five plausible candidate genes per QTL were identified based on known expression profiles (organ, stage, stress factors) and function. Meiosis-, development- and flowering-related genes were frequent, as well a stress signaling kinases and transcription factors. Putative epigenetic factors such as DNA methylases or histone-related genes were frequent in cold-acclimation QTLs, and positive-effect alleles were frequent in cold-tolerant highland rice from Madagascar. The results indicate that epigenetic control of acclimation may be important in indica rice genotypes adapted to cool environments.

miR156-Targeted SBP-Box Transcription Factors Interact with DWARF53 to Regulate TEOSINTE BRANCHED1 and BARREN STALK1 Expression in Bread Wheat

Genetic and environmental factors affect bread wheat (Triticum aestivum) plant architecture, which determines grain yield. In this study, we demonstrate that miR156 controls bread wheat plant architecture. We show that overexpression of tae-miR156 in bread wheat cultivar Kenong199 leads to increased tiller number and severe defects in spikelet formation, probably due to the tae-miR156-mediated repression of a group of SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes. Furthermore, we found that the expression of two genes TEOSINTE BRANCHED1 (TaTB1) and BARREN STALK1 (TaBA1), whose orthologous genes in diverse plant species play conserved roles in regulating plant architecture, is markedly reduced in the tae-miR156-OE bread wheat plants. Significantly, we demonstrate that the strigolactone (SL) signaling repressor DWARF53 (TaD53), which physically associates with the transcriptional corepressor TOPLESS, can directly interact with the N-terminal domains of miR156-controlled TaSPL3/17. Most importantly, TaSPL3/17-mediated transcriptional activation of TaBA1 and TaTB1 can be largely repressed by TaD53 in the transient expression system. Our results reveal potential association between miR156-TaSPLs and SL signaling pathways during bread wheat tillering and spikelet development.

MiR529a modulates panicle architecture through regulating SQUAMOSA PROMOTER BINDING-LIKE genes in rice (Oryza sativa)

MiR529a affects rice panicle architecture by targeting OsSPL2,OsSPL14 and OsSPL17 genes that could regulate their downstream panicle related genes. The panicle architecture determines the grain yield and quality of rice, which could be regulated by many transcriptional factors. The SQUAMOSA PROMOTER BINDING-LIKE (SPL) transcription factors are involved in the regulation of panicle development, which are targeted by miR156 and miR529. The expression profile demonstrated that miR529a is preferentially expressed in the early panicle of rice and it might regulate panicle development in rice. However, the regulation mechanism of miR529-SPL is still not clear. In this study, we predicted five miR529a putative target genes, OsSPL2, OsSPL14, OsSPL16, OsSPL17 and OsSPL18, while only the expression of OsSPL2, OsSPL14, and OsSPL17 was regulated by miR529a in the rice panicle. Overexpression of miR529a dramatically affected panicle architecture, which was regulated by OsSPL2, OsSPL14, and OsSPL17. Furthermore, the 117, 35, and 25 pathway genes associated with OsSPL2, OsSPL14 and OsSPL17, respectively, were predicted, and they shared 20 putative pathway genes. Our results revealed that miR529a could play a vital role in the regulation of panicle architecture through regulating OsSPL2, OsSPL14, OsSPL17 and the complex networks formed by their pathway and downstream genes. These findings will provide new genetic resources for reshaping ideal plant architecture and breeding high yield rice varieties.

Transcriptome Profiling of Wheat Inflorescence Development from Spikelet Initiation to Floral Patterning Identified Stage-Specific Regulatory Genes

Early reproductive development in cereals is crucial for final grain number per spike and hence the yield potential of the crop. To date, however, no systematic analyses of gene expression profiles during this important process have been conducted for common wheat (Triticum aestivum). Here, we studied the transcriptome profiles at four stages of early wheat reproductive development, from spikelet initiation to floral organ differentiation. K-means clustering and stage-specific transcript identification detected dynamically expressed homeologs of important transcription regulators in spikelet and floral meristems that may be involved in spikelet initiation, floret meristem specification, and floral organ patterning, as inferred from their homologs in model plants. Small RNA transcriptome sequencing discovered key microRNAs that were differentially expressed during wheat inflorescence development alongside their target genes, suggesting that miRNA-mediated regulatory mechanisms for floral development may be conserved in cereals and Arabidopsis. Our analysis was further substantiated by the functional characterization of the ARGONAUTE1d (AGO1d) gene, which was initially expressed in stamen primordia and later in the tapetum during anther maturation. In agreement with its stage-specific expression pattern, the loss of function of the predominantly expressed B homeolog of AGO1d in a tetraploid durum wheat mutant resulted in smaller anthers with more infertile pollens than the wild type and a reduced grain number per spike. Together, our work provides a first glimpse of the gene regulatory networks in wheat inflorescence development that may be pivotal for floral and grain development, highlighting potential targets for genetic manipulation to improve future wheat yields.

Genome-wide identification and expression analysis of SBP-like transcription factor genes in Moso Bamboo (Phyllostachys edulis)

BACKGROUND

The SQUAMOSA promoter binding protein-like (SPL) proteins are plant-specific transcription factors (TFs) that function in a variety of developmental processes including growth, flower development, and signal transduction. SPL proteins are encoded by a gene family, and these genes have been characterized in two model grass species, Zea mays and Oryza sativa. The SPL gene family has not been well studied in moso bamboo (Phyllostachys edulis), a woody grass species.

RESULTS

We identified 32 putative PeSPL genes in the P. edulis genome. Phylogenetic analysis arranged the PeSPL protein sequences in eight groups. Similarly, phylogenetic analysis of the SBP-like and SBP proteins from rice and maize clustered them into eight groups analogous to those from P. edulis. Furthermore, the deduced PeSPL proteins in each group contained very similar conserved sequence motifs. Our analyses indicate that the PeSPL genes experienced a large-scale duplication event ~15 million years ago (MYA), and that divergence between the PeSPL and OsSPL genes occurred 34 MYA. The stress-response expression profiles and tissue-specificity of the putative PeSPL gene promoter regions showed that SPL genes in moso bamboo have potential biological functions in stress resistance as well as in growth and development. We therefore examined PeSPL gene expression in response to different plant hormone and drought (polyethylene glycol-6000; PEG) treatments to mimic biotic and abiotic stresses. Expression of three (PeSPL10, -12, -17), six (PeSPL1, -10, -12, -17, -20, -31), and nine (PeSPL5, -8, -9, -14, -15, -19, -20, -31, -32) genes remained relatively stable after treating with salicylic acid (SA), gibberellic acid (GA), and PEG, respectively, while the expression patterns of other genes changed. In addition, analysis of tissue-specific expression of the moso bamboo SPL genes during development showed differences in their spatiotemporal expression patterns, and many were expressed at high levels in flowers and leaves.

CONCLUSIONS

The PeSPL genes play important roles in plant growth and development, including responses to stresses, and most of the genes are expressed in different tissues. Our study provides a comprehensive understanding of the PeSPL gene family and may enable future studies on the function and evolution of SPL genes in moso bamboo.

Genome-wide identification and characterization of the SPL gene family in Ziziphus jujuba

SQUAMOSA Promoter-Binding Protein-Likes (SPLs) are plant specific transcription factors playing important roles in plant growth and development. The SPL gene family has been studied in various plant species; however, there is no report about SPLs in Zizyphus jujuba. In this study, we identified 18 putative ZjSPL genes in Z. jujuba using a genome-wide analysis. Sequence features, gene structures, conserved domains and motifs were analyzed. The phylogenetic relationships of SPLs in Z. jujuba and A. thaliana were revealed. A total of 5 pairs of ZjSPLs were identified, suggesting the importance of gene duplication in SPL gene expansion in Z. jujuba. In addition, 11 of the 18 ZjSPLs, belonging to G1, G2 and G5 subgroups, were found to be targets of miR156, suggesting the conservation of miR156-mediated posttranscriptional regulation in plants. Expression analysis revealed that eight ZjSPL genes were responsive to the infection of witches'-broom phytoplasma. Our results provide a basis for the further elucidation of the biological function of ZjSPLs and their regulation in witches'-broom disease.

Identification and positional distribution analysis of transcription factor binding sites for genes from the wheat fl-cDNA sequences

The binding sites of transcription factors (TFs) in upstream DNA regions are called transcription factor binding sites (TFBSs). TFBSs are important elements for regulating gene expression. To date, there have been few studies on the profiles of TFBSs in plants. In total, 4,873 sequences with 5' upstream regions from 8530 wheat fl-cDNA sequences were used to predict TFBSs. We found 4572 TFBSs for the MADS TF family, which was twice as many as for bHLH (1951), B3 (1951), HB superfamily (1914), ERF (1820), and AP2/ERF (1725) TFs, and was approximately four times higher than the remaining TFBS types. The percentage of TFBSs and TF members showed a distinct distribution in different tissues. Overall, the distribution of TFBSs in the upstream regions of wheat fl-cDNA sequences had significant difference. Meanwhile, high frequencies of some types of TFBSs were found in specific regions in the upstream sequences. Both TFs and fl-cDNA with TFBSs predicted in the same tissues exhibited specific distribution preferences for regulating gene expression. The tissue-specific analysis of TFs and fl-cDNA with TFBSs provides useful information for functional research, and can be used to identify relationships between tissue-specific TFs and fl-cDNA with TFBSs. Moreover, the positional distribution of TFBSs indicates that some types of wheat TFBS have different positional distribution preferences in the upstream regions of genes.

Global gene expression profiling related to temperature-sensitive growth abnormalities in interspecific crosses between tetraploid wheat and Aegilops tauschii

Triploid wheat hybrids between tetraploid wheat and Aegilops tauschii sometimes show abnormal growth phenotypes, and the growth abnormalities inhibit generation of wheat synthetic hexaploids. In type II necrosis, one of the growth abnormalities, necrotic cell death accompanied by marked growth repression occurs only under low temperature conditions. At normal temperature, the type II necrosis lines show grass-clump dwarfism with no necrotic symptoms, excess tillers, severe dwarfism and delayed flowering. Here, we report comparative expression analyses to elucidate the molecular mechanisms of the temperature-dependent phenotypic plasticity in the triploid wheat hybrids. We compared gene and small RNA expression profiles in crown tissues to characterize the temperature-dependent phenotypic plasticity. No up-regulation of defense-related genes was observed under the normal temperature, and down-regulation of wheat APETALA1-like MADS-box genes, considered to act as flowering promoters, was found in the grass-clump dwarf lines. Some microRNAs, including miR156, were up-regulated, whereas the levels of transcripts of the miR156 target genes SPLs, known to inhibit tiller and branch number, were reduced in crown tissues of the grass-clump dwarf lines at the normal temperature. Unusual expression of the miR156/SPLs module could explain the grass-clump dwarf phenotype. Dramatic alteration of gene expression profiles, including miRNA levels, in crown tissues is associated with the temperature-dependent phenotypic plasticity in type II necrosis/grass-clump dwarf wheat hybrids.

Virus-Based MicroRNA Silencing and Overexpressing in Common Wheat (Triticum aestivum L.)

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that arise from large RNA precursors with a stem-loop structure and play important roles in plant development and responses to environmental stress. Although a hundred and nineteen wheat miRNAs have been identified and registered in the miRBase (Release 21.0, June, 2014; http://www.mirbase.org), the functional characterization of these miRNAs in wheat growth and development is lagging due to lack of effective techniques to investigate endogenous miRNA functions in wheat. Here we report barley stripe mosaic virus(BSMV)-based miRNA overexpression and silence systems that can be applied to study miRNA functions in wheat. By utilizing the BSMV system, we successfully knocked down endogenous miR156 and miR166 levels and over-expressed endogenous miR156 and artificial miRNA against phytoene desaturase gene PDS (amiR-PDS) in wheat. amiR-PDS expression caused a great reduction in endogenous mRNA abundance of PDS gene in wheat plant, leading to leaf obviously photobleaching. miR156 silencing led to a great increase in the mRNA level of its target gene SPL2, resulting in a leaf-curl phenotype in wheat seedlings. In contrast, overexpression of miR156 led to a significant reduction in the mRNA level of SPL2 in wheat seedlings, resulting in serious delay of the vegetative phase transitions as well as booting and flowering in wheat. These confirmed that miR156 regulates wheat development and booting time through SPL genes. In summary, the BSMV-based miRNA overexpression and silence systems have extraordinary potential not only for functional study of protein-encoding genes but also for miRNA genes in wheat.

Molecular mechanisms of flowering under long days and stem growth habit in soybean

Precise timing of flowering is critical to crop adaptation and productivity in a given environment. A number of classical E genes controlling flowering time and maturity have been identified in soybean [Glycine max (L.) Merr.]. The public availability of the soybean genome sequence has accelerated the identification of orthologues of Arabidopsis flowering genes and their functional analysis, and has allowed notable progress towards understanding the molecular mechanisms of flowering in soybean. Great progress has been made particularly in identifying genes and modules that inhibit flowering in long-day conditions, because a reduced or absent inhibition of flowering by long daylengths is an essential trait for soybean, a short-day (SD) plant, to expand its adaptability toward higher latitude environments. In contrast, the molecular mechanism of floral induction by SDs remains elusive in soybean. Here we present an update on recent work on molecular mechanisms of flowering under long days and of stem growth habit, outlining the progress in the identification of genes responsible, the interplay between photoperiod and age-dependent miRNA-mediated modules, and the conservation and divergence in photoperiodic flowering and stem growth habit in soybean relative to other legumes, Arabidopsis, and rice (Oryza sativa L.).

Genomic organization, phylogenetic comparison, and expression profiles of the SPL family genes and their regulation in soybean

SQUAMOSA Promoter-Binding Protein-Like (SPL) genes form a major family of plant-specific transcription factors and play an important role in plant growth and development. In this study, we report the identification of 41 SPL genes (GmSPLs) in the soybean genome. Phylogenetic analysis revealed that these genes were divided into five groups (groups 1-5). Further, exon/intron structure and motif composition revealed that the GmSPL genes are conserved within their same group. The N-terminal zinc finger 1 (Zn1) of the SBP domain was a CCCH (Cys3His1) and the C terminus zinc finger 2 (Zn2) was a CCHC (Cys2HisCys) type. The 41 GmSPL genes were distributed unevenly on 17 of the 20 chromosomes, with tandem and segmental duplication events. We found that segmental duplication has made an important contribution to soybean SPL gene family expansion. The Ka/Ks ratios revealed that the duplicated GmSPL genes evolved under the effect of purifying selection. In addition, 17 of the 41 GmSPLs were found as targets of miR156; these might be involved in their posttranscriptional regulation through miR156. Importantly, RLM-RACE analysis confirmed the GmmiR156-mediated cleavage of GmSPL2a transcript in 2-4 mm stage of soybean seed. Alternative splicing events in 9 GmSPLs were detected which produces transcripts and proteins of different lengths that may modulate protein signaling, binding, localization, stability, and other properties. Expression analysis of the soybean SPL genes in various tissues and different developmental stages of seed suggested distinct spatiotemporal patterns. Differences in the expression patterns of miR156-targeted and miR156-non-targeted soybean SPL genes suggest that miR156 plays key functions in soybean development. Our results provide an important foundation for further uncovering the crucial roles of GmSPLs in the development of soybean and other biological processes.

Heteroblastic Development of Transfer Cells Is Controlled by the microRNA miR156/SPL Module

We report that wall ingrowth deposition in phloem parenchyma (PP) transfer cells (TCs) in leaf veins of Arabidopsis (Arabidopsis thaliana) represents a novel trait of heteroblasty. Development of PP TCs involves extensive deposition of wall ingrowths adjacent to cells of the sieve element/companion cell complex. These PP TCs potentially facilitate phloem loading by enhancing efflux of symplasmic Suc for subsequent active uptake into cells of the sieve element/companion cell complex. PP TCs with extensive wall ingrowths are ubiquitous in mature cotyledons and juvenile leaves, but dramatically less so in mature adult leaves, an observation consistent with PP TC development reflecting vegetative phase change (VPC) in Arabidopsis. Consistent with this conclusion, the abundance of PP TCs with extensive wall ingrowths varied across rosette development in three ecotypes displaying differing durations of juvenile phase, and extensive deposition of wall ingrowths was observed in rejuvenated leaves following prolonged defoliation. PP TC development across juvenile, transition, and adult leaves correlated positively with levels of miR156, a major regulator of VPC in plants, and corresponding changes in wall ingrowth deposition were observed when miR156 was overexpressed or its activity suppressed by target mimicry. Analysis of plants carrying miR156-resistant forms of SQUAMOSA PROMOTER BINDING PROTEIN LIKE (SPL) genes showed that wall ingrowth deposition was increased in SPL9-group but not SPL3-group genes, indicating that SPL9-group genes may function as negative regulators of wall ingrowth deposition in PP TCs. Collectively, our results point to wall ingrowth deposition in PP TCs being under control of the genetic program regulating VPC.

Combined analysis of mRNA and miRNA identifies dehydration and salinity responsive key molecular players in citrus roots

Citrus is one of the most economically important fruit crops around world. Drought and salinity stresses adversely affected its productivity and fruit quality. However, the genetic regulatory networks and signaling pathways involved in drought and salinity remain to be elucidated. With RNA-seq and sRNA-seq, an integrative analysis of miRNA and mRNA expression profiling and their regulatory networks were conducted using citrus roots subjected to dehydration and salt treatment. Differentially expressed (DE) mRNA and miRNA profiles were obtained according to fold change analysis and the relationships between miRNAs and target mRNAs were found to be coherent and incoherent in the regulatory networks. GO enrichment analysis revealed that some crucial biological processes related to signal transduction (e.g. ‘MAPK cascade’), hormone-mediated signaling pathways (e.g. abscisic acid- activated signaling pathway’), reactive oxygen species (ROS) metabolic process (e.g. ‘hydrogen peroxide catabolic process’) and transcription factors (e.g., ‘MYB, ZFP and bZIP’) were involved in dehydration and/or salt treatment. The molecular players in response to dehydration and salt treatment were partially overlapping. Quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) analysis further confirmed the results from RNA-seq and sRNA-seq analysis. This study provides new insights into the molecular mechanisms how citrus roots respond to dehydration and salt treatment.

Azacytidine and miR156 promote rooting in adult but not in juvenile Arabidopsis tissues

Poor adventitious root (AR) formation is a major obstacle in micropropagation and conventional vegetative propagation of many crops. It is affected by many endogenous and exogenous factors. With respect to endogenous factors, the phase change from juvenile to adult has a major influence on AR formation and rooting is usually much reduced or even fully inhibited in adult tissues. It has been reported that the phase change is characterized by an increase in DNA-methylation and a decrease in the expression of microRNA156 (miR156). In this paper, we examined the effect of azacytidine (AzaC) and miR156 on AR formation in adult and juvenile Arabidopsis tissues. To identify the ontogenetic state researchers have used flowering or leaf morphology. We have used the rootability which allows - in contrast with both other characteristics- to examine the ontogenetic state at the cellular level. Overexpression of miR156 promoted only the rooting of adult tissues indicating that the phase change-associated loss in tissues' competence to develop ARs is also under the control of miR156. Azacytidine inhibits DNA methylation during DNA replication. Azacytidine treatment also promoted AR formation in nonjuvenile tissues but had no or little effect in juvenile tissues. Its addition during seedling growth (by which all tissues become hypomethylated) or during the rooting treatment (by which only those cells become hypomethylated that are generated after taking the explant) are both effective in the promotion of rooting. An AzaC treatment may be useful in tissue culture for crops that are recalcitrant to root.

Field-grown miR156 transgenic switchgrass reproduction, yield, global gene expression analysis, and bioconfinement

BACKGROUND

Genetic engineering has been effective in altering cell walls for biofuel production in the bioenergy crop, switchgrass (Panicum virgatum). However, regulatory issues arising from gene flow may prevent commercialization of engineered switchgrass in the eastern United States where the species is native. Depending on its expression level, microRNA156 (miR156) can reduce, delay, or eliminate flowering, which may serve to decrease transgene flow. In this unique field study of transgenic switchgrass that was permitted to flower, two low (T14 and T35) and two medium (T27 and T37) miR156-overexpressing 'Alamo' lines with the transgene under the control of the constitutive maize (Zea mays) ubiquitin 1 promoter, along with nontransgenic control plants, were grown in eastern Tennessee over two seasons.

RESULTS

miR156 expression was positively associated with decreased and delayed flowering in switchgrass. Line T27 did not flower during the 2-year study. Line T37 did flower, but not all plants produced panicles. Flowering was delayed in T37, resulting in 70.6% fewer flowers than controls during the second field year with commensurate decreased seed yield: 1205 seeds per plant vs. 18,539 produced by each control. These results are notable given that line T37 produced equivalent vegetative aboveground biomass to the controls. miR156 transcript abundance of field-grown plants was congruent with greenhouse results. The five miR156 SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) target genes had suppressed expression in one or more of the transgenic lines. Line T27, which had the highest miR156 overexpression, showed significant downregulation for all five SPL genes. On the contrary, line T35 had the lowest miR156 overexpression and had no significant change in any of the five SPL genes.

CONCLUSIONS

Because of the research field's geographical features, this study was the first instance of any genetically engineered trait in switchgrass, in which experimental plants were allowed to flower in the field in the eastern U.S.; USDA-APHIS-BRS regulators allowed open flowering. We found that medium overexpression of miR156, e.g., line T37, resulted in delayed and reduced flowering accompanied by high biomass production. We propose that induced miR156 expression could be further developed as a transgenic switchgrass bioconfinement tool to enable eventual commercialization.

MicroRNA and Transcription Factor: Key Players in Plant Regulatory Network

Recent achievements in plant microRNA (miRNA), a large class of small and non-coding RNAs, are very exciting. A wide array of techniques involving forward genetic, molecular cloning, bioinformatic analysis, and the latest technology, deep sequencing have greatly advanced miRNA discovery. A tiny miRNA sequence has the ability to target single/multiple mRNA targets. Most of the miRNA targets are transcription factors (TFs) which have paramount importance in regulating the plant growth and development. Various families of TFs, which have regulated a range of regulatory networks, may assist plants to grow under normal and stress environmental conditions. This present review focuses on the regulatory relationships between miRNAs and different families of TFs like; NF-Y, MYB, AP2, TCP, WRKY, NAC, GRF, and SPL. For instance NF-Y play important role during drought tolerance and flower development, MYB are involved in signal transduction and biosynthesis of secondary metabolites, AP2 regulate the floral development and nodule formation, TCP direct leaf development and growth hormones signaling. WRKY have known roles in multiple stress tolerances, NAC regulate lateral root formation, GRF are involved in root growth, flower, and seed development, and SPL regulate plant transition from juvenile to adult. We also studied the relation between miRNAs and TFs by consolidating the research findings from different plant species which will help plant scientists in understanding the mechanism of action and interaction between these regulators in the plant growth and development under normal and stress environmental conditions.

miR156/SPL10 Modulates Lateral Root Development, Branching and Leaf Morphology in Arabidopsis by Silencing AGAMOUS-LIKE 79

The developmental functions of miR156-SPL regulatory network have been extensively studied in Arabidopsis, but the downstream genes regulated by each SPL have not been well characterized. In this study, Next Generation Sequencing-based transcriptome analysis was performed on roots of wild type (WT) and miR156 overexpression (miR156OE) plants. One of the SPL genes, SPL10, which represses lateral root growth in Arabidopsis, was significantly downregulated in miR156OE plants. A transcription factor, AGAMOUS-like MADS box protein 79 (AGL79), was also significantly downregulated in the miR156OE plants, but was upregulated in the SPL10 overexpression (SPL10OE) Arabidopsis plants. In addition, SPL10 was found to bind to the core consensus SPL binding sequences in AGL79 gene. Moreover, analyses of complementation lines revealed a linear relationship between SPL10 and AGL79 in regulating Arabidopsis plant development. In addition, it was observed that plant phenotypes are AGL79 dose-dependent, with higher expression causing narrow leaf shape, less number of leaves and early flowering time, whereas relatively lower AGL79 overexpression produce plants with more rosette leaves and more lateral branches. Our findings revealed direct binding of SPL10 to AGL79 promoter, which further suggests a role for miR156/SPL10 module in plant lateral root growth by directly regulating AGL79.

Variation of DNA methylation patterns associated with gene expression in rice (Oryza sativa) exposed to cadmium

We report genome-wide single-base resolution maps of methylated cytosines and transcriptome change in Cd-exposed rice. Widespread differences were identified in CG and non-CG methylation marks between Cd-exposed and Cd-free rice genomes. There are 2320 non-redundant differentially methylated regions detected in the genome. RNA sequencing revealed 2092 DNA methylation-modified genes differentially expressed under Cd exposure. More genes were found hypermethylated than those hypomethylated in CG, CHH and CHG (where H is A, C or T) contexts in upstream, gene body and downstream regions. Many of the genes were involved in stress response, metal transport and transcription factors. Most of the DNA methylation-modified genes were transcriptionally altered under Cd stress. A subset of loss of function mutants defective in DNA methylation and histone modification activities was used to identify transcript abundance of selected genes. Compared with wide type, mutation of MET1 and DRM2 resulted in general lower transcript levels of the genes under Cd stress. Transcripts of OsIRO2, OsPR1b and Os09g02214 in drm2 were significantly reduced. A commonly used DNA methylation inhibitor 5-azacytidine was employed to investigate whether DNA demethylation affected physiological consequences. 5-azacytidine provision decreased general DNA methylation levels of selected genes, but promoted growth of rice seedlings and Cd accumulation in rice plant.

Dynamic changes of small RNAs in rice spikelet development reveal specialized reproductive phasiRNA pathways

Dissection of the genetic pathways and mechanisms by which anther development occurs in grasses is crucial for both a basic understanding of plant development and for examining traits of agronomic importance such as male sterility. In rice, MULTIPLE SPOROCYTES1 (MSP1), a leucine-rich-repeat receptor kinase, plays an important role in anther development by limiting the number of sporocytes. OsTDL1a (a TPD1-like gene in rice) encodes a small protein that acts as a cofactor of MSP1 in the same regulatory pathway. In this study, we analyzed small RNA and mRNA changes in different stages of spikelets from wild-type rice, and from msp1 and ostdl1a mutants. Analysis of the small RNA data identified miRNAs demonstrating differential abundances. miR2275 was depleted in the two rice mutants; this miRNA is specifically enriched in anthers and functions to trigger the production of 24-nt phased secondary siRNAs (phasiRNAs) from PHAS loci. We observed that the 24-nt phasiRNAs as well as their precursor PHAS mRNAs were also depleted in the two mutants. An analysis of co-expression identified three Argonaute-encoding genes (OsAGO1d, OsAGO2b, and OsAGO18) that accumulate transcripts coordinately with phasiRNAs, suggesting a functional relationship. By mRNA in situ analysis, we demonstrated a strong correlation between the spatiotemporal pattern of these OsAGO transcripts and phasiRNA accumulations.