The R2R3 MYB transcription factors FOUR LIPS and MYB88 regulate female reproductive development

Gene Regulatory Network Controlling Flower Development in Spinach (Spinacia oleracea L.)

Spinach (Spinacia oleracea L.) is a dioecious, diploid, wind-pollinated crop cultivated worldwide. Sex determination plays an important role in spinach breeding. Hence, this study aimed to understand the differences in sexual differentiation and floral organ development of dioecious flowers, as well as the differences in the regulatory mechanisms of floral organ development of dioecious and monoecious flowers. We compared transcriptional-level differences between different genders and identified differentially expressed genes (DEGs) related to spinach floral development, as well as sex-biased genes to investigate the flower development mechanisms in spinach. In this study, 9189 DEGs were identified among the different genders. DEG analysis showed the participation of four main transcription factor families, MIKC_MADS, MYB, NAC, and bHLH, in spinach flower development. In our key findings, abscisic acid (ABA) and gibberellic acid (GA) signal transduction pathways play major roles in male flower development, while auxin regulates both male and female flower development. By constructing a gene regulatory network (GRN) for floral organ development, core transcription factors (TFs) controlling organ initiation and growth were discovered. This analysis of the development of female, male, and monoecious flowers in spinach provides new insights into the molecular mechanisms of floral organ development and sexual differentiation in dioecious and monoecious plants in spinach.

Transcriptional chronology reveals conserved genes involved in pennate diatom sexual reproduction

Sexual reproduction is a major driver of adaptation and speciation in eukaryotes. In diatoms, siliceous microalgae with a unique cell size reduction-restitution life cycle and among the world's most prolific primary producers, sex also acts as the main mechanism for cell size restoration through the formation of an expanding auxospore. However, the molecular regulators of the different stages of sexual reproduction and size restoration are poorly explored. Here, we combined RNA sequencing with the assembly of a 55 Mbp reference genome for Cylindrotheca closterium to identify patterns of gene expression during different stages of sexual reproduction. These were compared with a corresponding transcriptomic time series of Seminavis robusta to assess the degree of expression conservation. Integrative orthology analysis revealed 138 one-to-one orthologues that are upregulated during sex in both species, among which 56 genes consistently upregulated during cell pairing and gametogenesis, and 11 genes induced when auxospores are present. Several early, sex-specific transcription factors and B-type cyclins were also upregulated during sex in other pennate and centric diatoms, pointing towards a conserved core regulatory machinery for meiosis and gametogenesis across diatoms. Furthermore, we find molecular evidence that the pheromone-induced cell cycle arrest is short-lived in benthic diatoms, which may be linked to their active mode of mate finding through gliding. Finally, we exploit the temporal resolution of our comparative analysis to report the first marker genes for auxospore identity called AAE1-3 ("Auxospore-Associated Expression"). Altogether, we introduce a multi-species model of the transcriptional dynamics during size restoration in diatoms and highlight conserved gene expression dynamics during different stages of sexual reproduction.

An R2R3 MYB transcription factor PsFLP regulates the symmetric division of guard mother cells during stomatal development in Pisum sativum

MYB transcriptional regulators belong to one of the most significant transcription factors families in plants, among which R2R3-MYB transcription factors are involved in plant growth and development, hormone signal transduction, and stress response. Two R2R3-MYB transcription factors, FLP and its paralogous AtMYB88, redundantly regulate the symmetrical division of guard mother cells (GMCs), and abiotic stress response in Arabidopsis thaliana. Only one orthologue gene of FLP was identified in pea (Pisum sativum FLP; PsFLP). In this study, we explored the gene function of PsFLP by virus-induced gene silencing (VIGS) technology. The phenotypic analysis displayed that the silencing of PsFLP expression led to the abnormal development of stomata and the emergence of multiple guard cells tightly united. In addition, the abnormal stomata of flp could be fully rescued by PsFLP driven by the FLP promoter. In conclusion, the results showed that PsFLP plays a conservative negative role in regulating the symmetric division of GMC during stomatal development. Based on real-time quantitative PCR, the relative expressions of AAO3, NCED3, and SnRK2.3 significantly increased in the flp pFLP::PsFLP plants compared to mutant, indicating that PsFLP might be involved in drought stress response. Thus, PsFLP regulates the genes related to cell cycle division during the stomatal development of peas and participates in response to drought stress. The study provides a basis for further research on its function and application in leguminous crop breeding.

FOUR LIPS plays a role in meristemoid-to-GMC fate transition during stomatal development in Arabidopsis

Meristemoids, which are stomatal precursor cells, exhibit self-renewal and differentiation abilities. However, the only known core factor associated with meristemoid division termination and fate transition is the heterodimer formed by the basic helix-loop-helix proteins MUTE and SCREAMs (SCRMs). FOUR LIPS (FLP), a well-known transcription factor that restricts guard mother cell (GMC) division, is a direct target of MUTE. Whether FLP involves in meristemoid differentiation is unknown. Through sensitized genetic screening of flp-1, we identified a mute-like (mutl) mutant with arrested meristemoids. The mutant carried a novel allele of the MUTE locus, i.e., mute-4. Intriguingly, mute-4 is a hypomorphic allele that exhibits wild-type appearance with slightly delayed meristemoid-to-GMC transition, whereas it renders an unexpected mutl epidermis with most meristemoids arrested and very few stomata when combined with flp (flp mute-4), suggesting that FLP is a positive regulator during this transition process. Consistently, the expression of FLP increased during GMC commitment, and the number of cells at this stage was markedly increased in flp. flp scrm double mutants produced arrested meristemoids similar to mute, and FLP was able to interact physically with SCRM. Taken together, our results demonstrate that FLP functions together with MUTE and SCRMs to direct meristemoid-to-GMC fate transition.

