The Arabidopsis GAMYB-like genes, MYB33 and MYB65, are microRNA-regulated genes that redundantly facilitate anther development

Integration of mRNA and miRNA Analysis Reveals the Molecular Mechanism of Cotton Response to Salt Stress

To identify the regulatory network of known and novel microRNAs (miRNAs) and their targets responding to salt stress, a combined analysis of mRNA libraries, small RNA libraries, and degradome libraries were performed. In this study, we used unique molecular identifiers (UMIs), which are more sensitive, accurate, and reproducible than traditional methods of sequencing, to quantify the number of molecules and correct for amplification bias. We identified a total of 312 cotton miRNAs using seedlings at 0, 1, 3, and 6 h after NaCl treatment, including 80 known ghr-miRNAs and 232 novel miRNAs and found 155 miRNAs that displayed significant differential expression under salt stress. Among them, fifty-nine differentially expressed miRNAs were simultaneously induced in two or three tissues, while 66, 11, and 19 were specifically expressed in the roots, leaves, and stems, respectively. It is indicated there were different populations of miRNAs against salt stress in roots, leaves and stems. 399 candidate targets of salt-induced miRNAs showed significant differential expression before and after salt treatment, and 72 targets of 25 miRNAs were verified by degradome sequencing data. Furthermore, the regulatory relationship of miRNA-target gene was validated experimentally via 5'RLM-RACE, proving our data reliability. Gene ontology and KEGG pathway analysis found that salt-responsive miRNA targets among the differentially expressed genes were significantly enriched, and mainly involved in response to the stimulus process and the plant hormone signal transduction pathway. Furthermore, the expression levels of newly identified miRNA mir1 and known miRNAs miR390 and miR393 gradually decreased when subjected to continuous salt stress, while overexpression of these miRNAs both increased sensitivity to salt stress. Those newly identified miRNAs and mRNA pairs were conducive to genetic engineering and better understanding the mechanisms responding to salt stress in cotton.

Whole-Transcriptome RNA Sequencing Reveals the Global Molecular Responses and CeRNA Regulatory Network of mRNAs, lncRNAs, miRNAs and circRNAs in Response to Salt Stress in Sugar Beet (Beta vulgaris)

Sugar beet is an important sugar-yielding crop with some tolerance to salt, but the mechanistic basis of this tolerance is not known. In the present study, we have used whole-transcriptome RNA-seq and degradome sequencing in response to salt stress to uncover differentially expressed (DE) mRNAs, microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) in both leaves and roots. A competitive endogenous RNA (ceRNA) network was constructed with the predicted DE pairs, which revealed regulatory roles under salt stress. A functional analysis suggests that ceRNAs are implicated in copper redistribution, plasma membrane permeability, glycometabolism and energy metabolism, NAC transcription factor and the phosphoinositol signaling system. Overall, we conducted for the first time a full transcriptomic analysis of sugar beet under salt stress that involves a potential ceRNA network, thus providing a basis to study the potential functions of lncRNAs/circRNAs.

The Crosstalk between MicroRNAs and Gibberellin Signaling in Plants

Gibberellin (GA) is an integral phytohormone that plays prominent roles in controlling seed germination, stem elongation, leaf development and floral induction. It has been shown that GA regulates these diverse biological processes mainly through overcoming the suppressive effects of the DELLA proteins, a family of nuclear repressors of GA response. MicroRNAs (miRNAs), which have been identified as master regulators of gene expression in eukaryotes, are also involved in a wide range of plant developmental events through the repression of their target genes. The pathways of GA biosynthesis and signaling, as well as the pathways of miRNA biogenesis and regulation, have been profoundly delineated in the past several decades. Growing evidence has shown that miRNAs and GAs are coordinated in regulating plant development, as several components in GA pathways are targeted by miRNAs, and GAs also regulate the expression of miRNAs or their target genes vice versa. Here, we review the recent advances in our understanding of the molecular connections between miRNAs and GA, with an emphasis on the two miRNAs, miR156 and miR159.

AtMYB109 negatively regulates stomatal closure under osmotic stress in Arabidopsis thaliana

Osmotic stress, caused by drought, salinity, or PEG (polyethylene glycol), is one of the most important abiotic factors that hinder plant growth and development. In Arabidopsis, more than 100 R2R3-MYB transcription factors (TFs) have been identified, and many of them are involved in the transcriptional regulation of a variety of biological processes related to growth and development, as well as responses to biotic and abiotic stresses. However, the MYB TF involving in both plant development and stress response has rarely been reported. We report here that Arabidopsis AtMYB109, a R2R3-MYB TF, functions as a negative regulator of stomatal closure under osmotic stress as well as of pollen tube elongation. Under PEG-induced osmotic stress, whole leaves of AtMYB109-OXs were intensely wilted, while leaves of the wild-type (WT) and myb109 were weakly affected. Moreover, we confirmed that the wilting in AtMYB109-OXs was more severe than in WT and myb109 under drought conditions, and that after re-watering, WT and myb109 plants promptly recovered, while AtMYB109-OXs failed to survive. In addition, stomatal closure was delayed in the AtMYB109-OXs compared to the WT and myb109. However, proline content and the expression of stress-induced and proline synthesis genes were higher in the overexpression lines than in WT and myb109. Then, we observed that the expression of ICS1, a key gene in SA biosynthesis, was greatly suppressed in AtMYB109-OXs. In addition, we found that AtMYB109 expression gradually increased until the flowers were fully opened and thereafter dramatically decreased during silique development. The pollen tube growth was significantly suppressed in AtMYB109-OXs compared to the WT and myb109. Using EMSA and ChIP-qPCR, we confirmed that AtMYB109 bound to the promoter of RABA4D, a gene encoding a pollen development regulator. Taken together, we suggest the delayed stomatal closing and vulnerable phenotypes in the AtMYB109-OXs under osmotic stress are possibly directly or indirectly associated with a SA-mediated mechanism, and that AtMYB109 suppresses RABA4D that modulates pollen tube growth.

Deficiencies in the formation and regulation of anther cuticle and tryphine contribute to male sterility in cotton PGMS line

BACKGROUND

Male sterility is a simple and efficient pollination control system that is widely exploited in hybrid breeding. In upland cotton, CCRI9106, a photosensitive genetic male sterile (PGMS) mutant isolated from CCRI040029, was reported of great advantages to cotton heterosis. However, little information concerning the male sterility of CCRI9106 is known. Here, comparative transcriptome analysis of CCRI9106 (the mutant, MT) and CCRI040029 (the wild type, WT) anthers in Anyang (long-day, male sterile condition to CCRI9106) was performed to reveal the potential male sterile mechanism of CCRI9106.

RESULTS

Light and electron microscopy revealed that the male sterility phenotype of MT was mainly attributed to irregularly exine, lacking tryphine and immature anther cuticle. Based on the cytological characteristics of MT anthers, anther RNA libraries (18 in total) of tetrad (TTP), late uninucleate (lUNP) and binucleate (BNP) stages in MT and WT were constructed for transcriptomic analysis, therefore revealing a total of 870,4 differentially expressed genes (DEGs). By performing gene expression pattern analysis and protein-protein interaction (PPI) networks construction, we found down-regulation of DEGs, which enriched by the lipid biosynthetic process and the synthesis pathways of several types of secondary metabolites such as terpenoids, flavonoids and steroids, may crucial to the male sterility phenotype of MT, and resulting in the defects of anther cuticle and tryphine, even the irregularly exine. Furthermore, several lipid-related genes together with ABA-related genes and MYB transcription factors were identified as hub genes via weighted gene co-expression network analysis (WGCNA). Additionally, the ABA content of MT anthers was reduced across all stages when compared with WT anthers. At last, genes related to the formation of anther cuticle and tryphine could activated in MT under short-day condition.

CONCLUSIONS

We propose that the down-regulation of genes related to the assembly of anther cuticle and tryphine may lead to the male sterile phenotype of MT, and MYB transcription factors together with ABA played key regulatory roles in these processes. The conversion of fertility in different photoperiods may closely relate to the functional expression of these genes. These findings contribute to elucidate the mechanism of male sterility in upland cotton.

Deficiencies in the formation and regulation of anther cuticle and tryphine contribute to male sterility in cotton PGMS line

Background

Male sterility is a simple and efficient pollination control system that is widely exploited in hybrid breeding. In upland cotton, CCRI9106, a photosensitive genetic male sterile (PGMS) mutant isolated from CCRI040029, was reported of great advantages to cotton heterosis. However, little information concerning the male sterility of CCRI9106 is known. Here, comparative transcriptome analysis of CCRI9106 (the mutant, MT) and CCRI040029 (the wild type, WT) anthers in Anyang (long-day, male sterile condition to CCRI9106) was performed to reveal the potential male sterile mechanism of CCRI9106.

Results

Light and electron microscopy revealed that the male sterility phenotype of MT was mainly attributed to irregularly exine, lacking tryphine and immature anther cuticle. Based on the cytological characteristics of MT anthers, anther RNA libraries (18 in total) of tetrad (TTP), late uninucleate (lUNP) and binucleate (BNP) stages in MT and WT were constructed for transcriptomic analysis, therefore revealing a total of 870,4 differentially expressed genes (DEGs). By performing gene expression pattern analysis and protein-protein interaction (PPI) networks construction, we found down-regulation of DEGs, which enriched by the lipid biosynthetic process and the synthesis pathways of several types of secondary metabolites such as terpenoids, flavonoids and steroids, may crucial to the male sterility phenotype of MT, and resulting in the defects of anther cuticle and tryphine, even the irregularly exine. Furthermore, several lipid-related genes together with ABA-related genes and MYB transcription factors were identified as hub genes via weighted gene co-expression network analysis (WGCNA). Additionally, the ABA content of MT anthers was reduced across all stages when compared with WT anthers. At last, genes related to the formation of anther cuticle and tryphine could activated in MT under short-day condition.

Conclusions

We propose that the down-regulation of genes related to the assembly of anther cuticle and tryphine may lead to the male sterile phenotype of MT, and MYB transcription factors together with ABA played key regulatory roles in these processes. The conversion of fertility in different photoperiods may closely relate to the functional expression of these genes. These findings contribute to elucidate the mechanism of male sterility in upland cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07250-1.

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.

Cell wall Glycine-rich Protein2 is involved in tapetal differentiation and pollen maturation

The tapetum plays important roles in anther development by providing materials for pollen-wall formation and nutrients for pollen development. Here, we report the characterization of a male-sterile mutant of glycine-rich protein 2 (OsGRP2), which exhibits irregular cell division and dysfunction of the tapetum. GRP is a cellwall structural protein present in the cell walls of diverse plant species, but its function is unclear in pollen development. We found that few GRP genes are expressed in rice and thus focused on one highly expressed gene, OsGRP2. The tapetal cell walls of an OsGRP2 mutant did not thicken at the pollen mothercell stage, as a result, pollen maturation and fertility rate decreased. High OsGRP2 expression was detected in male-floral organs, and OsGRP2 was distributed in the tapetum. OsGRP2 participated in establishment of the cellwall network during early tapetum development. In conclusion, our results indicate that OsGRP2 plays important roles in the differentiation and function of the tapetum.

A study of male fertility control in Medicago truncatula uncovers an evolutionarily conserved recruitment of two tapetal bHLH subfamilies in plant sexual reproduction

Male sterility is an important tool for plant breeding and hybrid seed production. Male-sterile mutants are largely due to an abnormal development of either the sporophytic or gametophytic anther tissues. Tapetum, a key sporophytic tissue, provides nutrients for pollen development, and its delayed degeneration induces pollen abortion. Numerous bHLH proteins have been documented to participate in the degeneration of the tapetum in angiosperms, but relatively little attention has been given to the evolution of the involved developmental pathways across the phylogeny of land plants. A combination of cellular, molecular, biochemical and evolutionary analyses was used to investigate the male fertility control in Medicago truncatula. We characterized the male-sterile mutant empty anther1 (ean1) and identified EAN1 as a tapetum-specific bHLH transcription factor necessary for tapetum degeneration. Our study uncovered an evolutionarily conserved recruitment of bHLH subfamily II and III(a + c)1 in the regulation of tapetum degeneration. EAN1 belongs to the subfamily II and specifically forms heterodimers with the subfamily III(a + c)1 members, which suggests a heterodimerization mechanism conserved in angiosperms. Our work suggested that the pathway of two tapetal-bHLH subfamilies is conserved in all land plants, and likely was established before the divergence of the spore-producing land plants.

Molecular Control and Application of Male Fertility for Two-Line Hybrid Rice Breeding

The significance of the climate change may involve enhancement of plant growth as well as utilization of the environmental alterations in male fertility (MF) regulation via male sterility (MS) systems. We described that MS systems provide a fundamental platform for improvement in agriculture production and have been explicated for creating bulk germplasm of the two-line hybrids (EGMS) in rice as compared to the three-line, to gain production sustainability and exploit its immense potential. Environmental alterations such as photoperiod and/or temperature and humidity regulate MS in EGMS lines via genetic and epigenetic changes, regulation of the noncoding RNAs, and RNA-metabolism including the transcriptional factors (TFs) implication. Herein, this article enlightens a deep understanding of the molecular control of MF in EGMS lines and exploring the regulatory driving forces that function efficiently during plant adaption under a changing environment. We highlighted a possible solution in obtaining more stable hybrids through apomixis (single-line system) for seed production.

Slow Development Restores the Fertility of Photoperiod-Sensitive Male-Sterile Plant Lines

Photoperiod- and thermosensitive genic male sterility (P/TGMS) lines are widely used in crop breeding. The fertility conversion of Arabidopsis (Arabidopsis thaliana) TGMS lines including cals5-2, which is defective in callose wall formation, relies on slow development under low temperatures. In this study, we discovered that cals5-2 also exhibits PGMS. Fertility of cals5-2 was restored when pollen development was slowed under short-day photoperiods or low light intensity, suggesting that slow development restores the fertility of cals5-2 under these conditions. We found that several other TGMS lines with defects in pollen wall formation also exhibited PGMS characteristics. This similarity indicates that slow development is a general mechanism of PGMS fertility restoration. Notably, slow development also underlies the fertility recovery of TGMS lines. Further analysis revealed the pollen wall features during the formation of functional pollens of these P/TGMS lines under permissive conditions. We conclude that slow development is a general mechanism for fertility restoration of P/TGMS lines and allows these plants to take different strategies to overcome pollen formation defects.

The GAMYB-like gene SlMYB33 mediates flowering and pollen development in tomato

GAMYBs are positive GA signaling factors that exhibit essential functions in reproductive development, particularly in anther and pollen development. However, there is no direct evidence of the regulation of any GAMYB in these biological processes in tomato (Solanum lycopersicum). Here, we identified a tomato GAMYB-like gene, SlMYB33, and characterized its specific roles. SlMYB33 is predominately expressed in the stamens and pistils. During flower development, high mRNA abundance of SlMYB33 is detected in both male and female organs, such as microspore mother cells, anthers, pollen grains, and ovules. Silencing of SlMYB33 leads to delayed flowering, aberrant pollen viability, and poor fertility in tomato. Histological analyses indicate that SlMYB33 exerts its function in pollen development in the mature stage. Further transcriptomic analyses imply that the knockdown of SlMYB33 significantly inhibits the expression of genes related to flowering in shoot apices, and alters the transcription of genes controlling sugar metabolism in anthers. Taken together, our study suggests that SlMYB33 regulates tomato flowering and pollen maturity, probably by modulating the expression of genes responsible for flowering and sugar metabolism, respectively.

Global Quantitative Proteomics Studies Revealed Tissue-Preferential Expression and Phosphorylation of Regulatory Proteins in Arabidopsis