Transcriptomic analysis reveals the role of FOUR LIPS in response to salt stress in rice

An R2R3-MYB transcription factor FOUR LIPS associated with B-type Cyclin-Dependent Kinase 1;1 confers salt tolerance in rice. The Arabidopsis FOUR LIPS (AtFLP), an R2R3 MYB transcription factor, acts as an important stomatal development regulator. Only one orthologue protein of AtFLP, Oryza sativa FLP (OsFLP), was identified in rice. However, the function of OsFLP is largely unknown. In this study, we conducted RNA-seq and ChIP-seq to investigate the potential role of OsFLP in rice. Our results reveal that OsFLP is probably a multiple functional regulator involved in many biological processes in growth development and stress responses in rice. However, we mainly focus on the role of OsFLP in salt stress response. Consistently, phenotypic analysis under salt stress conditions showed that osflp exhibited significant sensitivity to salt stress, while OsFLP over-expression lines displayed obvious salt tolerance. Additionally, Yeast one-hybrid assay and electrophoretic mobility shift assay (EMSA) showed that OsFLP directly bound to the promoter region of Oryza sativa B-type Cyclin-Dependent Kinase 1;1 (OsCDKB1;1), and the expression of OsCDKB1;1 was repressed in osflp. Disturbing the expression of OsCDKB1;1 remarkably enhanced the tolerance to salt stress. Taken together, our findings reveal a crucial function of OsFLP regulating OsCDKB1;1 in salt tolerance and largely extend the knowledge about the role of OsFLP in rice.

MdFLP enhances drought tolerance by regulating MdNAC019 in self-rooted apple stocks

Self-rooted apple stocks are widely used for the production of apples worldwide. However, self-rooted apple stocks are weak due to shallow roots and poor grounding, resulting in poor drought resistance. Therefore, it is essential to understand the molecular mechanisms to develop self-rooted apple stock cultivars with drought resistance. We reported that MdFLP, an R2R3-MYB transcription factor, directly binds to the promoter of MdNAC019, activating its transcription and consequently enhancing drought tolerance in self-rooted apple stocks. In addition, MdFLP indirectly activates the transcriptional expression of abiotic stress-related genes, namely, MdERF6 and MdZAT10. The plants overexpressing MdFLP displayed stronger drought tolerance, whereas MdFLP-RNAi plants showed weak drought tolerance compared with non-transgenic "Gala" plants, indicating that MdFLP regulates drought tolerance in self-rooted apple stocks. Altogether, we believe that our findings provide novel insights into the functions of MdFLP in the regulation of drought tolerance in self-rooted apple stocks.

Comprehensive analysis of 1R- and 2R-MYBs reveals novel genic and protein features, complex organisation, selective expansion and insights into evolutionary tendencies

Myeloblastosis (MYB) family, the largest plant transcription factor family, has been subcategorised based on the number and type of repeats in the MYB domain. In spite of several reports, evolution of MYB genes and repeats remains enigmatic. Brassicaceae members are endowed with complex genomes, including dysploidy because of its unique history with multiple rounds of polyploidisation, genomic fractionations and rearrangements. The present study is an attempt to gain insights into the complexities of MYB family diversity, understand impacts of genome evolution on gene families and develop an evolutionary framework to understand the origin of various subcategories of MYB gene family. We identified and analysed 1129 MYBs that included 1R-, 2R-, 3R- and atypical-MYBs across sixteen species representing protists, fungi, animals and plants and exclude MYB identified from Brassicaceae except Arabidopsis thaliana; in addition, a total of 1137 2R-MYB genes from six Brassicaceae species were also analysed. Comparative analysis revealed predominance of 1R-MYBs in protists, fungi, animals and lower plants. Phylogenetic reconstruction and analysis of selection pressure suggested ancestral nature of R1-type repeat containing 1R-MYBs that might have undergone intragenic duplication to form multi-repeat MYBs. Distinct differences in gene structure between 1R-MYB and 2R-MYBs were observed regarding intron number, the ratio of gene length to coding DNA sequence (CDS) length and the length of exons encoding the MYB domain. Conserved as well as novel and lineage-specific intron phases were identified. Analyses of physicochemical properties revealed drastic differences indicating functional diversification in MYBs. Phylogenetic reconstruction of 1R- and 2R-MYB genes revealed a shared structure-function relationship in clades which was supported when transcriptome data was analysed in silico. Comparative genomics to study distribution pattern and mapping of 2R-MYBs revealed congruency and greater degree of synteny and collinearity among closely related species. Micro-synteny analysis of genomic segments revealed high conservation of genes that are immediately flanking the surrounding tandemly organised 2R-MYBs along with instances of local duplication, reorganisations and genome fractionation. In summary, polyploidy, dysploidy, reshuffling and genome fractionation were found to cause loss or gain of 2R-MYB genes. The findings need to be supported with functional validation to understand gene structure-function relationship along the evolutionary lineage and adaptive strategies based on comparative functional genomics in plants.

miR390 family of Cymbidium goeringii is involved in the development of reproductive organs in transgenic Arabidopsis

BACKGROUND