Organogenesis in plants occurs across all stages of the life cycle. Although previous studies have identified many genes as important for either vegetative or reproductive development at the RNA level, global information on translational and post-translational levels remains limited. In this study, six Arabidopsis stages/organs were analyzed using quantitative proteomics and phosphoproteomics, identifying 2187 non-redundant proteins and evidence for 1194 phosphoproteins. Compared to the expression observed in cauline leaves, the expression of 1445, 1644, and 1377 proteins showed greater than 1.5-fold alterations in stage 1–9 flowers, stage 10–12 flowers, and open flowers, respectively. Among these, 294 phosphoproteins with 472 phosphorylation sites were newly uncovered, including 275 phosphoproteins showing differential expression patterns, providing molecular markers and possible candidates for functional studies. Proteins encoded by genes preferentially expressed in anther (15), meiocyte (4), or pollen (15) were enriched in reproductive organs, and mutants of two anther-preferentially expressed proteins, acos5 and mee48, showed obviously reduced male fertility with abnormally organized pollen exine. In addition, more phosphorylated proteins were identified in reproductive stages (1149) than in the vegetative organs (995). The floral organ-preferential phosphorylation of GRP17, CDC2/CDKA.1, and ATSK11 was confirmed with western blot analysis. Moreover, phosphorylation levels of CDPK6 and MAPK6 and their interacting proteins were elevated in reproductive tissues. Overall, our study yielded extensive data on protein expression and phosphorylation at six stages/organs and provides an important resource for future studies investigating the regulatory mechanisms governing plant development.

Identification of microRNAs and their targets in inflorescences of an Ogura-type cytoplasmic male-sterile line and its maintainer fertile line of turnip (Brassica rapa ssp. rapifera) via high-throughput sequencing and degradome analysis

Cytoplasmic male sterility (CMS) is a widely used trait in angiosperms caused by perturbations in nucleus-mitochondrion interactions that suppress the production of functional pollen. MicroRNAs (miRNAs) are small non-coding RNAs that act as regulatory molecules of transcriptional or post-transcriptional gene silencing in plants. The discovery of miRNAs and their possible implications in CMS induction provides clues for the intricacies and complexity of this phenomenon. Previously, we characterized an Ogura-CMS line of turnip (Brassica rapa ssp. rapifera) that displays distinct impaired anther development with defective microspore production and premature tapetum degeneration. In the present study, high-throughput sequencing was employed for a genome-wide investigation of miRNAs. Six small RNA libraries of inflorescences collected from the Ogura-CMS line and its maintainer fertile (MF) line of turnip were constructed. A total of 120 pre-miRNAs corresponding to 89 mature miRNAs were identified, including 87 conversed miRNAs and 33 novel miRNAs. Among these miRNAs, the expression of 10 differentially expressed mature miRNAs originating from 12 pre-miRNAs was shown to have changed by more than two-fold between inflorescences of the Ogura-CMS line and inflorescences of the MF line, including 8 down- and 2 up-regulated miRNAs. The expression profiles of the differentially expressed miRNAs were confirmed by stem-loop quantitative real-time PCR. In addition, to identify the targets of the identified miRNAs, a degradome analysis was performed. A total of 22 targets of 25 miRNAs and 17 targets of 28 miRNAs were identified as being involved in the reproductive development for Ogura-CMS and MF lines of turnip, respectively. Negative correlations of expression patterns between partial miRNAs and their targets were detected. Some of these identified targets, such as squamosa promoter-binding-like transcription factor family proteins, auxin response factors and pentatricopeptide repeat-containing proteins, were previously reported to be involved in reproductive development in plants. Taken together, our results can help improve the understanding of miRNA-mediated regulatory pathways that might be involved in CMS occurrence in turnip.

Lipid Metabolism: Critical Roles in Male Fertility and Other Aspects of Reproductive Development in Plants

Fatty acids and their derivatives are essential building blocks for anther cuticle and pollen wall formation. Disruption of lipid metabolism during anther and pollen development often leads to genic male sterility (GMS). To date, many lipid metabolism-related GMS genes that are involved in the formation of anther cuticle, pollen wall, and subcellular organelle membranes in anther wall layers have been identified and characterized. In this review, we summarize recent progress on characterizing lipid metabolism-related genes and their roles in male fertility and other aspects of reproductive development in plants. On the basis of cloned GMS genes controlling biosynthesis and transport of anther cutin, wax, sporopollenin, and tryphine in Arabidopsis, rice, and maize as well as other plant species, updated lipid metabolic networks underlying anther cuticle development and pollen wall formation were proposed. Through bioinformatics analysis of anther RNA-sequencing datasets from three maize inbred lines (Oh43, W23, and B73), a total of 125 novel lipid metabolism-related genes putatively involved in male fertility in maize were deduced. More, we discuss the pathways regulating lipid metabolism-related GMS genes at the transcriptional and post-transcriptional levels. Finally, we highlight recent findings on lipid metabolism-related genes and their roles in other aspects of plant reproductive development. A comprehensive understanding of lipid metabolism, genes involved, and their roles in plant reproductive development will facilitate the application of lipid metabolism-related genes in gene editing, haploid and callus induction, molecular breeding and hybrid seed production in crops.

MYB Superfamily in Brassica napus: Evidence for Hormone-Mediated Expression Profiles, Large Expansion, and Functions in Root Hair Development

MYB proteins are involved in diverse important biological processes in plants. Herein, we obtained the MYB superfamily from the allotetraploid Brassica napus, which contains 227 MYB-related (BnMYBR/Bn1R-MYB), 429 R2R3-MYB (Bn2R-MYB), 22 R1R2R3-MYB (Bn3R-MYB), and two R1R2R2R1/2-MYB (Bn4R-MYB) genes. Phylogenetic analysis classified the Bn2R-MYBs into 43 subfamilies, and the BnMYBRs into five subfamilies. Sequence characteristics and exon/intron structures within each subfamily of the Bn2R-MYBs and BnMYBRs were highly conserved. The whole superfamily was unevenly distributed on 19 chromosomes and underwent unbalanced expansion in B. napus. Allopolyploidy between B. oleracea and B. rapa mainly contributed to the expansion in their descendent B. napus, in which B. rapa-derived genes were more retained. Comparative phylogenetic analysis of 2R-MYB proteins from nine Brassicaceae and seven non-Brassicaceae species identified five Brassicaceae-specific subfamilies and five subfamilies that are lacking from the examined Brassicaceae species, which provided an example for the adaptive evolution of the 2R-MYB gene family alongside angiosperm diversification. Ectopic expression of four Bn2R-MYBs under the control of the viral CaMV35S and/or native promoters could rescue the lesser root hair phenotype of the Arabidopsis thaliana wer mutant plants, proving the conserved negative roles of the 2R-MYBs of the S15 subfamily in root hair development. RNA-sequencing data revealed that the Bn2R-MYBs and BnMYBRs had diverse transcript profiles in roots in response to the treatments with various hormones. Our findings provide valuable information for further functional characterizations of B. napusMYB genes.

The lre-miR159a-LrGAMYB pathway mediates resistance to grey mould infection in Lilium regale

Grey mould is one of the most determinative factors of lily growth and plays a major role in limiting lily productivity. MicroRNA159 (miR159) is a highly conserved microRNA in plants, and participates in the regulation of plant development and stress responses. Our previous studies revealed that lre-miR159a participates in the response of Lilium regale to Botrytis elliptica according to deep sequencing analyses; however, the response mechanism remains unknown. Here, lre-miR159a and its target LrGAMYB gene were isolated from L. regale. Transgenic Arabidopsis overexpressing lre-MIR159a exhibited larger leaves and smaller necrotic spots on inoculation with Botrytis than those of wild-type and overexpressing LrGAMYB plants. The lre-MIR159a overexpression also led to repressed expression of two targets of miR159, AtMYB33 and AtMYB65, and enhanced accumulation of hormone-related genes, including AtPR1, AtPR2, AtNPR1, AtPDF1.2, and AtLOX for both the jasmonic acid and salicylic acid pathways. Moreover, lower levels of H₂ O₂ and O2- were observed in lre-MIR159a transgenic Arabidopsis, which reduced the damage from reactive oxygen species accumulation. Taken together, these results indicate that lre-miR159a positively regulates resistance to grey mould by repressing the expression of its target LrGAMYB gene and activating a defence response.

Insights into the Diversification and Evolution of R2R3-MYB Transcription Factors in Plants

As one of the largest families of transcription factors (TFs) in plants, R2R3-MYB proteins play crucial roles in regulating a series of plant-specific biological processes. Although the diversity of plant R2R3-MYB TFs has been studied previously, the processes and mechanisms underlying the expansion of these proteins remain unclear. Here, we performed evolutionary analyses of plant R2R3-MYB TFs with dense coverage of streptophyte algae and embryophytes. Our analyses revealed that ancestral land plants exhibited 10 subfamilies of R2R3-MYB proteins, among which orthologs of seven subfamilies were present in chlorophytes and charophycean algae. We found that asymmetric gene duplication events in different subfamilies account for the expansion of R2R3-MYB proteins in embryophytes. We further discovered that the largest subfamily of R2R3-MYBs in land plants, subfamily VIII, emerged in the common ancestor of Zygnematophyceae and embryophytes. During plant terrestrialization, six duplication events gave rise to seven clades of subfamily VIII. Subsequently, this TF subfamily showed a tendency for expansion in bryophytes, lycophytes, and ferns and extensively diversified in ancestral gymnosperms and angiosperms in clades VIII-A-1, VIII-D, and VIII-E. In contrast to subfamily VIII, other subfamilies of R2R3-MYB TFs have remained less expanded across embryophytes. The findings regarding phylogenetic analyses, auxiliary motifs, and DNA-binding specificities provide insight into the evolutionary history of plant R2R3-MYB TFs and shed light on the mechanisms underlying the extensive expansion and subsequent sub- and neofunctionalization of these proteins.

ALTERED MERISTEM PROGRAM1 regulates leaf identity independently of miR156-mediated translational repression

In Arabidopsis, loss of the carboxypeptidase ALTERED MERISTEM PROGRAM1 (AMP1) produces an increase in the rate of leaf initiation, an enlarged shoot apical meristem and an increase in the number of juvenile leaves. This phenotype is also observed in plants with reduced levels of miR156-targeted SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors, suggesting that AMP1 might promote SPL activity. However, we found that the amp1 mutant phenotype is only partially corrected by elevated SPL gene expression, and that amp1 has no significant effect on SPL transcript levels, or on the level or the activity of miR156. Although AMP1 has been reported to promote miRNA-mediated translational repression, amp1 did not prevent the translational repression of the miR156 target SPL9 or the miR159 target MYB33. These results suggest that AMP1 regulates vegetative phase change downstream of, or in parallel to, the miR156/SPL pathway, and that it is not universally required for miRNA-mediated translational repression.

Cloning and Functional Analysis of BcMYB101 Gene Involved in Leaf Development in Pak Choi (Brassica rapa ssp. Chinensis)

As one of the largest transcription factor families, MYB transcription factors are widely present, and they are involved in a diverse range of physiological activities in plants, such as leaf development. GAMYB genes belong to the R2R3-MYB subfamily, which includes the MYB33/65/101 gene, and these genes are studied well in seed germination and flowering, but their roles in leaf development are poorly understood. In the current study, we isolated a GAMYB transcription factor from pak choi, BcMYB101, and analyzed its characteristics and function. The sequence structure analysis indicated that BcMYB101 has a highly conserved R2R3 DNA-binding domain in the N-terminal region and three GAMYB-specific motifs (Box1, Box2, and Box3). The expression pattern of diverse tissues revealed that BcMYB101 has a higher transcript level in the petiole, leaf, root, and floral organs. Furthermore, the expression level was significantly elevated after GA (gibberellin) treatment, suggesting that the BcMYB101 response was positively regulated by GA. Subcellular localization exhibited that BcMYB101 was only present in the nuclear region, consistent with the characterization of the transcription factor. The overexpression of BcMYB101 elucidated that BcMYB101 increased leaf number and resulted in downward-curling cauline leaves. Moreover, the virus-induced BcMYB101 silencing displayed that BcMYB101 is involved in the regulation of curly leaves. Furthermore, we discovered that BcMYB101 has two trans-activation activities and one interaction protein, BcTCH4, using a trans-activation activity assay and a yeast two-hybrid assay, respectively. In this study, we firstly isolated the BcMYB101 gene and explored its function in leaf development, thereby providing a solid foundation for further research on the regulatory mechanism of leaf shape in Brassica or other species.

MYB94 and MYB96 additively inhibit callus formation via directly repressing LBD29 expression in Arabidopsis thaliana

Plant somatic cells can be reprogrammed during in vitro culture. Callus induction is the initial step of a typical plant regeneration system. Recent studies showed that auxin-induced callus formation in multiple organs occurs from the pericycle or pericycle-like cells via a root developmental pathway. However, the molecular control of callus formation is largely unknown. Here, two MYB transcription factors, MYB94 and MYB96, were shown to play negative roles in auxin-induced callus formation in Arabidopsis. MYB94 and MYB96 were expressed in the newly formed callus. myb96, myb94, and myb94 myb96 generated more calli than the WT, with myb94 myb96 producing the most. MYB94 and MYB96 repressed expression of LATERAL ORGAN BOUNDARIES-DOMAIN 29 (LBD29) via directly binding to the gene's promoter. The loss of function of LBD29 partly rescued the callus formation defect of myb94 myb96. Our findings found MYB94 and MYB96 to be important repressors of callus formation and MYB94/96-LBD29 as a new regulatory pathway acting in parallel with ARF7/19-LBDs' pathway to modulate in vitro callus formation.

Slowing development restores the fertility of thermo-sensitive male-sterile plant lines

Temperature-sensitive genic male sterility (TGMS) lines are widely used in the breeding of hybrid crops^1,2, but by what means temperature as a general environmental factor reverses the fertility of different TGMS lines remains unknown. Here, we identified an Arabidopsis TGMS line named reversible male sterile (rvms) that is fertile at low temperature (17 °C) and encodes a GDSL lipase. Cytological observations and statistical analysis showed that low temperature slows pollen development. Further screening of restorers of rvms, as well as crossing with a slow-growth line at normal temperature (24 °C), demonstrate that slowing of development overcomes the defects of rvms microspores and allows them to develop into functional pollen. Several other Arabidopsis TGMS lines were identified, and their fertility was also restored by slowing of development. Given that male reproductive development is conserved³, we propose that slowing of development is a general mechanism applicable to the sterility-fertility conversion of TGMS lines from different plant species.

miR159 Represses a Constitutive Pathogen Defense Response in Tobacco

MicroR159 (miR159) regulation of GAMYB expression is highly conserved in terrestrial plants; however, its functional role remains poorly understood. In Arabidopsis (Arabidopsis thaliana), although GAMYB-like genes are constitutively transcribed during vegetative growth, their effects are suppressed by strong and constitutive silencing by miR159. GAMYB expression occurs only if miR159 function is inhibited, which results in detrimental pleiotropic defects, questioning the purpose of the miR159-GAMYB pathway. Here, miR159 function was inhibited in tobacco (Nicotiana tabacum) and rice (Oryza sativa) using miRNA MIM159 technology. Similar to observations in Arabidopsis, inhibition of miR159 in tobacco and rice resulted in pleiotropic defects including stunted growth, implying functional conservation of the miR159-GAMYB pathway among angiosperms. In MIM159 tobacco, transcriptome profiling revealed that genes associated with defense and programmed cell death were strongly activated, including a suite of 22 PATHOGENESIS-RELATED PROTEIN (PR) genes that were 100- to 1,000-fold upregulated. Constitutive expression of a miR159-resistant GAMYB transgene in tobacco resulted in phenotypes similar to that of MIM159 tobacco and activated PR gene expression, verifying the dependence of the above-mentioned changes on GAMYB expression. Consistent with the broad defense response, MIM159 tobacco appeared immune to Phytophthora infection. These findings suggest that the tobacco miR159-GAMYB pathway functions in the biotic defense response, which becomes activated upon miR159 inhibition. However, PR gene expression was not upregulated in Arabidopsis or rice when miR159 was inhibited, suggesting that miR159-GAMYB pathway functional differences exist between species, or factors in addition to miR159 inhibition are required in Arabidopsis and rice to activate this broad defense response.

Research advance in regulation of fruit quality characteristics by microRNAs

MicroRNAs (miRNAs) are short (19–24 nucleotides in length) noncoding RNAs that have a profound effect on gene expression. By completely or almost perfectly base-pairing with their individual target mRNAs they cause mRNA cleavage or repression of translation. As important regulators, miRNAs plays an important role in the regulation of fruit quality. Extensive studies have been reported in fruits, however current studies are mostly focused on the identification of miRNAs and the prediction and validation of target genes. This review summarizes research progress on the role of miRNAs in regulating fruit ripening and senescence and quality characteristics, such as coloration, flavor metabolism, and texture for providing information for future research.