miR390s is an ancient family with a high level of conservation among plant miRNAs. Through the auxin signal transduction pathway, miR390 participates in diverse biological processes of plant growth and development. As an important Chinese traditional orchid, Cymbidium goeringii has unique flower shape and elegant fragrance. But its development has been greatly restricted because of the low flower bud differentiation and the difficult reproduction. This study aims to provide guidance for the role of cgo-miR390 in reproductive organ development to enhance the ornamental and economic value of Cymbidium.

RESULTS

MIR390a, MIR390b and MIR390c of C. goeringii were cloned, and their length ranged from 130 to 150 nt. Each precursor sequence of cgo-miR390 contains 2 to 3 mature miRNAs. Three kinds of cgo-miR390s displayed distinct temporal and spatial expression patterns during floral development in C. goeringii. The overexpression of MIR390s alters morphology and function of stamens and pistils in Arabidopsis, such as enlargement of anther aspect ratio and separation of stylar and stigmas, which affects the development of fruits and seeds. In particular, the pollen amount decreased and the seed abortion rate increased in cgo-MIR390c-overexpressed plants.

CONCLUSIONS

cgo-miR390 family affected the development of reproductive organs in transgenic Arabidopsis. The study provides references for the genetic improvement for orchid with potentially great economic benefit.

Arabidopsis F-BOX STRESS INDUCED 4 is required to repress excessive divisions in stomatal development

During the terminal stage of stomatal development, the R2R3-MYB transcription factors FOUR LIPS (FLP/MYB124) and MYB88 limit guard mother cell division by repressing the transcript levels of multiple cell-cycle genes. In Arabidopsis thaliana possessing the weak allele flp-1, an extra guard mother cell division results in two stomata having direct contact. Here, we identified an ethylmethane sulfonate-mutagenized mutant, flp-1 xs01c, which exhibited more severe defects than flp-1 alone, producing giant tumor-like cell clusters. XS01C, encoding F-BOX STRESS-INDUCED 4 (FBS4), is preferentially expressed in epidermal stomatal precursor cells. Overexpressing FBS4 rescued the defective stomatal phenotypes of flp-1 xs01c and flp-1 mutants. The deletion or substitution of a conserved residue (Proline166) within the F-box domain of FBS4 abolished or reduced, respectively, its interaction with Arabidopsis Skp1-Like1 (ASK1), the core subunit of the Skp1/Cullin/F-box E3 ubiquitin ligase complex. Furthermore, the FBS4 protein physically interacted with CYCA2;3 and induced its degradation through the ubiquitin-26S proteasome pathway. Thus, in addition to the known transcriptional pathway, the terminal symmetric division in stomatal development is ensured at the post-translational level, such as through the ubiquitination of target proteins recognized by the stomatal lineage F-box protein FBS4.

Capsicum chinense MYB Transcription Factor Genes: Identification, Expression Analysis, and Their Conservation and Diversification With Other Solanaceae Genomes

Myeloblastosis (MYB) genes are important transcriptional regulators of plant growth, development, and secondary metabolic biosynthesis pathways, such as capsaicinoid biosynthesis in Capsicum. Although MYB genes have been identified in Capsicum annuum, no comprehensive study has been conducted on other Capsicum species. We identified a total of 251 and 240 MYB encoding genes in Capsicum chinense MYBs (CcMYBs) and Capsicum baccatum MYBs (CbMYBs). The observation of twenty tandem and 41 segmental duplication events indicated expansion of the MYB gene family in the C. chinense genome. Five CcMYB genes, i.e., CcMYB101, CcMYB46, CcMYB6, CcPHR8, and CcRVE5, and two CaMYBs, i.e., CaMYB3 and CaHHO1, were found within the previously reported capsaicinoid biosynthesis quantitative trait loci. Based on phylogenetic analysis with tomato MYB proteins, the Capsicum MYBs were classified into 24 subgroups supported by conserved amino acid motifs and gene structures. Also, a total of 241 CcMYBs were homologous with 225 C. annuum, 213 C. baccatum, 125 potato, 79 tomato, and 23 Arabidopsis MYBs. Synteny analysis showed that all 251 CcMYBs were collinear with C. annuum, C. baccatum, tomato, potato, and Arabidopsis MYBs spanning over 717 conserved syntenic segments. Using transcriptome data from three fruit developmental stages, a total of 54 CcMYBs and 81 CaMYBs showed significant differential expression patterns. Furthermore, the expression of 24 CcMYBs from the transcriptome data was validated by quantitative real-time (qRT) PCR analysis. Eight out of the 24 CcMYBs validated by the qRT-PCR were highly expressed in fiery hot C. chinense than in the lowly pungent C. annuum. Furthermore, the co-expression analysis revealed several MYB genes clustered with genes from the capsaicinoid, anthocyanin, phenylpropanoid, carotenoid, and flavonoids biosynthesis pathways, and related to determining fruit shape and size. The homology modeling of 126 R2R3 CcMYBs showed high similarity with that of the Arabidopsis R2R3 MYB domain template, suggesting their potential functional similarity at the proteome level. Furthermore, we have identified simple sequence repeat (SSR) motifs in the CcMYB genes, which could be used in Capsicum breeding programs. The functional roles of the identified CcMYBs could be studied further so that they can be manipulated for Capsicum trait improvement.

The regulatory module MdBT2-MdMYB88/MdMYB124-MdNRTs regulates nitrogen usage in apple