Identification and profiling of microRNAs and differentially expressed genes during anther development between a genetic male-sterile mutant and its wildtype cotton via high-throughput RNA sequencing

Genetic male sterility (GMS) facilitates hybrid seed production in crops including cotton (Gossypium hirsutum). However, the genetic and molecular mechanisms specifically involved in this developmental process are poorly understood. In this study, small RNA sequencing, degradome sequencing, and transcriptome sequencing were performed to analyze miRNAs and their target genes during anther development in a GMS mutant (‘Dong A’) and its fertile wildtype (WT). A total of 80 known and 220 novel miRNAs were identified, 71 of which showed differential expressions during anther development. A further degradome sequencing revealed a total of 117 candidate target genes cleaved by 16 known and 36 novel miRNAs. Based on RNA-seq, 24, 11, and 21 predicted target genes showed expression correlations with the corresponding miRNAs at the meiosis, tetrad and uninucleate stages, respectively. In addition, a large number of differentially expressed genes were identified, most of which were involved in sucrose and starch metabolism, carbohydrate metabolism, and plant hormone signal transduction based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The results of our study provide valuable information for further functional investigations of the important miRNAs and target genes involved in genetic male sterility and advance our understanding of miRNA regulatory functions during cotton anther development.

miRNA-mRNA Integrated Analysis Reveals Roles for miRNAs in a Typical Halophyte, Reaumuria soongorica, during Seed Germination under Salt Stress

MicroRNAs (miRNAs) are endogenous small RNAs that play a crucial role in plant growth, development, and environmental stress responses. Reaumuria soongorica is a typical halophyte that is widely distributed in saline–alkali desert regions. Under salt stress, R. soongorica can complete germination, a critical biological process in the life cycle of seed plants. To identify miRNAs and predict target mRNAs involved in seed germination during salt stress, nine small-RNA libraries were constructed and analyzed from R. soongorica seeds treated with various concentrations of NaCl. We also obtained transcriptome data under the same treatment conditions. Further analysis identified 88 conserved miRNAs representing 25 defined families and discovered 13 novel miRNAs from nine libraries. A co-expression analysis was performed on the same samples to identify putative miRNA–mRNA interactions that were responsive to salt stress. A comparative analysis of expression during germination under 273 (threshold) and 43 mM (optimal) NaCl treatments identified 13 differentially expressed miRNAs and 23 corresponding target mRNAs, while a comparison between 43 mM NaCl and non-salt-stress conditions uncovered one differentially expressed miRNA and one corresponding target mRNA. These results provide basic data for further study of molecular mechanisms involved in the germination of salt-stressed R. soongorica seeds, and also provide a reference for the improvement of salt tolerance during plant germination.

Deep-reaching acceleration of global mean ocean circulation over the past two decades

Ocean circulation redistributes Earth's energy and water masses and influences global climate. Under historical greenhouse warming, regional ocean currents show diverse tendencies, but whether there is an emerging trend of the global mean ocean circulation system is not yet clear. Here, we show a statistically significant increasing trend in the globally integrated oceanic kinetic energy since the early 1990s, indicating a substantial acceleration of global mean ocean circulation. The increasing trend in kinetic energy is particularly prominent in the global tropical oceans, reaching depths of thousands of meters. The deep-reaching acceleration of the ocean circulation is mainly induced by a planetary intensification of surface winds since the early 1990s. Although possibly influenced by wind changes associated with the onset of a negative Pacific decadal oscillation since the late 1990s, the recent acceleration is far larger than that associated with natural variability, suggesting that it is principally part of a long-term trend.

MYB81, a microspore-specific GAMYB transcription factor, promotes pollen mitosis I and cell lineage formation in Arabidopsis

Sexual reproduction in flowering plants relies on the production of haploid gametophytes that consist of germline and supporting cells. During male gametophyte development, the asymmetric mitotic division of an undetermined unicellular microspore segregates these two cell lineages. To explore genetic regulation underlying this process, we screened for pollen cell patterning mutants and isolated the heterozygous myb81-1 mutant that sheds ~50% abnormal pollen. Typically, myb81-1 microspores fail to undergo pollen mitosis I (PMI) and arrest at polarized stage with a single central vacuole. Although most myb81-1 microspores degenerate without division, a small fraction divides at later stages and fails to acquire correct cell fates. The myb81-1 allele is transmitted normally through the female, but rarely through pollen. We show that myb81-1 phenotypes result from impaired function of the GAMYB transcription factor MYB81. The MYB81 promoter shows microspore-specific activity and a MYB81-RFP fusion protein is only expressed in a narrow window prior to PMI. Ectopic expression of MYB81 driven by various promoters can severely impair vegetative or reproductive development, reflecting the strict microspore-specific control of MYB81. Our data demonstrate that MYB81 has a key role in the developmental progression of microspores, enabling formation of the two male cell lineages that are essential for sexual reproduction in Arabidopsis.

Tapetum-Dependent Male Meiosis Progression in Plants: Increasing Evidence Emerges

In higher plants, male meiosis is a key process during microsporogenesis and is crucial for male fertility and seed set. Meiosis involves a highly dynamic organization of chromosomes and cytoskeleton and specifically takes place within sexual cells. However, studies in multiple plant species have suggested that the normal development of tapetum, the somatic cell layer surrounding the developing male meiocytes, is indispensable for the completion of the male meiotic cell cycle. Disrupted tapetum development causes alterations in the expression of a large range of genes involved in male reproduction. Moreover, recent experiments suggest that small RNAs (sRNAs) present in the anthers, including microRNAs (miRNAs) and phased, secondary, small interfering RNAs (phasiRNAs), play a potential but important role in controlling male meiosis, either by influencing the expression of meiotic genes in the meiocytes or through other unclear mechanisms, supporting the hypothesis that male meiosis is non-cell autonomously regulated. In this mini review, we summarize the recorded meiotic defects that occur in plants with defective tapetum development in both Arabidopsis and crops. Thereafter, we outline the latest understanding on the molecular mechanisms that potentially underpin the tapetum-dependent regulation of male meiosis, and we especially discuss the regulatory role of sRNAs. At the end, we propose several outstanding questions that should be addressed in future studies.

Rice Female Meiosis: Genome-Wide mRNA, Small RNA, and DNA Methylation Analysis During Ovule Development

Meiosis is an essential process in sexual life cycle, not only for the genomic stability maintenance but also for the genetic diversity creation through recombination. In rice ovule, megaspore mother cells undergo meiosis to form megaspores; then the functional megaspore performs three rounds of mitoses to form female gametophyte. However, the mechanism of gene expression and regulation in female meiosis process is still poorly understood. As important gene regulatory factors, miRNAs and DNA methylation are widely involved in plant meiosis and ovule development. In order to systematically study the potential mechanism of gene expression and regulation in female meiosis, ovules at megaspore mother cell meiosis stage, functional megaspore mitosis stage, and mature female gametophytes are collected to perform genome-wide RNA sequencing, small RNA sequencing, and bisulfite sequencing. Through bioinformatics analysis, we obtained many differentially expressed genes, miRNAs, and differentially methylated genes related to female meiosis. These data may provide important clues for further revealing the mechanism of female meiosis in rice.

Abnormal tapetum development in hermaphrodites of an androdioecious tree, Tapiscia sinensis

Tapiscia sinensis Oliv. (Tapisciaceae) has been proven to be a functional androdioecious species with both male and hermaphroditic individuals, and the pollen viability of males is far higher than that of hermaphrodites. To better understand the causes of the low pollen viability in hermaphroditic flowers, different stages of anther development were observed. We found that hermaphroditic flowers exhibit abnormal tapetum development, resulting in low pollen viability. To clarify the underlying molecular mechanism of abnormal tapetum development in hermaphrodites, quantitative real-time PCR analyses were performed. The results revealed that the expression levels of an important transcription factor for tapetum development and function, T. sinensis DYSFUNCTIONAL TAPETUM1 (TsDYT1), and its potential downstream regulatory genes T. sinensis DEFECTIVE in TAPETAL DEVELOPMENT and FUNCTION1 (TsTDF1), T. sinensis ABORTED MICROSPORE (TsAMS) and T. sinensis MALE STERILITY 1 (TsMS1) were all significantly downregulated in hermaphrodites compared with males at some key stages of anther development. The amino acid sequence similarity, expression pattern, gene structure and subcellular localization of these genes were analyzed, and the results indicated functional conservation between T. sinensis and homologues in Arabidopsis thaliana. Next, rapid amplification of cDNA end and thermal asymmetric interlaced PCR were employed to clone the full-length cDNA and promoter sequences of these genes, respectively. In addition, results of yeast two-hybrid analysis showed that TsDYT1 can form heterodimers with TsAMS, and yeast one-hybrid analysis demonstrated that TsDYT1 directly binds to the promoter regions of TsTDF1 and TsMS1. TsTDF1 can directly regulate expression of TsAMS, suggesting that a functionally conserved pathway exists between A. thaliana and T. sinensis to regulate tapetum development. In conclusion, the results suggest that abnormal expression of core transcription factors for tapetum development, including TsDYT1, TsTDF1, TsAMS and TsMS1, plays an important role in the abnormal development of the tapetum in T. sinensis hermaphrodites. Furthermore, a hermaphroditic tapetum with abnormal function causes the low pollen viability of hermaphroditic trees. Our results provide new insight into our understanding of the underlying mechanism of why pollen viability is much higher in males than hermaphrodites of the androdioecious tree T. sinensis.

A Natural Variant of miR397 Mediates a Feedback Loop in Circadian Rhythm

MicroRNAs (miRNAs) are small noncoding RNAs of ∼21 nt in length, which have regulatory roles in many biological processes. In animals, proper functioning of the circadian clock, which is closely linked to the fitness of almost all living organisms, is regulated by miRNAs. However, to date, there have been no reports of the roles of miRNA in regulation of the plant circadian rhythm. Here, we report a natural variant of miR397 that lengthens the circadian period and controls flowering time in Arabidopsis (Arabidopsis thaliana). Highly conserved among angiosperms, the miRNA miR397 has two members in Arabidopsis: miR397a and miR397b. However, only miR397b significantly delayed flowering. Our results suggest that miR397b controls flowering by targeting CASEIN KINASE II SUBUNIT BETA3 (CKB3), in turn modulating the circadian period of CIRCADIAN CLOCK ASSOCIATED1 (CCA1). We further demonstrated that CCA1 directly bound to the promoter of MIR397B and suppressed its expression, forming a miR397b-CKB3-CCA1 circadian regulation feedback circuit. Evolutionary analysis revealed that miR397b is a newly evolved genetic variant in Arabidopsis, and the miR397b targeting mode may have a role in enhancing plant fitness. Our results provide evidence for miRNA-mediated circadian regulation in plants and suggest the existence of a feedback loop to manipulate plant flowering through the regulation of circadian rhythm.

Integrated mRNA, MicroRNA Transcriptome and Degradome Analyses Provide Insights into Stamen Development in Moso Bamboo

A flower is an essential organ for sexual reproduction in flowering plants, which has been extensively studied in model plants. In this study, we used transcriptomic, small RNA and degradome analyses to characterize key microRNAs (miRNAs) and their targets in floral organs of moso bamboo. In total, we identified 13,051 differentially expressed genes and 109 known miRNAs from 26 miRNA families. We aligned the miRNAs to known miRNA databases and revealed some conserved as well as novel miRNAs. Sixteen conserved miRNAs were specifically and highly expressed in stamens, including miRNA159 and miRNA166. In situ hybridization shows that miRNA159 plays a key role in the regulation of stamen development, and the expression levels of its targets PheMYB98 and PheMYB42 were low. Furthermore, Phe-MIRNA159 partially recovers phenotypes of mir159ab double mutant. Overexpression of Phe-MIR159 could cause failure in anther dehisce, and the mature pollens could not be dispersed and further reduce fertility in Arabidopsis. Semi-thin section result shows that anther endothelial layer of Phe-MIRNA159 overexpressing lines is lack of secondary thickening, resulting in limited force for anther opening. Phe-miR159 may regulate the expression of genes related to secondary thickening through negative regulation of AtMYB33, affecting the anther dehiscence. Taken together, this study provides insights regarding molecular networks underlying floral organs development of moso bamboo.

Discovery and Profiling of microRNAs at the Critical Period of Sex Differentiation in Xanthoceras sorbifolium Bunge

Research Highlights: The critical period of sex differentiation in Xanthoceras sorbifolium was investigated. Multiple microRNAs (miRNAs) were identified to influence female and male flower development, with some complementary functions. Background and Objectives: Xanthoceras sorbifolium Bunge is widely cultivated owing to its multipurpose usefulness. However, as a monoecious plant, the low female–male flowers ratio and consequent low seed yield are the main bottlenecks for industrial-scale development of seed utilization. MiRNAs play crucial regulatory roles in flower development and sex differentiation; therefore, we evaluated the roles of miRNAs in the critical period of sex differentiation in X. sorbifolium. Materials and Methods: Four small RNA libraries for female and male flower buds of the critical period of sex differentiation were constructed from paraffin-embedded sections. The miRNAs were characterized by high-throughput sequencing, and differentially expressed miRNAs were validated by reverse transcription-quantitative polymerase chain reaction. Results: There were obvious differences in male and female pistil and stamen flower buds, with elongated inflorescence and clear separation of flower buds marking the critical period of sex differentiation. A total of 1619 conserved miRNAs (belonging to 34 families) and 219 novel miRNAs were identified. Among these, 162 conserved and 14 novel miRNAs exhibited significant differential expression in the four libraries, and 1677 putative target genes of 112 differentially expressed miRNAs were predicted. These target genes were involved in diverse developmental and metabolic processes, including 17 miRNAs directly associated with flower and gametophyte development, mainly associated with carbohydrate metabolism and glycan biosynthesis and metabolism pathways. Some miRNA functions were confirmed, and others were found to be complemented. Conclusions: Multiple miRNAs closely related to sex differentiation in X. sorbifolium were identified. The theoretical framework presented herein might guide sex ratio regulation to enhance seed yield.

Identification and characterization of microRNAs in phloem and xylem from ramie (Boehmeria nivea)

Ramie (Boehmeria nivea) is a widely cropped species in southern China due to its high economic value of natural fiber for industry. Development of phloem and xylem is key evidence for generating fiber. However, the MicroRNA (miRNA) profiles of phloem and xylem in ramie have not been reported yet. miRNA belong to a small RNA family which has been recognized as an important regulator for various biological processes. In the present study, we aimed to identify differently expressed miRNAs between phloem and xylem in adult ramie. The results showed that 137 and 122 unique conserved miRNAs were identified from phloem and xylem libraries, respectively. Meanwhile, 4 novel miRNAs were identified from ramie by miRDeep2. Of these miRNAs, 77 conserved miRNAs in ramie were differentially expressed. Among the differentially expressed miRNAs, 44 miRNAs and 33 miRNAs were up-regulated and down-regulated in phloem compared to that in xylem, respectively. The functions of differentially expressed miRNAs were associated with regulating the development and differentiation of phloem and xylem. The present study provides a glance of miRNA profiles for further understanding of miRNA role in ramie development.

Characterization of watermelon anther and its programmed cell death-associated events during dehiscence under cold stress

The 'neglected' thermophile fruit crop of watermelon was first used as a model crop to study the PCD associated with anther dehiscence in cold-exposed condition during anther development. Anther dehiscence ensures normal pollen release and successful fertilization at fruit-setting stages in flowering plants. However, most researches pertinent to anther dehiscence are centered on model plant and/or major field crops under optimal growth condition. Due to anther indehiscence in cold condition, crop plants of thermophile tropical or subtropical fruit crops fail to accomplish timely pollination and fertilization, resulting in a great yield loss annually. Herein, we developed an ideal model crop for studying the programmed cell death (PCD) associated with anther dehiscence under low-temperature stress using the S-shaped spiral anther in watermelon as instead. Our results revealed that, including the tapetal cell layers, both cells of the interlocular septum and the stomium were blocked in PCD associated with anther dehiscence at 15 °C. Likewise, TUNEL assays visualized the evidence that low temperature at 15 °C interferes with not only the PCD of tapetal cells, but also the PCD of interlocular septum and stomium. Furthermore, the expressions of genes correlated with PCD of tapetum and stomium were significantly inhibited at 15 °C, suggesting that low temperature affects anther dehiscence by inhibiting PCD of sporophytic tissue-related gene expressions. The findings of the current research provide mechanistic insights into anther indehiscence leading to poor fruit-setting for thermophile fruit crop such as watermelon under adverse cold condition in flowering.