Less than 40% of the nitrogen (N) fertilizer applied to soil is absorbed by crops. Thus, improving the N use efficiency of crops is critical for agricultural development. However, the underlying regulation of these processes remains largely unknown, particularly in woody plants. By conducting yeast two-hybrid assays, we identified one interacting protein of MdMYB88 and MdMYB124 in apple (Malus × domestica), namely BTB and TAZ domain protein 2 (MdBT2). Ubiquitination and protein stabilization analysis revealed that MdBT2 ubiquitinates and degrades MdMYB88 and MdMYB124 via the 26S proteasome pathway. MdBT2 negatively regulates nitrogen usage as revealed by the reduced fresh weight, dry weight, N concentration, and N usage index of MdBT2 overexpression calli under low-N conditions. In contrast, MdMYB88 and MdMYB124 increase nitrate absorption, allocation, and remobilization by regulating expression of MdNRT2.4, MdNRT1.8, MdNRT1.7, and MdNRT1.5 under N limitation, thereby regulating N usage. The results obtained illustrate the mechanism of a regulatory module comprising MdBT2-MdMYB88/MdMYB124-MdNRTs, through which plants modulate N usage. These data contribute to a molecular approach to improve the N usage of fruit crops under limited N acquisition.

Mating type specific transcriptomic response to sex inducing pheromone in the pennate diatom Seminavis robusta

Sexual reproduction is a fundamental phase in the life cycle of most diatoms. Despite its role as a source of genetic variation, it is rarely reported in natural circumstances and its molecular foundations remain largely unknown. Here, we integrate independent transcriptomic datasets to prioritize genes responding to sex inducing pheromones (SIPs) in the pennate diatom Seminavis robusta. We observe marked gene expression changes associated with SIP treatment in both mating types, including an inhibition of S phase progression, chloroplast division, mitosis, and cell wall formation. Meanwhile, meiotic genes are upregulated in response to SIP, including a sexually induced diatom specific cyclin. Our data further suggest an important role for reactive oxygen species, energy metabolism, and cGMP signaling during the early stages of sexual reproduction. In addition, we identify several genes with a mating type specific response to SIP, and link their expression pattern with physiological specialization, such as the production of the attraction pheromone diproline in mating type − (MT−) and mate-searching behavior in mating type + (MT+). Combined, our results provide a model for early sexual reproduction in pennate diatoms and significantly expand the suite of target genes to detect sexual reproduction events in natural diatom populations.

Expansin gene TaEXPA2 positively regulates drought tolerance in transgenic wheat (Triticum aestivum L.)

Expansins loosen plant cell walls and are involved in cell enlargement and various abiotic stresses. In previous studies, we cloned the expansin gene TaEXPA2 from the wheat cultivar HF9703. Here, we studied its function and regulation in wheat drought stress tolerance. The results indicated that TaEXPA2-overexpressing wheat plants (OE) exhibited drought tolerant phenotypes, whereas down-regulation of TaEXPA2 by RNA interference (RNAi) resulted in elevated drought sensitivity, as measured by survival rate, photosynthetic rate and water containing ability under drought stress. Overexpression of TaEXPA2 enhanced the antioxidant capacity in wheat plants, via elevation of antioxidant enzyme activity and the increase of the transcripts of some ROS scavenging enzyme-related genes. Further investigation revealed that TaEXPA2 positively influenced lateral root formation under drought conditions. A MYB transcription factor of wheat named TaMPS activates TaEXPA2 expression directly by binding to its promoter. Overexpression of TaMPS in Arabidopsis conferred drought tolerance associated with improved lateral root number, and the close homolog genes of TaEXPA2 were up-regulated in Arabidopsis roots overexpressing TaMPS, which suggest that TaMPS may function as one of the regulator of TaEXPA2 gene expression in the root lateral development under drought stress. These findings suggest that TaEXPA2 positively regulates drought stress tolerance in wheat.

Phylogenomic Analysis of R2R3 MYB Transcription Factors in Sorghum and their Role in Conditioning Biofuel Syndrome

Background

Large scale cultivation of sorghum for food, feed, and biofuel requires concerted efforts for engineering multipurpose cultivars with optimised agronomic traits. Due to their vital role in regulating the biosynthesis of phenylpropanoid-derived compounds, biomass composition, biotic, and abiotic stress response, R2R3-MYB family transcription factors are ideal targets for improving environmental resilience and economic value of sorghum.

Methods

We used diverse computational biology tools to survey the sorghum genome to identify R2R3-MYB transcription factors followed by their structural and phylogenomic analysis. We used in-house generated as well as publicly available high throughput expression data to analyse the R2R3 expression patterns in various sorghum tissue types.

Results

We have identified a total of 134 R2R3-MYB genes from sorghum and developed a framework to predict gene functions. Collating information from the physical location, duplication, structural analysis, orthologous sequences, phylogeny, and expression patterns revealed the role of duplications in clade-wise expansion of the R2R3-MYB family as well as intra-clade functional diversification. Using publicly available and in-house generated RNA sequencing data, we provide MYB candidates for conditioning biofuel syndrome by engineering phenylpropanoid biosynthesis and sugar signalling pathways in sorghum.

Conclusion

The results presented here are pivotal to prioritize MYB genes for functional validation and optimize agronomic traits in sorghum.

Molecular Mechanisms of the Floral Biology of Jatropha curcas: Opportunities and Challenges as an Energy Crop

Fossil fuel sources are a limited resource and could eventually be depleted. Biofuels have emerged as a renewable alternative to fossil fuels. Jatropha has grown in significance as a potential bioenergy crop due to its high content of seed oil. However, Jatropha's lack of high-yielding seed genotypes limits its potential use for biofuel production. The main cause of lower seed yield is the low female to male flower ratio (1:25-10), which affects the total amount of seeds produced per plant. Here, we review the genetic factors responsible for floral transitions, floral organ development, and regulated gene products in Jatropha. We also summarize potential gene targets to increase seed production and discuss challenges ahead.

Genome-Wide Association Analysis Identifies Candidate Genes Regulating Seed Number Per Silique in Arabidopsis thaliana

Seed weight and number ultimately determine seed yield. Arabidopsis seed number comprised of silique number and seed number per silique (SNS). Comparing seed development and weight, determinants of seed number remain largely uncharacterized. In this study, taking advantage of 107 available Arabidopsis accessions, genome-wide association analysis (GWAS) was employed to identify the candidate genes regulating SNS. GWAS-based genotype and phenotype association analysis identified 38 most significant SNPs marker sites that were mapped to specific chromosomal positions and allowed us to screen for dozens of candidate genes. One of them (PIN3) was selected for functional validation based on gene expression analysis. It is a positive regulator of Arabidopsis SNS. Although silique length of PIN3 loss of function mutant was not significantly changed, its SNS and seed density (SD) were significantly reduced as compared with the wild type. Notably, PIN3 overexpression lines driven by a placenta-specific promoter STK exhibited significantly shorter siliques, slightly reduced SNS, but significant increased SD compared with wild type, suggesting that PIN3 positively regulates SD through inducing ovule primordia initiation regardless of the placenta size. Ovule initiation determines the maximal possibility of SNS, and new genes and mechanism regulating SNS through modulating ovule initiation is worth further investigated.

Comparative gene expression profile analysis of ovules provides insights into Jatropha curcas L. ovule development

Jatropha curcas, an economically important biofuel feedstock with oil-rich seeds, has attracted considerable attention among researchers in recent years. Nevertheless, valuable information on the yield component of this plant, particularly regarding ovule development, remains scarce. In this study, transcriptome profiles of anther and ovule development were established to investigate the ovule development mechanism of J. curcas. In total, 64,325 unigenes with annotation were obtained, and 1723 differentially expressed genes (DEGs) were identified between different stages. The DEG analysis showed the participation of five transcription factor families (bHLH, WRKY, MYB, NAC and ERF), five hormone signaling pathways (auxin, gibberellic acid (GA), cytokinin, brassinosteroids (BR) and jasmonic acid (JA)), five MADS-box genes (AGAMOUS-2, AGAMOUS-1, AGL1, AGL11, and AGL14), SUP and SLK3 in ovule development. The role of GA and JA in ovule development was evident with increases in flower buds during ovule development: GA was increased approximately twofold, and JA was increased approximately sevenfold. In addition, the expression pattern analysis using qRT-PCR revealed that CRABS CLAW and AGAMOUS-2 were also involved in ovule development. The upregulation of BR signaling genes during ovule development might have been regulated by other phytohormone signaling pathways through crosstalk. This study provides a valuable framework for investigating the regulatory networks of ovule development in J. curcas.

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.

Transcriptional regulation of MdPIN3 and MdPIN10 by MdFLP during apple self-rooted stock adventitious root gravitropism

BACKGROUND

The close planting of dwarfing self-rooted rootstocks is currently a widely used method for apple production; however, self-rooted rootstocks are weak with shallow roots and poor grounding. Therefore, understanding the molecular mechanisms that establish the gravitropic set-point angles (GSAs) of the adventitious roots of self-rooted apple stocks is important for developing self-rooted apple rootstock cultivars with deep roots.

RESULTS

We report that the apple FOUR LIPS (MdFLP), an R2R3-MYB transcription factor (TF), functions in establishing the GSA of the adventitious roots of self-rooted apple stocks in response to gravity. Biochemical analyses demonstrate that MdFLP directly binds to the promoters of two auxin efflux carriers, MdPIN3 and MdPIN10, that are involved in auxin transport, activates their transcriptional expression, and thereby promotes the development of adventitious roots in self-rooted apple stocks. Additionally, the apple auxin response factor MdARF19 influences the expression of those auxin efflux carriers and the establishment of the GSA of adventitious roots of apple in response to gravity by directly activating the expression of MdFLP.

CONCLUSION

Our findings provide new insights into the transcriptional regulation of MdFLP by the auxin response factor MdARF19 in the regulation of the GSA of adventitious roots of self-rooted apple stocks in response to gravity.

Identification of lncRNAs involved in rice ovule development and female gametophyte abortion by genome-wide screening and functional analysis

Background

As important female reproductive tissues, the rice (Oryza sativa L.) ovule and female gametophyte is significant in terms of their fertility. Long noncoding RNAs (lncRNAs) play important and wide-ranging roles in the growth and development of plants and have become a major research focus in recent years. Therefore, we explored the characterization and expression change of lncRNAs during ovule development and female gametophytic abortion.

Results

In our study, whole-transcriptome strand-specific RNA sequencing (ssRNA-seq) was performed in the ovules of a high-frequency female-sterile rice line (fsv1) and a wild-type rice line (Gui99) at the megaspore mother cell meiosis stage (stage 1), functional megaspore mitosis stage (stage 2) and female gametophyte mature stage (stage 3). By comparing two rice lines, we identified 152, 233, and 197 differentially expressed lncRNAs at the three ovule developmental stages. Functional analysis of the coherent target genes of these differentially expressed lncRNAs indicated that many lncRNAs participate in multiple pathways such as hormone and cellular metabolism and signal transduction. Moreover, there were many differentially expressed lncRNAs acting as the precursors of some miRNAs that are involved in the development of ovules and female gametophytes. In addition, we have found that lncRNAs can act as decoys, competing with mRNAs for binding to miRNAs to maintain the normal expression of genes related to ovule and female gametophyte development.