Evolution of MIR159/319 genes in Brassica campestris and their function in pollen development

MIR159/319 have conserved evolution and diversified function after WGT in Brassica campestris, both of them can lead pollen vitality and germination abnormality, Bra-MIR319c also can function in flower development. MiR159 and miR319 are extensively studied highly conserved microRNAs which play roles in vegetative development, reproduction, and hormone regulation. In this study, the effects of whole-genome triplication (WGT) on the evolution of the MIR159/319 family and the functional diversification of the genes were comprehensively investigated in Brassica campestris. We identified 11 MIR159/319 genes in B. campestris, which produced five mature sequences. After analyzing the precursor sequences and phylogenetic tree, we found that Bra-MIR159/319 have evolutionary conservatism. Furthermore, Bra-MIR159/319 show functional diversification after WGT, as indicated by their expression patterns and the cis-element in their promoter. GUS signal showed that Bra-MIR159a and Bra-MIR319c can be expressed in anther but in different development stages. In B. campestris, overexpressed MIR159a and MIR319c contribute to late anther development and promote pollen abortion. Moreover, Bra-MIR319c can partially assume the function of MIR319a in flower development.

Identification of microRNAs and relative target genes in Moringa oleifera leaf and callus

MicroRNAs, a class of small, non-coding RNAs, play important roles in plant growth, development and stress response by negatively regulating gene expression. Moringa oleifera Lam. plant has many medical and nutritional uses; however, little attention has been dedicated to its potential for the bio production of active compounds. In this study, 431 conserved and 392 novel microRNA families were identified and 9 novel small RNA libraries constructed from leaf, and cold stress treated callus, using high-throughput sequencing technology. Based on the M. oleifera genome, the microRNA repertoire of the seed was re-evaluated. qRT-PCR analysis confirmed the expression pattern of 11 conserved microRNAs in all groups. MicroRNA159 was found to be the most abundant conserved microRNA in leaf and callus, while microRNA393 was most abundantly expressed in the seed. The majority of predicted microRNA target genes were transcriptional factors involved in plant reproduction, growth/development and abiotic/biotic stress response. In conclusion, this is the first comprehensive analysis of microRNAs in M. oleifera leaf and callus which represents an important addition to the existing M. oleifera seed microRNA database and allows for possible exploitation of plant microRNAs induced with abiotic stress, as a tool for bio-enrichment with pharmacologically important phytochemicals.

Using Transcriptome Analysis to Screen for Key Genes and Pathways Related to Cytoplasmic Male Sterility in Cotton (Gossypium hirsutum L.)

Cotton (Gossypium hirsutum L.) is one of the most important cash crops worldwide. Cytoplasmic male sterility (CMS) is an excellent breeding system for exploitation of heterosis, which has great potential to increase crop yields. To understand the molecular mechanism of CMS in cotton, we compared transcriptome, cytomorphological, physiological and bioinformatics data between the CMS line C2P5A and its maintainer line C2P5B. By using high-throughput sequencing technology, 178,166 transcripts were assembled and 2013 differentially expression genes (DEGs) were identified at three different stages of C2P5A anther development. In this study, we identified DEGs associated with reactive oxygen species (ROS), peroxisomes, aldehyde dehydrogenases (ALDH), cytochrome oxidase subunit VI, and cytochrome P450, and DEGs associated with tapetum development, Jojoba acyl-CoA reductase-related male sterility protein, basic helix-loop-helix (bHLH) and MYB transcription factors. The abnormal expression of one of these genes may be responsible for the CMS C2P5A line. In gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, DEGs were mainly related to carbohydrate metabolism, amino acid metabolism, transport and catabolism, and signal transduction. Carbohydrate metabolism provides energy for anther development, starch and sucrose metabolism, fatty acid biosynthesis and metabolism and ascorbate and aldarate metabolism. These results showed that numerous genes and multiple complex metabolic pathways regulate cotton anther development. Weighted correlation network analysis (WGCNA) indicated that three modules, ‘turquoise,’ ‘blue,’ and ‘green,’ were specific for the CMS C2P5A line. The ‘turquoise’ and ‘blue’ modules were mainly related to carbohydrate metabolism, amino acid metabolism, energy metabolism, peroxisomes, pyruvate metabolism as well as fatty acid degradation. The ‘green’ module was mainly related to energy metabolism, carbon metabolism, translation, and lipid metabolism. RNA-sequencing and WGCNA polymerization modules were screened for key genes and pathways related to CMS in cotton. This study presents a new perspective for further research into the metabolic pathways of pollen abortion in the CMS C2P5A line and also provides a theoretical basis for its breeding and production.

Comprehensive Insight into Gibberellin- and Jasmonate-Mediated Stamen Development

In flowering plants, proper development of male generative organs is required for successful sexual reproduction. Stamen primordia arise in the third whorl of floral organs and subsequently differentiate into filaments and anthers. The early phase of stamen development, in which meiosis occurs, is followed by a late developmental phase, which consists of filament elongation coordinated with pollen maturation, anther dehiscence and finally viable pollen grain release. Stamen development and function are modulated by phytohormones, with a key role of gibberellins (GAs) and jasmonates (JAs). Long-term, extensive investigations, mainly involving GA/JA-deficient and GA/JA-response mutants, have led to a better understanding of the hormone-dependent molecular mechanisms of stamen development. In several species, the principal functions of GAs are to stimulate filament elongation through increased cell elongation and to promote anther locule opening. In the GA-dependent regulation of early stamen development, both the tapetum and developing pollen were identified as major targets. JAs mainly control the late stages of stamen development, such as filament elongation, viable pollen formation and anther dehiscence. A hierarchical relationship between GAs and JAs was recognized mainly in the control of late stamen development. By repressing DELLA proteins, GAs modulate the transcriptional activity of JA biosynthesis genes to promote JA production. A high level of JAs induces a complex of transcription factors crucial for normal stamen development.

Different expression pattern of flowering pathway genes contribute to male or female organ development during floral transition in the monoecious weed Ambrosia artemisiifolia L. (Asteraceae)

The highly allergenic and invasive weed Ambrosia artemisiifolia L. is a monoecius plant with separated male and female flowers. The genetic regulation of floral morphogenesis is a less understood field in the reproduction biology of this species. Therefore the objective of this work was to investigate the genetic control of sex determination during floral organogenesis. To this end, we performed a genome-wide transcriptional profiling of vegetative and generative tissues during the plant development comparing wild-growing and in vitro cultivated plants. RNA-seq on Illumina NextSeq 500 platform with an integrative bioinformatics analysis indicated differences in 80 floral gene expressions depending on photoperiodic and endogenous initial signals. Sex specificity of genes was validated based on RT-qPCR experiments. We found 11 and 16 uniquely expressed genes in female and male transcriptomes that were responsible particularly to maintain fertility and against abiotic stress. High gene expression of homologous such as FD, FT, TFL1 and CAL, SOC1, AP1 were characteristic to male and female floral meristems during organogenesis. Homologues transcripts of LFY and FLC were not found in the investigated generative and vegetative tissues. The repression of AP1 by TFL1 homolog was demonstrated in male flowers resulting exclusive expression of AP2 and PI that controlled stamen and carpel formation in the generative phase. Alterations of male and female floral meristem differentiation were demonstrated under photoperiodic and hormonal condition changes by applying in vitro treatments.

Small RNA profiling from meiotic and post-meiotic anthers reveals prospective miRNA-target modules for engineering male fertility in sorghum

Understanding male gametophyte development is essential to augment hybrid production in sorghum. Although small RNAs are known to critically influence anther/pollen development, their roles in sorghum reproduction have not been deciphered yet. Here, we report small RNA profiling and high-confidence annotation of microRNAs (miRNAs) from meiotic and post-meiotic anthers in sorghum. We identified 262 miRNAs (82 known and 180 novel), out of which 58 (35 known and 23 novel) exhibited differential expression between two stages. Out of 35 differentially expressed known miRNAs, 13 are known to regulate anther/pollen development in other plant species. We also demonstrated conserved spatiotemporal patterns of 21- and 24-nt phasiRNAs and their respective triggers, miR2118 and miR2275, in sorghum anthers as evidenced in other monocots. miRNA target identification yielded 5622 modules, of which 46 modules comprising 16 known and 8 novel miRNA families with 38 target genes are prospective candidates for engineering male fertility in grasses.

Comparative Transcriptome Analysis of Gene Expression Patterns in Tomato Under Dynamic Light Conditions

Plants grown under highly variable natural light regimes differ strongly from plants grown under constant light (CL) regimes. Plant phenotype and adaptation responses are important for plant biomass and fitness. However, the underlying regulatory mechanisms are still poorly understood, particularly from a transcriptional perspective. To investigate the influence of different light regimes on tomato plants, three dynamic light (DL) regimes were designed, using a CL regime as control. Morphological, photosynthetic, and transcriptional differences after five weeks of treatment were compared. Leaf area, plant height, shoot /root weight, total chlorophyll content, photosynthetic rate, and stomatal conductance all significantly decreased in response to DL regimes. The biggest expression difference was found between the treatment with the highest light intensity at the middle of the day with a total of 1080 significantly up-/down-regulated genes. A total of 177 common differentially expressed genes were identified between DL and CL conditions. Finally, significant differences were observed in the levels of gene expression between DL and CL treatments in multiple pathways, predominantly of plant–pathogen interactions, plant hormone signal transductions, metabolites, and photosynthesis. These results expand the understanding of plant development and photosynthetic regulations under DL conditions by multiple pathways.

Overexpression of a stamen-specific R2R3-MYB gene BcMF28 causes aberrant stamen development in transgenic Arabidopsis

In flowering plants, stamen development is a complex multistage process, which is highly regulated by a series of transcription factors. In this study, BcMF28, which encodes a R2R3-MYB transcription factor, was isolated from Brassica campestris. BcMF28 is localized in the nucleus and cytoplasm, and acts as a transcriptional activator. Quantitative real-time PCR and promoter activity analysis revealed that BcMF28 was predominately expressed in inflorescences. The expression of BcMF28 was specifically detected in tapetum, developing microspores, anther endothecium, and filaments during late stamen development. The overexpression of BcMF28 in Arabidopsis resulted in aberrant stamen development, including filament shortening, anther indehiscence, and pollen abortion. Detailed analysis of anther development in transgenic plants revealed that the degeneration of septum and stomium did not occur, and endothecium lignification was affected. Furthermore, the expression levels of genes involved in the phenylpropanoid metabolism pathway were altered in BcMF28-overexpressing transgenic plants. Our results suggest that BcMF28 plays an important regulatory role during late stamen development.

Comparative Analysis of miRNA Abundance Revealed the Function of Vvi-miR828 in Fruit Coloring in Root Restriction Cultivation Grapevine (Vitis vinifera L.)

Root restriction cultivation leads to early maturation and quality improvement, especially in the anthocyanin content in grapevine. However, the molecular mechanisms that underlie these changes have not been thoroughly elucidated. In this study, four small RNA libraries were constructed, which included the green soft stage (GS) and ripe stage (RS) of ‘Muscat’ (Vitis vinifera L.) grape berries that were grown under root restriction (RR) and in traditional cultivation (no root restriction, CK). A total of 162 known miRNAs and 14 putative novel miRNAs were detected from the four small RNA libraries by high-throughput sequencing. An analysis of differentially expressed miRNAs (DEMs) revealed that 13 miRNAs exhibited significant differences in expression between RR and CK at the GS and RS stages, respectively. For different developmental stages of fruit, 23 and 34 miRNAs showed expression differences between the GS and RS stages in RR and CK, respectively. The expression patterns of the eight DEMs and their targets were verified by qRT-PCR, and the expression profiles of target genes were confirmed to be complementary to the corresponding miRNAs in RR and CK. The function of Vvi-miR828, which showed the down regulated expression in the RS stage under root restriction, was identified by gene transformation in Arabidopsis. The anthocyanin content significantly decreased in transgenic lines, which indicates the regulatory capacity of Vvi-miR828 in fruit coloration. The miRNA expression pattern comparison between RR and CK might provide a means of unraveling the miRNA-mediated molecular process regulating grape berry development under root restricted cultivation.

Joint RNA-Seq and miRNA Profiling Analyses to Reveal Molecular Mechanisms in Regulating Thickness of Pod Canopy in Brassica napus

Oilseed rape (Brassica napus) is the second largest oilseed crop worldwide. As an architecture component of B. napus, thickness of pod canopy (TPC) plays an important role in yield formation, especially under high-density cultivation conditions. However, the mechanisms underlying the regulation of TPC remain unclear. RNA and microRNA (miRNA) profiling of two groups of B. napus lines with significantly different TPC at the bolting with a tiny bud stage revealed differential expressions of numerous genes involved in nitrogen-related pathways. Expression of several nitrogen-related response genes, including ASP5, ASP2, ASN3, ATCYSC1, PAL2, APT2, CRTISO, and COX15, was dramatically changed in the thick TPC lines compared to those in the thin TPC lines. Differentially expressed miRNAs also included many involved in nitrogen-related pathways. Expression of most target genes was negatively associated with corresponding miRNAs, such as miR159, miR6029, and miR827. In addition, 12 (including miR319, miR845, and miR158) differentially expressed miRNAs between two plant tissues sampled (stem apex and flower bud) were identified, implying that they might have roles in determining overall plant architecture. These results suggest that nitrogen signaling may play a pivotal role in regulating TPC in B. napus.

Biology and Function of miR159 in Plants

MicroR159 (miR159) is ancient, being present in the majority of land plants where it targets a class of regulatory genes called GAMYB or GAMYB-like via highly conserved miR159-binding sites. These GAMYB genes encode R2R3 MYB domain transcription factors that transduce the gibberellin (GA) signal in the seed aleurone and the anther tapetum. Here, GAMYB plays a conserved role in promoting the programmed cell death of these tissues, where miR159 function appears weak. By contrast, GAMYB is not involved in GA-signaling in vegetative tissues, but rather its expression is deleterious, leading to the inhibition of growth and development. Here, the major function of miR159 is to mediate strong silencing of GAMYB to enable normal growth. Highlighting this requirement of strong silencing are conserved RNA secondary structures associated with the miR159-binding site in GAMYB mRNA that promotes miR159-mediated repression. Although the miR159-GAMYB pathway in vegetative tissues has been implicated in a number of different functions, presently no conserved role for this pathway has emerged. We will review the current knowledge of the different proposed functions of miR159, and how this ancient pathway has been used as a model to help form our understanding of miRNA biology in plants.

Systematic Analysis of MYB Family Genes in Potato and Their Multiple Roles in Development and Stress Responses

The MYB proteins represent a large family of transcription factors and play important roles in development, senescence, and stress responses in plants. In the current study, 233 MYB transcription factor-encoding genes were identified and analyzed in the potato genome, including 119 R1-MYB, 112 R2R3-MYB, and two R1R2R3-MYB members. R2R3-MYB is the most abundant MYB subclass and potato R2R3-MYB members together with their Arabidopsis homologs were divided into 35 well-supported subgroups as the result of phylogenetic analyses. Analyses on gene structure and protein motif revealed that members from the same subgroup shared similar exon/intron and motif organization, further supporting the results of phylogenetic analyses. Evolution of the potato MYB family was studied via syntenic analysis. Forty-one pairs of StMYB genes were predicted to have arisen from tandem or segmental duplication events, which played important roles in the expansion of the StMYB family. Expression profiling revealed that the StMYB genes were expressed in various tissues and several StMYB genes were identified to be induced by different stress conditions. Notably, StMYB030 was found to act as the homolog of AtMYB44 and was significantly up-regulated by salt and drought stress treatments. Furthermore, overexpression of StMYB030 in Arabidopsis enhanced salt stress tolerance of transgenic plants. The results from this study provided information for further functional analysis and for crop improvements through genetic manipulation of these StMYB genes.