Conclusion

These results provide important clues for elucidating the female gametophyte abortion mechanism in rice. This study also expands our understanding about the biological functions of lncRNAs and the annotation of the rice genome.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5442-6) contains supplementary material, which is available to authorized users.

MdMYB88 and MdMYB124 Enhance Drought Tolerance by Modulating Root Vessels and Cell Walls in Apple

Water deficit is one of the main limiting factors in apple (Malus × domestica Borkh.) cultivation. Root architecture plays an important role in the drought tolerance of plants; however, research efforts to improve drought tolerance of apple trees have focused on aboveground targets. Due to the difficulties associated with visualization and data analysis, there is currently a poor understanding of the genetic players and molecular mechanisms involved in the root architecture of apple trees under drought conditions. We previously observed that MdMYB88 and its paralog MdMYB124 regulate apple tree root morphology. In this study, we found that MdMYB88 and MdMYB124 play important roles in maintaining root hydraulic conductivity under long-term drought conditions and therefore contribute toward adaptive drought tolerance. Further investigation revealed that MdMYB88 and MdMYB124 regulate root xylem development by directly binding MdVND6 and MdMYB46 promoters and thus influence expression of their target genes under drought conditions. In addition, MdMYB88 and MdMYB124 were shown to regulate the deposition of cellulose and lignin root cell walls in response to drought. Taken together, our results provide novel insights into the importance of MdMYB88 and MdMYB124 in root architecture, root xylem development, and secondary cell wall deposition in response to drought in apple trees.

Ectopic Overexpression of bol-miR171b Increases Chlorophyll Content and Results in Sterility in Broccoli ( Brassica oleracea L var. italica)

MiR171 plays pleiotropic roles in the growth and development of several plant species. However, the mechanism underlying the miR171-mediated regulation of organ development in broccoli remains unknown. In this study, bol-miR171b was characterized and found to be differentially expressed in various broccoli organs. The ectopic overexpression of bol-miR171b in Arabidopsis affected the leaf and silique development of transgenic lines. In particular, the chlorophyll content of leaves from overexpressed bol-miR171b transgenic Arabidopsis was higher than that of the vector controls. The fertility and seed yield of Arabidopsis with overexpressed bol-miR171b were markedly lower than those of the vector controls. Similarly, overexpressed bol-miR171b transgenic broccoli exhibited dark green leaves with high chlorophyll content, and nearly all of the flowers were sterile. These results demonstrated that overexpression of bol-miR171b could increase the chlorophyll content of transgenic plants. Degradome sequencing was conducted to identify the targets of bol-miR171b. Two members of the GRAS gene family, BolSCL6 and BolSCL27, were cleaved by bol-miR171b-3p in broccoli. In addition to the genes targeted by bol-miR171b-3p, adenylylsulfate reductase 3 ( APSR3), which played important roles in plant sulfate assimilation and reduction, was speculated to be cleaved by bol-miR171b-5p, suggesting that the star sequence of bol-miR171b may also have functions in broccoli. Comparative transcriptome analysis further revealed that the genes involved in chloroplast development and sulfate homeostasis should participate in the bol-miR171b -mediated regulatory network. Taken together, these findings provided new insights into the function and regulation of bol-miR171b in broccoli and indicated the potential of bol-miR171b as a small RNA molecule that increased leaf chlorophyll in plants by genetic engineering.

RNA-seq Analysis Reveals Gene Expression Profiling of Female Fertile and Sterile Ovules of PinusTabulaeformis Carr. during Free Nuclear Mitosis of the Female Gametophyte

The development of the female gametophyte (FG) is one of the key processes of life cycle alteration between the haploid gametophyte and the diploid sporophytes in plants and it is required for successful seed development after fertilization. It is well demonstrated that free nuclear mitosis (FNM) of FG is crucial for the development of the ovule. However, studies of the molecular mechanism of ovule and FG development focused mainly on angiosperms, such as Arabidopsis thaliana and further investigation of gymnosperms remains to be completed. Here, Illumina sequencing of six transcriptomic libraries obtained from developing and abortive ovules at different stages during free nuclear mitosis of magagametophyte (FNMM) was used to acquire transcriptome data and gene expression profiles of Pinus tabulaeformis. Six cDNA libraries generated a total of 71.0 million high-quality clean reads that aligned with 63,449 unigenes and the comparison between developing and abortive ovules identified 7174 differentially expressed genes (DEGs). From the functional annotation results, DEGs involved in the cell cycle and phytohormone regulation were highlighted to reveal their biological importance in ovule development. Furthermore, validation of DEGs from the phytohormone signal transduction pathway was performed using quantitative real-time PCR analysis, revealing the dynamics of transcriptional networks and potential key components in the regulation of FG development in P. tabulaeformis were identified. These findings provide new insights into the regulatory mechanisms of ovule development in woody gymnosperms.

An atypical R2R3 MYB transcription factor increases cold hardiness by CBF-dependent and CBF-independent pathways in apple

Apple (Malus × domestica) trees are vulnerable to freezing temperatures. However, there has been only limited success in developing cold-hardy cultivars. This lack of progress is due at least partly to lack of understanding of the molecular mechanisms of freezing tolerance in apple. In this study, we evaluated the potential roles for two R2R3 MYB transcription factors (TFs), MYB88 and the paralogous FLP (MYB124), in cold stress in apple and Arabidopsis. We found that MYB88 and MYB124 positively regulate freezing tolerance and cold-responsive gene expression in both apple and Arabidopsis. Chromatin-Immunoprecipitation-qPCR and electrophoretic mobility shift assays showed that MdMYB88/MdMYB124 act as direct regulators of the COLD SHOCK DOMAIN PROTEIN 3 (MdCSP3) and CIRCADIAN CLOCK ASSOCIATED 1 (MdCCA1) genes. Dual luciferase reporter assay indicated that MdCCA1 but not MdCSP3 activated the expression of MdCBF3 under cold stress. Moreover, MdMYB88 and MdMYB124 promoted anthocyanin accumulation and H₂ O₂ detoxification in response to cold. Taken together, our results suggest that MdMYB88 and MdMYB124 positively regulate cold hardiness and cold-responsive gene expression under cold stress by C-REPEAT BINDING FACTOR (CBF)-dependent and CBF-independent pathways.

GmMYB181, a Soybean R2R3-MYB Protein, Increases Branch Number in Transgenic Arabidopsis

Soybean (Glycine max) is an important economic crop that provides abundant oil and high quality protein for human beings. As the process of reproductive growth directly determines the crop seed yield and quality, we initiated studies to identify genes that regulate soybean floral organ development. One R2R3-MYB transcription factor gene, designated as GmMYB181, was found to be enriched in flowers based on microarray analysis and was further functionally investigated in transgenic Arabidopsis. GmMYB181 protein contains two MYB domains, which localized to the nucleus and displayed transcriptional activation in yeast hybrid system. Real-time quantitative PCR (qRT-PCR) results suggested GmMYB181 exclusively expressed in flower tissue. In Arabidopsis, overexpression of GmMYB181 altered the morphology of floral organs, fruit size and plant architecture, including outward curly sepals, smaller siliques, increased lateral branches and reduced plant height, indicating that GmMYB181 is involved in the development of reproductive organs and plays an important role in controlling plant architecture. Further, microarray analysis revealed that overexpressing GmMYB181 in Arabidopsis affected the expression of 3450 genes in mature flowers, including those involved in floral organ, seed/fruit development, and responded to different hormone signals.

Genome-Wide Analysis of DNA Methylation During Ovule Development of Female-Sterile Rice fsv1

The regulation of female fertility is an important field of rice sexual reproduction research. DNA methylation is an essential epigenetic modification that dynamically regulates gene expression during development processes. However, few reports have described the methylation profiles of female-sterile rice during ovule development. In this study, ovules were continuously acquired from the beginning of megaspore mother cell meiosis until the mature female gametophyte formation period, and global DNA methylation patterns were compared in the ovules of a high-frequency female-sterile line (fsv1) and a wild-type rice line (Gui99) using whole-genome bisulfite sequencing (WGBS). Profiling of the global DNA methylation revealed hypo-methylation, and 3471 significantly differentially methylated regions (DMRs) were observed in fsv1 ovules compared with Gui99. Based on functional annotation and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis of differentially methylated genes (DMGs), we observed more DMGs enriched in cellular component, reproduction regulation, metabolic pathway, and other pathways. In particular, many ovule development genes and plant hormone-related genes showed significantly different methylation patterns in the two rice lines, and these differences may provide important clues for revealing the mechanism of female gametophyte abortion.

De novo comparative transcriptome analysis of genes involved in fruit morphology of pumpkin cultivars with extreme size difference and development of EST-SSR markers

The genetic basis of fruit size and shape was investigated for the first time in Cucurbita species and genetic loci associated with fruit morphology have been identified. Although extensive genomic resources are available at present for tomato (Solanum lycopersicum), cucumber (Cucumis sativus), melon (Cucumis melo) and watermelon (Citrullus lanatus), genomic databases for Cucurbita species are limited. Recently, our group reported the generation of pumpkin (Cucurbita pepo) transcriptome databases from two contrasting cultivars with extreme fruit sizes. In the current study we used these databases to perform comparative transcriptome analysis in order to identify genes with potential roles in fruit morphology and fruit size. Differential Gene Expression (DGE) analysis between cv. 'Munchkin' (small-fruit) and cv. 'Big Moose' (large-fruit) revealed a variety of candidate genes associated with fruit morphology with significant differences in gene expression between the two cultivars. In addition, we have set the framework for generating EST-SSR markers, which discriminate different C. pepo cultivars and show transferability to related Cucurbitaceae species. The results of the present study will contribute to both further understanding the molecular mechanisms regulating fruit morphology and furthermore identifying the factors that determine fruit size. Moreover, they may lead to the development of molecular marker tools for selecting genotypes with desired morphological traits.

Transcriptional regulation of PIN genes by FOUR LIPS and MYB88 during Arabidopsis root gravitropism

PIN proteins are auxin export carriers that direct intercellular auxin flow and in turn regulate many aspects of plant growth and development including responses to environmental changes. The Arabidopsis R2R3-MYB transcription factor FOUR LIPS (FLP) and its paralogue MYB88 regulate terminal divisions during stomatal development, as well as female reproductive development and stress responses. Here we show that FLP and MYB88 act redundantly but differentially in regulating the transcription of PIN3 and PIN7 in gravity-sensing cells of primary and lateral roots. On the one hand, FLP is involved in responses to gravity stimulation in primary roots, whereas on the other, FLP and MYB88 function complementarily in establishing the gravitropic set-point angles of lateral roots. Our results support a model in which FLP and MYB88 expression specifically determines the temporal-spatial patterns of PIN3 and PIN7 transcription that are closely associated with their preferential functions during root responses to gravity.

Plants respond to reorientation by altering the distribution of the plant hormone auxin causing roots to bend towards gravity. Here Wang et al. find that expression of the PIN3 and PIN7 auxin transporters in gravity sensing cells is controlled by concerted action of the FLP and MYB88 transcription factors.

The FOUR LIPS and MYB88 transcription factor genes are widely expressed in Arabidopsis thaliana during development

PREMISE OF THE STUDY

The FOUR LIPS (FLP) and MYB88 transcription factors, which are closely related in structure and function, control the development of stomata, as well as entry into megasporogenesis in Arabidopsis thaliana. However, other locations where these transcription factors are expressed are poorly described. Documenting additional locations where these genes are expressed might define new functions for these genes.

METHODS

Expression patterns were examined throughout vegetative and reproductive development. The expression from two transcriptional-reporter fusions were visualized with either β-glucuronidase (GUS) or green fluorescence protein (GFP).

KEY RESULTS

Both flp and myb88 genes were expressed in many, previously unreported locations, consistent with the possibility of additional functions for FLP and MYB88. Moreover, expression domains especially of FLP display sharp cutoffs or boundaries.

CONCLUSIONS

In addition to stomatal and reproductive development, FLP and MYB88, which are R2R3 MYB transcription factor genes, are expressed in many locations in cells, tissues, and organs.

Transcriptomic analysis reveals the roles of microtubule-related genes and transcription factors in fruit length regulation in cucumber (Cucumis sativus L.)

Cucumber (Cucumis sativus L.) fruit is a type of fleshy fruit that is harvested immaturely. Early fruit development directly determines the final fruit length and diameter, and consequently the fruit yield and quality. Different cucumber varieties display huge variations of fruit length, but how fruit length is determined at the molecular level remains poorly understood. To understand the genes and gene networks that regulate fruit length in cucumber, high throughout RNA-Seq data were used to compare the transcriptomes of early fruit from two near isogenic lines with different fruit lengths. 3955 genes were found to be differentially expressed, among which 2368 genes were significantly up-regulated and 1587 down-regulated in the line with long fruit. Microtubule and cell cycle related genes were dramatically activated in the long fruit, and transcription factors were implicated in the fruit length regulation in cucumber. Thus, our results built a foundation for dissecting the molecular mechanism of fruit length control in cucumber, a key agricultural trait of significant economic importance.

MULTIPASS, a rice R2R3-type MYB transcription factor, regulates adaptive growth by integrating multiple hormonal pathways

Growth regulation is an important aspect of plant adaptation during environmental perturbations. Here, the role of MULTIPASS (OsMPS), an R2R3-type MYB transcription factor of rice, was explored. OsMPS is induced by salt stress and expressed in vegetative and reproductive tissues. Over-expression of OsMPS reduces growth under non-stress conditions, while knockdown plants display increased biomass. OsMPS expression is induced by abscisic acid and cytokinin, but is repressed by auxin, gibberellin and brassinolide. Growth retardation caused by OsMPS over-expression is partially restored by auxin application. Expression profiling revealed that OsMPS negatively regulates the expression of EXPANSIN (EXP) and cell-wall biosynthesis as well as phytohormone signaling genes. Furthermore, the expression of OsMPS-dependent genes is regulated by auxin, cytokinin and abscisic acid. Moreover, we show that OsMPS is a direct upstream regulator of OsEXPA4, OsEXPA8, OsEXPB2, OsEXPB3, OsEXPB6 and the endoglucanase genes OsGLU5 and OsGLU14. The multiple responses of OsMPS and its target genes to various hormones suggest an integrative function of OsMPS in the cross-talk between phytohormones and the environment to regulate adaptive growth.

Transcript profile analyses of maize silks reveal effective activation of genes involved in microtubule-based movement, ubiquitin-dependent protein degradation, and transport in the pollination process

Pollination is the first crucial step of sexual reproduction in flowering plants, and it requires communication and coordination between the pollen and the stigma. Maize (Zea mays) is a model monocot with extraordinarily long silks, and a fully sequenced genome, but little is known about the mechanism of its pollen-stigma interactions. In this study, the dynamic gene expression of silks at four different stages before and after pollination was analyzed. The expression profiles of immature silks (IMS), mature silks (MS), and silks at 20 minutes and 3 hours after pollination (20MAP and 3HAP, respectively) were compared. In total, we identified 6,337 differentially expressed genes in silks (SDEG) at the four stages. Among them, the expression of 172 genes were induced upon pollination, most of which participated in RNA binding, processing and transcription, signal transduction, and lipid metabolism processes. Genes in the SDEG dataset could be divided into 12 time-course clusters according to their expression patterns. Gene Ontology (GO) enrichment analysis revealed that many genes involved in microtubule-based movement, ubiquitin-mediated protein degradation, and transport were predominantly expressed at specific stages, indicating that they might play important roles in the pollination process of maize. These results add to current knowledge about the pollination process of grasses and provide a foundation for future studies on key genes involved in the pollen-silk interaction in maize.