MADS-complexes regulate transcriptome dynamics during pollen maturation

Whole-Genome Sequencing and Analysis of Tumour-Forming Radish (Raphanus sativus L.) Line

Spontaneous tumour formation in higher plants can occur in the absence of pathogen invasion, depending on the plant genotype. Spontaneous tumour formation on the taproots is consistently observed in certain inbred lines of radish (Raphanus sativus var. radicula Pers.). In this paper, using Oxford Nanopore and Illumina technologies, we have sequenced the genomes of two closely related radish inbred lines that differ in their ability to spontaneously form tumours. We identified a large number of single nucleotide variants (amino acid substitutions, insertions or deletions, SNVs) that are likely to be associated with the spontaneous tumour formation. Among the genes involved in the trait, we have identified those that regulate the cell cycle, meristem activity, gene expression, and metabolism and signalling of phytohormones. After identifying the SNVs, we performed Sanger sequencing of amplicons corresponding to SNV-containing regions to validate our results. We then checked for the presence of SNVs in other tumour lines of the radish genetic collection and found the ERF118 gene, which had the SNVs in the majority of tumour lines. Furthermore, we performed the identification of the CLAVATA3/ESR (CLE) and WUSCHEL (WOX) genes and, as a result, identified two unique radish CLE genes which probably encode proteins with multiple CLE domains. The results obtained provide a basis for investigating the mechanisms of plant tumour formation and also for future genetic and genomic studies of radish.

Genome-wide analysis of MADS-box transcription factor gene family in wild emmer wheat (Triticum turgidum subsp. dicoccoides)

The members of MADS-box gene family have important roles in regulating the growth and development of plants. MADS-box genes are highly regarded for their potential to enhance grain yield and quality under shifting global conditions. Wild emmer wheat (Triticum turgidum subsp. dicoccoides) is a progenitor of common wheat and harbors valuable traits for wheat improvement. Here, a total of 117 MADS-box genes were identified in the wild emmer wheat genome and classified to 90 MIKCC, 3 MIKC*, and 24 M-type. Furthermore, a phylogenetic analysis and expression profiling of the emmer wheat MADS-box gene family was presented. Although some MADS-box genes belonging to SOC1, SEP1, AGL17, and FLC groups have been expanded in wild emmer wheat, the number of MIKC-type MADS-box genes per subgenome is similar to that of rice and Arabidopsis. On the other hand, M-type genes of wild emmer wheat is less frequent than that of Arabidopsis. Gene expression patterns over different tissues and developmental stages agreed with the subfamily classification of MADS-box genes and was similar to common wheat and rice, indicating their conserved functionality. Some TdMADS-box genes are also differentially expressed under drought stress. The promoter region of each of the TdMADS-box genes harbored 6 to 48 responsive elements, mainly related to light, however hormone, drought, and low-temperature related cis-acting elements were also present. In conclusion, the results provide detailed information about the MADS-box genes of wild emmer wheat. The present work could be useful in the functional genomics efforts toward breeding for agronomically important traits in T. dicoccoides.

Transcription Factors and Their Regulatory Roles in the Male Gametophyte Development of Flowering Plants

Male gametophyte development in plants relies on the functions of numerous genes, whose expression is regulated by transcription factors (TFs), non-coding RNAs, hormones, and diverse environmental stresses. Several excellent reviews are available that address the genes and enzymes associated with male gametophyte development, especially pollen wall formation. Growing evidence from genetic studies, transcriptome analysis, and gene-by-gene studies suggests that TFs coordinate with epigenetic machinery to regulate the expression of these genes and enzymes for the sequential male gametophyte development. However, very little summarization has been performed to comprehensively review their intricate regulatory roles and discuss their downstream targets and upstream regulators in this unique process. In the present review, we highlight the research progress on the regulatory roles of TF families in the male gametophyte development of flowering plants. The transcriptional regulation, epigenetic control, and other regulators of TFs involved in male gametophyte development are also addressed.

Cell-type-specific alternative splicing in the Arabidopsis germline

During sexual reproduction in flowering plants, the two haploid sperm cells (SCs) embedded within the cytoplasm of a growing pollen tube are carried to the embryo sac for double fertilization. Pollen development in flowering plants is a dynamic process that encompasses changes at transcriptome and epigenome levels. While the transcriptome of pollen and SCs in Arabidopsis (Arabidopsis thaliana) is well documented, previous analyses have mostly been based on gene-level expression. In-depth transcriptome analysis, particularly the extent of alternative splicing (AS) at the resolution of SC and vegetative nucleus (VN), is still lacking. Therefore, we performed RNA-seq analysis to generate a spliceome map of Arabidopsis SCs and VN isolated from mature pollen grains. Based on our de novo transcriptome assembly, we identified 58,039 transcripts, including 9,681 novel transcripts, of which 2,091 were expressed in SCs and 3,600 in VN. Four hundred and sixty-eight genes were regulated both at gene and splicing levels, with many having functions in mRNA splicing, chromatin modification, and protein localization. Moreover, a comparison with egg cell RNA-seq data uncovered sex-specific regulation of transcription and splicing factors. Our study provides insights into a gamete-specific AS landscape at unprecedented resolution.

Transcriptome-wide identification of WRKY transcription factors and their expression profiles under different stress in Cynanchum thesioides

Cynanchum thesioides (Freyn) K. Schum. is an important economic and medicinal plant widely distributed in northern China. WRKY transcription factors (TFs) play important roles in plant growth, development and regulating responses. However, there is no report on the WRKY genes in Cynanchum thesioides. A total of 19 WRKY transcriptome sequences with complete ORFs were identified as WRKY transcriptome sequences by searching for WRKYs in RNA sequencing data. Then, the WRKY genes were classified by phylogenetic and conserved motif analysis of the WRKY family in Cynanchum thesioides and Arabidopsis thaliana. qRT-PCR was used to determine the expression patterns of 19 CtWRKY genes in different tissues and seedlings of Cynanchum thesioides under plant hormone (ABA and ETH) and abiotic stresses (cold and salt). The results showed that 19 CtWRKY genes could be divided into groups I-III according to their structure and phylogenetic characteristics, and group II could be divided into five subgroups. The prediction of CtWRKY gene protein interactions indicates that CtWRKY is involved in many biological processes. In addition, the CtWRKY gene was differentially expressed in different tissues and positively responded to abiotic stress and phytohormone treatment, among which CtWRKY9, CtWRKY18, and CtWRKY19 were significantly induced under various stresses. This study is the first to identify the WRKY gene family in Cynanchum thesioides, and the systematic analysis lays a foundation for further identification of the function of WRKY genes in Cynanchum thesioides.

Global transcriptome analysis reveals potential genes associated with genic male sterility of rapeseed (Brassica napus L.)

Rapeseed is the third leading source of edible oil in the world. Genic male sterility (GMS) lines provide crucial material for harnessing heterosis for rapeseed. GMS lines have been widely used successfully for rapeseed hybrid production. The physiological and molecular mechanism of pollen development in GMS lines of rapeseed (Brassica napus L.) need to be determined for the creation of hybrids and cultivation of new varieties. However, limited studies have focused on systematically mining genes that regulate the pollen development of GMS lines in B. napus. In the present study, to determine the stage at which pollen development begins to show abnormality in the GMS lines, we performed semi-thin section analysis of the anthers with five pollen development stages. The results indicated that the abnormal pollen development in DGMS lines might start at the meiotic stage, and abnormal pollen development in RGMS lines probably occurred before the tetrad stage. To investigate the critical genes and pathways involved in pollen development in GMS lines, we constructed and sequenced 24 transcriptome libraries for the flower buds from the fertile and sterile lines of two recessive GMS (RGMS) lines (6251AB and 6284AB) and two dominant GMS (DGMS) lines (4001AB and 4006AB). A total of 23,554 redundant DEGs with over two-fold change between sterile and fertile lines were obtained. A total of 346 DEGs were specifically related to DGMS, while 1,553 DEGs were specifically related to RGMS. A total of 1,545 DEGs were shared between DGMS and RGMS. And 253 transcription factors were found to be differentially expressed between the sterile and fertile lines of GMS. In addition, 6,099 DEGs possibly related to anther, pollen, and microspore development processes were identified. Many of these genes have been reported to be involved in anther and microspore developmental processes. Several DEGs were speculated to be key genes involved in the regulation of fertility. Three differentially expressed genes were randomly selected and their expression levels were verified by quantitative PCR (qRT-PCR). The results of qRT-PCR largely agreed with the transcriptome sequencing results. Our findings provide a global view of genes that are potentially involved in GMS occurrence. The expression profiles and function analysis of these DEGs were provided to expand our understanding of the complex molecular mechanism in pollen and sterility development in B. napus.

Cellulose accumulation mediated by PbrCSLD5, a cellulose synthase-like protein, results in cessation of pollen tube growth in Pyrus bretschneideri

Cellulose, a key component of the cell wall, plays an important role in maintaining the growth of pollen tubes. However, the molecular mechanism of cellulose participating in the cessation of pear pollen tube growth remains unclear. Here, we reported that at 15 h post-cultured (HPC), the slow-growth pear pollen tubes showed thickened cell walls and cellulose accumulation in the inner wall. Transcriptome data and quantitative real-time PCR analysis showed that PbrCSLD5, a cellulose synthesis-like gene, was highly expressed in the 15 HPC pear pollen tubes. Knockdown of PbrCSLD5 caused a decrease in cellulose content in pear pollen tubes. Moreover, PbrCSLD5 overexpression in Arabidopsis resulted in the accumulation of cellulose and disruption of normal pollen tube growth. Transcription factor PbrMADS52 was found to bind to the promoter of PbrCSLD5 and enhanced its expression. Our results suggested that the PbrMADS52-PbrCSLD5 signaling pathway led to increased cellulose content in the pear pollen tube cell wall, thereby inhibiting pollen tube growth. These results provided new insights into the regulation of pollen tube growth.

The poly(A) polymerase PAPS1 mediates pollen maturation by regulating sperm cell differentiation in plants

In flowering plants, a haploid microspore undergoes an asymmetric division to produce the male germline that encounters a mitotic division to produce two germ cells. The resulting germ cells undergo a series of specialization events to produce the two sperm cells required for double fertilization. These events include to upregulate male germline-specific while downregulating male germline-nonspecific regulon, but how these specializations events are regulated, are still unresolved. To know how plant sperm cell is specialized, we mutagenized Arabidopsis double homozygous transgenic line (MGH3p-MGH3::eGFP and ACTIN11p-H2B::mRFP) by an ethyl methane sulfonate (EMS) treatment and isolated a mutant with sperms identity loss, resulting in a completely male defective plant. Second-generation sequencing identified a point mutation G/A causing premature stop codon TGG/TGA in the poly(A) polymerase PAPS1 that is linked with phenotype. Further, we found that paps1 mutant fails to upregulate male germline-specific regulon and to downregulate male germline-nonspecific factors required for sperm cell differentiation and attaining pollen maturation. Previously, polyadenylation of pre-mRNAs by PAPS1 has been found crucial for both RNA-based silencing processes and the processing of pre-mRNAs into mature mRNAs ready for translation. This study concludes that PAPS1 mediates sperm cell differentiation through upregulating specific while silencing the nonspecific factors of male germlines.

Auxin boosts energy generation pathways to fuel pollen maturation in barley

Pollen grains become increasingly independent of the mother plant as they reach maturity through poorly understood developmental programs. We report that the hormone auxin is essential during barley pollen maturation to boost the expression of genes encoding almost every step of heterotrophic energy production pathways. Accordingly, auxin is necessary for the flux of sucrose and hexoses into glycolysis and to increase the levels of pyruvate and two tricarboxylic (TCA) cycle metabolites (citrate and succinate). Moreover, bioactive auxin is synthesized by the pollen-localized enzyme HvYUCCA4, supporting that pollen grains autonomously produce auxin to stimulate a specific cellular output, energy generation, that fuels maturation processes such as starch accumulation. Our results demonstrate that auxin can shift central carbon metabolism to drive plant cell development, which suggests a direct mechanism for auxin's ability to promote growth and differentiation.

Transcriptome Analysis of Triple Mutant for OsMADS62, OsMADS63, and OsMADS68 Reveals the Downstream Regulatory Mechanism for Pollen Germination in Rice (Oryza sativa)

The MADS (MCM1-AGAMOUS-DEFFICIENS-SRF) gene family has a preserved domain called MADS-box that regulates downstream gene expression as a transcriptional factor. Reports have revealed three MADS genes in rice, OsMADS62, OsMADS63, and OsMADS68, which exhibits preferential expression in mature rice pollen grains. To better understand the transcriptional regulation of pollen germination and tube growth in rice, we generated the loss-of-function homozygous mutant of these three OsMADS genes using the CRISPR-Cas9 (clustered regularly interspaced short palindromic repeats-CRISPR associated protein 9) system in wild-type backgrounds. Results showed that the triple knockout (KO) mutant showed a complete sterile phenotype without pollen germination. Next, to determine downstream candidate genes that are transcriptionally regulated by the three OsMADS genes during pollen development, we proceeded with RNA-seq analysis by sampling the mature anther of the mutant and wild-type. Two hundred and seventy-four upregulated and 658 downregulated genes with preferential expressions in the anthers were selected. Furthermore, downregulated genes possessed cell wall modification, clathrin coat assembly, and cellular cell wall organization features. We also selected downregulated genes predicted to be directly regulated by three OsMADS genes through the analyses for promoter sequences. Thus, this study provides a molecular background for understanding pollen germination and tube growth mediated by OsMADS62, OsMADS63, and OsMADS68 with mature pollen preferred expression.

Proteomic Analysis of Generative and Vegetative Nuclei Reveals Molecular Characteristics of Pollen Cell Differentiation in Lily

In plants, the cell fates of a vegetative cell (VC) and generative cell (GC) are determined after the asymmetric division of the haploid microspore. The VC exits the cell cycle and grows a pollen tube, while the GC undergoes further mitosis to produce two sperm cells for double fertilization. However, our understanding of the mechanisms underlying their fate differentiation remains limited. One major advantage of the nuclear proteome analysis is that it is the only method currently able to uncover the systemic differences between VC and GC due to GC being engulfed within the cytoplasm of VC, limiting the use of transcriptome. Here, we obtained pure preparations of the vegetative cell nuclei (VNs) and generative cell nuclei (GNs) from germinating lily pollens. Utilizing these high-purity VNs and GNs, we compared the differential nucleoproteins between them using state-of-the-art quantitative proteomic techniques. We identified 720 different amount proteins (DAPs) and grouped the results in 11 fate differentiation categories. Among them, we identified 29 transcription factors (TFs) and 10 cell fate determinants. Significant differences were found in the molecular activities of vegetative and reproductive nuclei. The TFs in VN mainly participate in pollen tube development. In comparison, the TFs in GN are mainly involved in cell differentiation and male gametogenesis. The identified novel TFs may play an important role in cell fate differentiation. Our data also indicate differences in nuclear pore complexes and epigenetic modifications: more nucleoporins synthesized in VN; more histone variants and chaperones; and structural maintenance of chromosome (SMC) proteins, chromatin remodelers, and DNA methylation-related proteins expressed in GN. The VC has active macromolecular metabolism and mRNA processing, while GC has active nucleic acid metabolism and translation. Moreover, the members of unfolded protein response (UPR) and programmed cell death accumulate in VN, and DNA damage repair is active in GN. Differences in the stress response of DAPs in VN vs. GN were also found. This study provides a further understanding of pollen cell differentiation mechanisms and also a sound basis for future studies of the molecular mechanisms behind cell fate differentiation.

Nepenthes × ventrata Transcriptome Profiling Reveals a Similarity Between the Evolutionary Origins of Carnivorous Traps and Floral Organs

The emergence of the carnivory syndrome and traps in plants is one of the most intriguing questions in evolutionary biology. In the present study, we addressed it by comparative transcriptomics analysis of leaves and leaf-derived pitcher traps from a predatory plant Nepenthes ventricosa × Nepenthes alata. Pitchers were collected at three stages of development and a total of 12 transcriptomes were sequenced and assembled de novo. In comparison with leaves, pitchers at all developmental stages were found to be highly enriched with upregulated genes involved in stress response, specification of shoot apical meristem, biosynthesis of sucrose, wax/cutin, anthocyanins, and alkaloids, genes encoding digestive enzymes (proteases and oligosaccharide hydrolases), and flowering-related MADS-box genes. At the same time, photosynthesis-related genes in pitchers were transcriptionally downregulated. As the MADS-box genes are thought to be associated with the origin of flower organs from leaves, we suggest that Nepenthes species could have employed a similar pathway involving highly conserved MADS-domain transcription factors to develop a novel structure, pitcher-like trap, for capture and digestion of animal prey during the evolutionary transition to carnivory. The data obtained should clarify the molecular mechanisms of trap initiation and development and may contribute to solving the problem of its emergence in plants.

The MADS-Box Transcription Factor EjAGL65 Controls Loquat Flesh Lignification via Direct Transcriptional Inhibition of EjMYB8

Loquat fruit accumulates lignin in its flesh when undergoing chilling injury during postharvest storage, making it a suitable model for the study of flesh lignification. Transcriptional regulation of lignin biosynthesis is principally controlled by the NAC-MYB transcriptional cascade in model plants. Previous research has demonstrated that EjMYB8 activates lignin biosynthesis through direct interaction with the promoter of Ej4CL1. However, the classic NAC-MYB gene regulation network has not been established. Here, the MADS-box gene EjAGL65 was discovered by screening a cDNA library using the EjMYB8 promoter as bait in yeast. A phylogenetic analysis and structural comparisons revealed that EjAGL65 belongs to the Mδ subgroup of the MADS-box family, whose members have not been reported as being involved in the regulation of lignin deposition. EjAGL65 transcription was downregulated at 0°C compared to 5°C, indicating a negative correlation with the change of lignin content. A dual-luciferase assay indicated that EjAGL65 is capable of inhibiting the promoter activity of EjMYB8 in vivo. These results showed that the Mδ MADS-box gene EjAGL65 transcriptionally regulates EjMYB8 during postharvest chilling induced flesh lignification, which differs from the classical regulation model of lignin biosynthesis that has been illustrated for developmental lignin accumulation.

Characteristics of banana B genome MADS-box family demonstrate their roles in fruit development, ripening, and stress

MADS-box genes are critical regulators of growth and development in flowering plants. Sequencing of the Musa balbisiana (B) genome has provided a platform for the systematic analysis of the MADS-box gene family in the important banana ancestor Musa balbisiana. Seventy-seven MADS-box genes, including 18 type I and 59 type II, were strictly identified from the banana (Pisang Klutuk Wulung, PKW, 2n = 2x = 22) B genome. These genes have been preferentially placed on the banana B genome. Evolutionary analysis suggested that M. balbisiana MCM1-AGAMOUS-DEFICIENS-SRF (MbMADS) might be organized into the MIKC^c, MIKC*, Mα, Mβ, and Mγ groups according to the phylogeny. MIKC^c was then further categorized into 10 subfamilies according to conserved motif and gene structure analyses. The well-defined MADS-box genes highlight gene birth and death in banana. MbMADSes originated from the same ancestor as MaMADSes. Transcriptome analysis in cultivated banana (ABB) revealed that MbMADSes were conserved and differentially expressed in several organs, in various fruit developing and ripening stages, and in stress treatments, indicating the participation of these genes in fruit development, ripening, and stress responses. Of note, SEP/AGL2 and AG, as well as other several type II MADS-box genes, including the STMADS11 and TM3/SOC1 subfamilies, indicated elevated expression throughout banana fruit development, ripening, and stress treatments, indicating their new parts in controlling fruit development and ripening. According to the co-expression network analysis, MbMADS75 interacted with bZIP and seven other transcription factors to perform its function. This systematic analysis reveals fruit development, ripening, and stress candidate MbMADSes genes for additional functional studies in plants, improving our understanding of the transcriptional regulation of MbMADSes genes and providing a base for genetic modification of MADS-mediated fruit development, ripening, and stress.

Depletion of KNL2 Results in Altered Expression of Genes Involved in Regulation of the Cell Cycle, Transcription, and Development in Arabidopsis

Centromeres contain specialized nucleosomes at which histone H3 is partially replaced by the centromeric histone H3 variant cenH3 that is required for the assembly, maintenance, and proper function of kinetochores during mitotic and meiotic divisions. Previously, we identified a KINETOCHORE NULL 2 (KNL2) of Arabidopsis thaliana that is involved in the licensing of centromeres for the cenH3 recruitment. We also demonstrated that a knockout mutant for KNL2 shows mitotic and meiotic defects, slower development, reduced growth rate, and fertility. To analyze an effect of KNL2 mutation on global gene transcription of Arabidopsis, we performed RNA-sequencing experiments using seedling and flower bud tissues of knl2 and wild-type plants. The transcriptome data analysis revealed a high number of differentially expressed genes (DEGs) in knl2 plants. The set was enriched in genes involved in the regulation of the cell cycle, transcription, development, and DNA damage repair. In addition to comprehensive information regarding the effects of KNL2 mutation on the global gene expression, physiological changes in plants are also presented, which provides an integrated understanding of the critical role played by KNL2 in plant growth and development.

Investigation of the genes associated with a male sterility mutant (msm) in Chinese cabbage (Brassica campestris ssp. pekinensis) using RNA-Seq

In Chinese cabbage, hybrid seed production is performed using male sterility lines, an important approach to heterosis utilization. In this study, a stably inherited male sterile mutant msm was obtained from the 'FT'-doubled haploid line of Chinese cabbage using isolated microspore culture combined with ⁶⁰Co γ-ray mutagenesis. The genetic backgrounds of 'FT' and msm were highly consistent; however, compared with wild-type 'FT', msm exhibited completely degenerated stamens and no pollen phenotype. Other characters showed no significant differences. Cytological observations revealed that stamen abortion in msm begins during the tetrad period and that tapetum cells were abnormally expanded and highly vacuolated, leading to microspore abortion. Genetic analysis indicated that the msm mutant phenotype is controlled by a single recessive nuclear gene. Comparative transcriptome analysis of 'FT' and msm flower buds using RNA-Seq technology revealed 1653 differentially expressed genes, among which, a large number associated with male sterility were detected, including 64 pollen development- and pollen tube growth-related genes, 94 pollen wall development-related genes, 11 phytohormone-related genes, and 16 transcription factor-related genes. An overwhelming majority of these genes were down-regulated in msm compared with 'FT'. Furthermore, KEGG pathway analysis indicated that a variety of carbohydrate metabolic and lipid metabolic pathways were significantly enriched, which may be related to pollen abortion. The expression patterns of 24 male sterility-related genes were analyzed using qRT-PCR. In addition, 24,476 single-nucleotide polymorphisms and 413,073 insertion-deletion events were specifically detected in msm. These results will facilitate elucidation of the regulatory mechanisms underlying male sterility in Chinese cabbage.

OsCPK21 is required for pollen late-stage development in rice

In flowering plants, pollen development is a critical step for reproductive success and necessarily involves complex genetic regulatory networks. Calcium-dependent protein kinases (CPKs) are plant-specific calcium sensors involved in the regulation of plant development and adaption to the environment; however, whether they play a role in regulating male reproduction remains elusive. Here, we found that the knockdown of spikelet-specific OsCPK21 causes pollen abortion in OsCPK21-RNAi transgenic plants. Severe defects in pollen development initiated at stage 10 of anther development and simultaneous cell death occurred in the pollen cells of OsCPK21-RNAi plants. Microarray analysis and qRT-PCR revealed that the transcription of OsCPK21 is coordinated with that of MIKC*-type MADS box transcription factors OsMADS62, OsMADS63, and OsMADS68 during rice anther development. We further showed that OsCPK21 indirectly up-regulates the transcription of OsMADS62, OsMADS63, and OsMADS68 through the potential MYB binding site, DRE/CRT element, and/or new ERF binding motif localised in the promoter region of these three MADS genes. These findings suggest that OsCPK21 plays an essential role in pollengenesis, possibly via indirectly regulating the transcription of MIKC*-type MADS box proteins.

MADS-box genes underground becoming mainstream: plant root developmental mechanisms

Plant growth is largely post-embryonic and depends on meristems that are active throughout the lifespan of an individual. Developmental patterns rely on the coordinated spatio-temporal expression of different genes, and the activity of transcription factors is particularly important during most morphogenetic processes. MADS-box genes constitute a transcription factor family in eukaryotes. In Arabidopsis, their proteins participate in all major aspects of shoot development, but their role in root development is still not well characterized. In this review we synthetize current knowledge pertaining to the function of MADS-box genes highly expressed in roots: XAL1, XAL2, ANR1 and AGL21, as well as available data for other MADS-box genes expressed in this organ. The role of Trithorax group and Polycomb group complexes on MADS-box genes' epigenetic regulation is also discussed. We argue that understanding the role of MADS-box genes in root development of species with contrasting architectures is still a challenge. Finally, we propose that MADS-box genes are key components of the gene regulatory networks that underlie various gene expression patterns, each one associated with the distinct developmental fates observed in the root. In the case of XAL1 and XAL2, their role within these networks could be mediated by regulatory feedbacks with auxin.

miR167 limits anther growth to potentiate anther dehiscence

In flowering plants, anther dehiscence and pollen release are essential for sexual reproduction. Anthers dehisce after cell wall degradation weakens stomium cell junctions in each anther locule, and desiccation creates mechanical forces that open the locules. Either effect or both together may break stomium cell junctions. The microRNA miR167 negatively regulates ARF6 and ARF8, which encode auxin response transcription factors. Arabidopsis mARF6 or mARF8 plants with mutated miR167 target sites have defective anther dehiscence and ovule development. Null mir167a mutations recapitulated mARF6 and mARF8 anther and ovule phenotypes, indicating that MIR167a is the main miR167 precursor gene that delimits ARF6 and ARF8 expression in these organs. Anthers of mir167a or mARF6/8 plants overexpressed genes encoding cell wall loosening functions associated with cell expansion, and grew larger than wild-type anthers did starting at flower stage 11. Experimental desiccation enabled dehiscence of miR167-deficient anthers, indicating competence to dehisce. Conversely, high humidity conditions delayed anther dehiscence in wild-type flowers. These results support a model in which miR167-mediated anther growth arrest permits anther dehiscence. Without miR167 regulation, excess anther growth delays dehiscence by prolonging desiccation.

Asymmetric birth and death of type I and type II MADS-box gene subfamilies in the rubber tree facilitating laticifer development

The rubber tree (Hevea brasiliensis Muell. Arg.) is a rubber producing crop and contains specialized laticifers. MADS-box genes are a family of transcription factor genes that regulate plant development, especially floral organ and gametophyte development. 97 MADS-box genes were identified in the rubber tree through transcriptomes and genome mining. 93.8% of the genes were mapped onto the genome scaffolds in correspondence to the coverage (93.8%) of current version of sequenced genome. Phylogenetic analysis indicates that type II MADS-box genes have been more actively duplicated than their orthologous genes in Arabidopsis and rice, so that most (70, 72.2%) of the MADS-box genes in the rubber tree belong to type II subfamily. This is a high percentage compared to those in Arabidopsis (43.7%) and rice (56.8%). Moreover, 69 out of 70 type II genes in the rubber tree are transcribed, and they are mostly predominantly expressed in flowers, but some genes are predominantly expressed in laticifers, suggesting their roles in both flower and laticifer development. The number of type I genes in the rubber tree is only 27 (27.8%), a much smaller number compared to their orthologous genes in Arabidopsis (56.3%) and rice (43.2%). At the same time, most of the type I genes (55.6%, 15) in the rubber tree are silent and are probably pseudogenes. The high birth rate and low death rate of type II genes and low birth rate and high death rate of type I genes may corresponds to special developmental requirements in the rubber tree, e.g. the development of laticifer system for biosynthesis of cis-polyisoprene, the rubber. Moreover, atypical MIKC* factors (e.g. HbMADS1 in S-clade, and HbMADS20 in P-clade) are identified. These genes are diverged to typical MIKC* genes in sequences and facilitate functions required in laticifer development and rubber biosynthesis, which is not necessary in Arabidopsis and rice.

Genome sequence of Jatropha curcas L., a non-edible biodiesel plant, provides a resource to improve seed-related traits

Jatropha curcas (physic nut), a non-edible oilseed crop, represents one of the most promising alternative energy sources due to its high seed oil content, rapid growth and adaptability to various environments. We report ~339 Mbp draft whole genome sequence of J. curcas var. Chai Nat using both the PacBio and Illumina sequencing platforms. We identified and categorized differentially expressed genes related to biosynthesis of lipid and toxic compound among four stages of seed development. Triacylglycerol (TAG), the major component of seed storage oil, is mainly synthesized by phospholipid:diacylglycerol acyltransferase in Jatropha, and continuous high expression of homologs of oleosin over seed development contributes to accumulation of high level of oil in kernels by preventing the breakdown of TAG. A physical cluster of genes for diterpenoid biosynthetic enzymes, including casbene synthases highly responsible for a toxic compound, phorbol ester, in seed cake, was syntenically highly conserved between Jatropha and castor bean. Transcriptomic analysis of female and male flowers revealed the up-regulation of a dozen family of TFs in female flower. Additionally, we constructed a robust species tree enabling estimation of divergence times among nine Jatropha species and five commercial crops in Malpighiales order. Our results will help researchers and breeders increase energy efficiency of this important oil seed crop by improving yield and oil content, and eliminating toxic compound in seed cake for animal feed.

MADS-Box Genes Are Key Components of Genetic Regulatory Networks Involved in Abiotic Stress and Plastic Developmental Responses in Plants

Plants, as sessile organisms, adapt to different stressful conditions, such as drought, salinity, extreme temperatures, and nutrient deficiency, via plastic developmental and growth responses. Depending on the intensity and the developmental phase in which it is imposed, a stress condition may lead to a broad range of responses at the morphological, physiological, biochemical, and molecular levels. Transcription factors are key components of regulatory networks that integrate environmental cues and concert responses at the cellular level, including those that imply a stressful condition. Despite the fact that several studies have started to identify various members of the MADS-box gene family as important molecular components involved in different types of stress responses, we still lack an integrated view of their role in these processes. In this review, we analyze the function and regulation of MADS-box gene family members in response to drought, salt, cold, heat, and oxidative stress conditions in different developmental processes of several plants. In addition, we suggest that MADS-box genes are key components of gene regulatory networks involved in plant responses to stress and plant developmental plasticity in response to seasonal changes in environmental conditions.

The evolution and patterning of male gametophyte development

The reproductive adaptations of land plants have played a key role in their terrestrial colonization and radiation. This encompasses mechanisms used for the production, dispersal and union of gametes to support sexual reproduction. The production of small motile male gametes and larger immotile female gametes (oogamy) in specialized multicellular gametangia evolved in the charophyte algae, the closest extant relatives of land plants. Reliance on water and motile male gametes for sexual reproduction was retained by bryophytes and basal vascular plants, but was overcome in seed plants by the dispersal of pollen and the guided delivery of non-motile sperm to the female gametes. Here we discuss the evolutionary history of male gametogenesis in streptophytes (green plants) and the underlying developmental biology, including recent advances in bryophyte and angiosperm models. We conclude with a perspective on research trends that promise to deliver a deeper understanding of the evolutionary and developmental mechanisms of male gametogenesis in plants.

Promoters of Arabidopsis Casein kinase I-like 2 and 7 confer specific high-temperature response in anther

KEY MESSAGE: (1) We systematically analyze the promoter activities of AtCKLs in various tissues; (2) AtCKL2 and AtCKL7 were expressed in early developmental anthers under high temperature (HT) conditions; (3) AtMYB24 may function as a positive regulator of AtCKL2 and AtCKL7 expression under HT. High temperature (HT) can seriously impede plant growth and development, causing severe loss of crop production. In Arabidopsis, AtCKL genes show high similarity to GhCKI, a gene reported to disrupt tapetal programmed cell death in cotton. However, most of AtCKL genes are not well characterized. Here, we systematically analyzed the expression patterns of AtCKLs in various tissues. The expression of AtCKL2 and AtCKL7 was induced in early anther development under HT, which is similar to the case of GhCKI. In silico analysis of AtCKL2 and AtCKL7 promoters indicated that four types of transcription factors (TFs) (MADS, NAC, WRKY and R2R3-MYB) might bind to AtCKL2 and AtCKL7 promoters. Furthermore, three MADS, three NAC, one WRKY, and three R2R3-MYB TFs were up-regulated in stage 1-8 anthers and three R2R3-MYB TFs were up-regulated in stage 9-10 anthers under HT, implying the important roles of R2R3-MYB genes in the response of anthers to HT. Among the R2R3-MYB genes, AtMYB24 showed the similar expression as AtCKL2 and AtCKL7 in the anthers under HT. Additionally, yeast one-hybrid and dual-luciferase reporter system assays verified that AtMYB24 could bind to AtCKL2 and AtCKL7 promoters and activate the expression of these two genes. In brief, this study provides the overall expression profiles of AtCKLs, useful information for unraveling the molecular mechanism of AtCKL2 and AtCKL7 gene expression in early anther development under HT, and important clues for elucidating the mechanism of transcriptional regulation of CKI genes in plant anther under HT, which are critical to the reduction of crop yield loss resulting from HT.

WRKY2/34-VQ20 Modules in Arabidopsis thaliana Negatively Regulate Expression of a Trio of Related MYB Transcription Factors During Pollen Development

Male gametogenesis in plants is tightly controlled and involves the complex and precise regulation of transcriptional reprogramming. Interactions between WRKY proteins and VQ motif-containing proteins are required to control these complicated transcriptional networks. However, our understanding of the mechanisms by which these complexes affect downstream gene expression is quite limited. In this study, we found that WRKY2 and WKRY34 repress MYB97, MYB101, and MYB120 expression during male gametogenesis. MYB expression was up-regulated in the wrky2-1 wrky34-1 vq20-1 triple mutant during male gametogenesis. The expression levels of six potential targets of the three MYBs increased the most in the wrky2-1 wrky34-1 vq20-1 triple mutant, followed by the wrky2-1 wrky34-1 double mutant, compared with in wild-type. Yeast one-hybrid and dual luciferase reporter assays indicated that WRKY2 and WRKY34 recognized the MYB97 promoter by binding to its W-boxes. MYB97 overexpression caused defects in pollen germination and pollen tube length, which impacted male fertility. Thus, WRKY2/34-VQ20 complexes appear to negatively regulate the expression of certain MYBs during plant male gametogenesis.

Molecular cloning, characterization and expression analysis of MADS-box genes associated with reproductive development in Momordica dioica Roxb

The repertoire and functions of MADS-box family transcription factors (TFs) largely remains unexplored with respect to floral organogenesis of Momordica dioica Roxb. Degenerative PCR followed by rapid amplification of cDNA ends was employed in the present study to clone and characterize 17 MADS-box genes (designated as MdMADS01 to MdMADS17) from the floral buds of M. dioica. The cloned genes were clustered into three subgroups (11 MIKC^C, 4 MIKC* and 2 Mα) based on phylogenetic relationships with the MADS-box genes from Cucumis sativus, Cucumis melo and Arabidopsis thaliana. Southern hybridization showed that all the isolated genes were represented by single copy locus in the M. dioica genome. Gene structure analysis revealed 1-8 exons in MdMADS-box genes with the number of exons in MIKC greatly exceeding from that in M-type genes. Motif elicitation of the MdMADS-box genes indicated the presence of additional domains with MIKC type, suggesting that they had more complex structures. Expression analysis of MdMADS genes in six M. dioica transcriptome suggested that, 11 MIKC^C-type genes are associated with floral homeotic functions, 4 MIKC*-type genes (MdMADS12 to MdMADS15) controlled the growth of male gametophyte, while the two M-type genes (MdMADS16 and MdMADS17) played significant role in female gametogenesis and seed development. Overall, these are the first set of MADS-box genes from M. dioica exhibiting a differential expression pattern during floral development. The results from this study will provide valuable information for further functional studies of candidate MADS-box genes in the sexual dimorphism of this economically important dioecious cucurbit.

Phylogenomic Synteny Network Analysis of MADS-Box Transcription Factor Genes Reveals Lineage-Specific Transpositions, Ancient Tandem Duplications, and Deep Positional Conservation

Conserved genomic context provides critical information for comparative evolutionary analysis. With the increase in numbers of sequenced plant genomes, synteny analysis can provide new insights into gene family evolution. Here, we exploit a network analysis approach to organize and interpret massive pairwise syntenic relationships. Specifically, we analyzed synteny networks of the MADS-box transcription factor gene family using 51 completed plant genomes. In combination with phylogenetic profiling, several novel evolutionary patterns were inferred and visualized from synteny network clusters. We found lineage-specific clusters that derive from transposition events for the regulators of floral development (APETALA3 and PI) and flowering time (FLC) in the Brassicales and for the regulators of root development (AGL17) in Poales. We also identified two large gene clusters that jointly encompass many key phenotypic regulatory Type II MADS-box gene clades (SEP1, SQUA, TM8, SEP3, FLC, AGL6, and TM3). Gene clustering and gene trees support the idea that these genes are derived from an ancient tandem gene duplication that likely predates the radiation of the seed plants and then expanded by subsequent polyploidy events. We also identified angiosperm-wide conservation of synteny of several other less studied clades. Combined, these findings provide new hypotheses for the genomic origins, biological conservation, and divergence of MADS-box gene family members.

Overexpression of a tomato miR171 target gene SlGRAS24 impacts multiple agronomical traits via regulating gibberellin and auxin homeostasis

In Arabidopsis, the miR171-GRAS module has been clarified as key player in meristem maintenance. However, the knowledge about its role in fruit crops like tomato (Solanum lycopersicum) remains scarce. We previously identified tomato SlGRAS24 as a target gene of Sly-miR171. To study the role of this probable transcription factor, we generated transgenic tomato plants underexpressing SlGRAS24, overexpressing SlGRAS24, overexpressing Sly-miR171 and expressing β-glucuronidase (GUS) under the SlGRAS24 promoter (proSlGRAS24-GUS). Plants overexpressing SlGRAS24 (SlGRAS24-OE) had pleiotropic phenotypes associated with multiple agronomical traits including plant height, flowering time, leaf architecture, lateral branch number, root length, fruit set and development. Many GA/auxin-related genes were down-regulated and altered responsiveness to exogenous IAA/NAA or GA₃ application was observed in SlGRAS24-OE seedlings. Moreover, compromised fruit set and development in SlGRAS24-OE was also observed. These newly identified phenotypes for SlGRAS24 homologs in tomato were later proved to be caused by impaired pollen sacs and fewer viable pollen grains. At anthesis, the comparative transcriptome results showed altered expression of genes involved in pollen development and hormone signalling. Taken together, our data demonstrate that SlGRAS24 participates in a series of developmental processes through modulating gibberellin and auxin signalling, which sheds new light on the involvement of hormone crosstalk in tomato development.

Characterization of pollen-expressed bZIP protein interactions and the role of ATbZIP18 in the male gametophyte

KEY MESSAGE : bZIP TF network in pollen. Transcriptional control of gene expression represents an important mechanism guiding organisms through developmental processes and providing plasticity towards environmental stimuli. Because of their sessile nature, plants require effective gene regulation for rapid response to variation in environmental and developmental conditions. Transcription factors (TFs) provide such control ensuring correct gene expression in spatial and temporal manner. Our work reports the interaction network of six bZIP TFs expressed in Arabidopsis thaliana pollen and highlights the potential functional role for AtbZIP18 in pollen. AtbZIP18 was shown to interact with three other pollen-expressed bZIP TFs-AtbZIP34, AtbZIP52, and AtbZIP61 in yeast two-hybrid assays. AtbZIP18 transcripts are highly expressed in pollen, and at the subcellular level, an AtbZIP18-GFP fusion protein was located in the nucleus and cytoplasm/ER. To address the role of AtbZIP18 in the male gametophyte, we performed phenotypic analysis of a T-DNA knockout allele, which showed slightly reduced transmission through the male gametophyte. Some of the phenotype defects in atbzip18 pollen, although observed at low penetrance, were similar to those seen at higher frequency in the T-DNA knockout of the interacting partner, AtbZIP34. To gain deeper insight into the regulatory role of AtbZIP18, we analysed atbzip18/- pollen microarray data. Our results point towards a potential repressive role for AtbZIP18 and its functional redundancy with AtbZIP34 in pollen.

Transcriptome Analysis of Hamelia patens (Rubiaceae) Anthers Reveals Candidate Genes for Tapetum and Pollen Wall Development

Studies of the anther transcriptome on non-model plants without a known genome are surprisingly scarce. RNA-Seq and digital gene expression (DGE) profiling provides a comprehensive approach to identify candidate genes contributing to developmental processes in non-model species. Here we built a transcriptome library of developing anthers of Hamelia patens and analyzed DGE profiles from each stage to identify genes that regulate tapetum and pollen development. In total 7,720 putative differentially expressed genes across four anther stages were identified. The number of putative stage-specific genes was: 776 at microspore mother cell stage, 807 at tetrad stage, 322 at uninucleate microspore stage, and the highest number (1,864) at bicellular pollen stage. GO enrichment analysis revealed 243 differentially expressed and 108 stage-specific genes that are potentially related to tapetum development, sporopollenin synthesis, and pollen wall. The number of expressed genes, their function and expression profiles were all significantly correlated with anther developmental processes. Overall comparisons of anther and pollen transcriptomes with those of rice and Arabidopsis together with the expression profiles of homologs of known anther-expressed genes, revealed conserved patterns and also divergence. The divergence may reflect taxon-specific differences in gene expression, the use RNA-seq as a more sensitive methodology, variation in tissue composition and sampling strategies. Given the lack of genomic sequence, this study succeeded in assigning putative identity to a significant proportion of anther-expressed genes and genes relevant to tapetum and pollen development in H. patens. The anther transcriptome revealed a molecular distinction between developmental stages, serving as a resource to unravel the functions of genes involved in anther development in H. patens and informing the analysis of other members of the Rubiaceae.

A Conserved cis-Regulatory Module Determines Germline Fate through Activation of the Transcription Factor DUO1 Promoter

The development of the male germline within pollen relies upon the activation of numerous target genes by the transcription factor DUO POLLEN1 (DUO1). The expression of DUO1 is restricted to the male germline and is first detected shortly after the asymmetric division that segregates the germ cell lineage. Transcriptional regulation is critical in controlling DUO1 expression, since transcriptional and translational fusions show similar expression patterns. Here, we identify key promoter sequences required for the germline-specific regulation of DUO1 transcription. Combining promoter deletion analyses with phylogenetic footprinting in eudicots and in Arabidopsis accessions, we identify a cis-regulatory module, Regulatory region of DUO1 (ROD1), which replicates the expression pattern of DUO1 in Arabidopsis (Arabidopsis thaliana). We show that ROD1 from the legume Medicago truncatula directs male germline-specific expression in Arabidopsis, demonstrating conservation of DUO1 regulation among eudicots. ROD1 contains several short conserved cis-regulatory elements, including three copies of the motif DNGTGGV, required for germline expression and tandem repeats of the motif YAACYGY, which enhance DUO1 transcription in a positive feedback loop. We conclude that a cis-regulatory module conserved in eudicots directs the spatial and temporal expression of the transcription factor DUO1 to specify male germline fate and sperm cell differentiation.

Mapping and Characterization of the fefe Gene That Controls Iron Uptake in Melon (Cucumis melo L.)

Iron (Fe) deficiency in plants limits crop growth and productivity. Molecular mechanisms that plants use to sense and respond to Fe deficiency by coordinated expression of Fe-uptake genes are not fully understood. The C940-fe chlorotic melon (Cucumis melo) mutant known as fefe is unable to upregulate Fe-uptake genes, however, the FeFe gene had not been identified. In this study, we used two F₂ mapping populations to map and identify the FeFe gene as bHLH38, a homolog of subgroup Ib bHLH genes from Arabidopsis thaliana that are involved in transcriptional regulation of Fe-uptake genes in partnership with the FIT gene. A Ty1-copia type retrotransposon insertion of 5.056 kb within bHLH38 is responsible for the defect in bHLH38 in fefe, based on sequencing and expression analysis. This retrotransposon insertion results in multiple non-functional transcripts expected to result in an altered and truncated protein sequence. Hairy root transformation of fefe plants using wild-type bHLH38 resulted in functional complementation of the chlorotic fefe phenotype. Using a yeast-2-hybrid assay, the transcription factor Fit interacted with the wild-type bHLH38 protein, but did not interact with the fefe bHLH38 protein, suggesting that heterodimer formation of Fit/bHLH38 to regulate Fe-uptake genes does not occur in fefe roots. The second subgroup Ib bHLH gene in the melon genome is not functionally redundant to bHLH38, in contrast to Arabidopsis where four subgroup Ib bHLH genes are functionally redundant. Whereas the Arabidopsis bHLH transcript levels are upregulated by Fe deficiency, melon bHLH38 was not regulated at the transcript level. Thus, the fefe mutant may provide a platform for studying bHLH38 genes and proteins from other plant species.

Expression and Protein Interaction Analyses Reveal Combinatorial Interactions of LBD Transcription Factors During Arabidopsis Pollen Development

LATERAL ORGAN BOUNDARIES DOMAIN (LBD) transcription factor gene family members play key roles in diverse aspects of plant development. LBD10 and LBD27 have been shown to be essential for pollen development in Arabidopsis thaliana. From the previous RNA sequencing (RNA-Seq) data set of Arabidopsis pollen, we identified the mRNAs of LBD22, LBD25 and LBD36 in addition to LBD10 and LBD27 in Arabidopsis pollen. Here we conducted expression and cellular analysis using GFP:GUS (green fluorescent protein:β-glucuronidase) reporter gene and subcellular localization assays using LBD:GFP fusion proteins expressed under the control of their own promoters in Arabidopsis. We found that these LBD proteins display spatially and temporally distinct and overlapping expression patterns during pollen development. Bimolecular fluorescence complementation and GST (glutathione S-transferase) pull-down assays demonstrated that protein-protein interactions occur among the LBDs exhibiting overlapping expression during pollen development. We further showed that LBD10, LBD22, LBD25, LBD27 and LBD36 interact with each other to form heterodimers, which are localized to the nucleus in Arabidopsis protoplasts. Taken together, these results suggest that combinatorial interactions among LBD proteins may be important for their function in pollen development in Arabidopsis.

Gene expression profile indicates involvement of NO in Camellia sinensis pollen tube growth at low temperature

BACKGROUND

Nitric oxide (NO) functions as a critical signaling molecule in the low-temperature stress responses in plants, including polarized pollen tube growth in Camellia sinensis. Despite this, the potential mechanisms underlying the participation of NO in pollen tube responses to low temperature remain unclear. Here, we investigate alterations to gene expression in C. sinensis pollen tubes exposed to low-temperature stress and NO using RNA-Seq technology, in order to find the potential candidate genes related to the regulation of pollen tube elongation by NO under low-temperature stress.

RESULTS

Three libraries were generated from C. sinensis cv. 'Longjingchangye' pollen tubes cultured at 25 °C (CsPT-CK) and 4 °C (CsPT-LT) or with 25 μM DEA NONOate (CsPT-NO). The number of unigenes found for the three biological replications were 39,726, 40,440 and 41,626 for CsPT-CK; 36,993, 39,070 and 39,439 for CsPT-LT; and 39,514, 38,298 and 39,061 for CsPT-NO. A total of 36,097 unique assembled and annotated sequences from C. sinensis pollen tube reads were found in a BLAST search of the following databases: NCBI non-redundant nucleotide, Swiss-prot protein, Kyoto Encyclopedia of Genes and Genomes, Cluster of Orthologous Groups of proteins, and Gene Ontology. The absolute values of log2Ratio > 1 and probability > 0.7 were used as the thresholds for significantly differential gene expression, and 766, 497 and 929 differentially expressed genes (DEGs) were found from the comparison analyses of the CK-VS-LT, CK-VS-NO and LT-VS-NO libraries, respectively. Genes related to metabolism and signaling pathways of plant hormones, transcription factors (TFs), vesicle polarized trafficking, cell wall biosynthesis, the ubiquitination machinery of the ubiquitin system and species-specific secondary metabolite pathways were mainly observed in the CK-VS-LT and CK-VS-NO libraries.

CONCLUSION

Differentially expressed unigenes related to the inhibition of C. sinensis pollen tube growth under low temperature and NO are identified in this study. The transcriptomic gene expression profiles present a valuable genomic tool to improve studying the molecular mechanisms underlying low-temperature tolerance in pollen tube.

Expression analysis of fertilization/early embryogenesis-associated genes in Phalaenopsis orchids

One of the distinct reproductive programs in orchid species is pollination-triggered ovule development and megasporogenesis. During sexual reproduction, fertilization occurs days to months after pollination. The molecular mechanisms evolved to carry out this strategic reproductive program remain unclear. In the August issue of Plant Physiology ¹ , we report comprehensive studies of comparative genome-wide gene expression in various reproductive tissues and the molecular events associated with developmental transitions unique to sexual reproduction of Phalaenopsis aphrodite. Transcriptional factors and signaling components whose expression is specifically enriched in interior ovary tissues when fertilization occurs and embryos start to develop have been identified. Here, we report verification of additional fertilization-associated genes, DOMAINS REARRANGED METHYLTRANSFERASE 1 (PaDRM1), CHROMOMETHYLTRANSFERASE 1 (PaCMT1), SU(VAR)3-9 RELATED PROTEIN 1 (PaSUVR1), INDOLE-3-ACETIC ACID inducible 30-like 1 (PaIAA30L1), and ETHYLENE INSENSITIVE 3-like 1 (PaEIN3L1), and discuss their potential roles in gametophyte development, epigenetic reprogramming, and hormone regulation during fertilization and establishment of embryo development in Phalaenopsis orchids.

TOMATO AGAMOUS1 and ARLEQUIN/TOMATO AGAMOUS-LIKE1 MADS-box genes have redundant and divergent functions required for tomato reproductive development

Within the tomato MADS-box gene family, TOMATO AGAMOUS1 (TAG1) and ARLEQUIN/TOMATO AGAMOUS LIKE1 (hereafter referred to as TAGL1) are, respectively, members of the euAG and PLE lineages of the AGAMOUS clade. They perform crucial functions specifying stamen and carpel development in the flower and controlling late fruit development. To gain insight into the roles of TAG1 and TAGL1 genes and to better understand their functional redundancy and diversification, we characterized single and double RNAi silencing lines of these genes and analyzed expression profiles of regulatory genes involved in reproductive development. Double RNAi lines did show cell abnormalities in stamens and carpels and produced extremely small fruit-like organs displaying some sepaloid features. Expression analyses indicated that TAG1 and TAGL1 act together to repress fourth whorl sepal development, most likely through the MACROCALYX gene. Results also proved that TAG1 and TAGL1 have diversified their functions in fruit development: while TAG1 controls placenta and seed formation, TAGL1 participates in cuticle development and lignin biosynthesis inhibition. It is noteworthy that both TAG1 and double RNAi plants lacked seed development due to abnormalities in pollen formation. This seedless phenotype was not associated with changes in the expression of B-class stamen identity genes Tomato MADS-box 6 and Tomato PISTILLATA observed in silencing lines, suggesting that other regulatory factors should participate in pollen formation. Taken together, results here reported support the idea that both redundant and divergent functions of TAG1 and TAGL1 genes are needed to control tomato reproductive development.

A Transcriptome Atlas of Physcomitrella patens Provides Insights into the Evolution and Development of Land Plants

Identifying the genetic mechanisms that underpin the evolution of new organ and tissue systems is an aim of evolutionary developmental biology. Comparative functional genetic studies between angiosperms and bryophytes can define those genetic changes that were responsible for developmental innovations. Here, we report the generation of a transcriptome atlas covering most phases in the life cycle of the model bryophyte Physcomitrella patens, including detailed sporophyte developmental progression. We identified a comprehensive set of sporophyte-specific transcription factors, and found that many of these genes have homologs in angiosperms that function in developmental processes such as flowering and shoot branching. Deletion of the PpTCP5 transcription factor results in development of supernumerary sporangia attached to a single seta, suggesting that it negatively regulates branching in the moss sporophyte. Given that TCP genes repress branching in angiosperms, we suggest that this activity is ancient. Finally, comparison of P. patens and Arabidopsis thaliana transcriptomes led us to the identification of a conserved core of transcription factors expressed in tip-growing cells. We identified modifications in the expression patterns of these genes that could account for developmental differences between P. patens tip-growing cells and A. thaliana pollen tubes and root hairs.

Structural and functional annotation of the MADS-box transcription factor family in grapevine

Background

MADS-box genes encode transcription factors that are involved in developmental control and signal transduction in eukaryotes. In plants, they are associated to numerous development processes most notably those related to reproductive development: flowering induction, specification of inflorescence and flower meristems, establishment of flower organ identity, as well as regulation of fruit, seed and embryo development. Genomic analyses of MADS-box genes in different plant species are providing new relevant information on the function and evolution of this transcriptional factor family. We have performed a true genome-wide analysis of the complete set of MADS-box genes in grapevine (Vitis vinifera), analyzed their expression pattern and establish their phylogenetic relationships (including MIKC* and type I MADS-box) with genes from 16 other plant species. This study was integrated to previous works on the family in grapevine.

Results

A total of 90 MADS-box genes were detected in the grapevine reference genome by completing current gene annotations with a genome-wide analysis based on sequence similarity. We performed a thorough in-depth curation of all gene models and combined the results with gene expression information including RNAseq data to clarifying the expression of newly identified genes and improve their functional characterization. Curated data were uploaded to the ORCAE database for grapevine in the frame of the grapevine genome curation effort. This approach resulted in the identification of 30 additional MADS box genes. Among them, ten new MIKC^C genes were identified, including a potential new group of short proteins similar to the SVP protein subfamily. The MIKC* subgroup contains six genes in grapevine that can be grouped in the S (4 genes) and P (2 genes) clades, showing less redundancy than that observed in Arabidopsis thaliana. Expression pattern of these genes in grapevine is compatible with a role in male gametophyte development. Most of the identified new genes belong to the type I MADS-box genes and were classified as members of the Mα and Mγ subclasses. Ours analyses indicate that only few members of type I genes in grapevine have homology in other species and that species-specific clades appeared both in the Mα and Mγ subclasses. On the other hand, as deduced from the phylogenetic analysis with other plant species, genes that can be crucial for development of central cell, endosperm and embryos seems to be conserved in plants.

Conclusions

The genome analysis of MADS-box genes in grapevine, the characterization of their pattern of expression and the phylogenetic analysis with other plant species allowed the identification of new MADS-box genes not yet described in other plant species as well as basic characterization of their possible role, particularly in the case of type I and MIKC* genes.

Electronic supplementary material

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

Anther and pollen development: A conserved developmental pathway

Pollen development is a critical step in plant development that is needed for successful breeding and seed formation. Manipulation of male fertility has proved a useful trait for hybrid breeding and increased crop yield. However, although there is a good understanding developing of the molecular mechanisms of anther and pollen anther development in model species, such as Arabidopsis and rice, little is known about the equivalent processes in important crops. Nevertheless the onset of increased genomic information and genetic tools is facilitating translation of information from the models to crops, such as barley and wheat; this will enable increased understanding and manipulation of these pathways for agricultural improvement.

A decade of pollen transcriptomics

Overview of pollen transcriptome studies. Pollen development is driven by gene expression, and knowledge of the molecular events underlying this process has undergone a quantum leap in the last decade through studies of the transcriptome. Here, we outline historical evidence for male haploid gene expression and review the wealth of pollen transcriptome data now available. Knowledge of the transcriptional capacity of pollen has progressed from genetic studies to the direct analysis of RNA and from gene-by-gene studies to analyses on a genomic scale. Microarray and/or RNA-seq data can now be accessed for all phases and cell types of developing pollen encompassing 10 different angiosperms. These growing resources have accelerated research and will undoubtedly inspire new directions and the application of system-based research into the mechanisms that govern the development, function and evolution of angiosperm pollen.

SYBR Green-activated sorting of Arabidopsis pollen nuclei based on different DNA/RNA content

Key message: Purification of pollen nuclei. Germ cell epigenetics is a critical topic in plants and animals. The male gametophyte (pollen) of flowering plants is an attractive model to study genetic and epigenetic reprogramming during sexual reproduction, being composed of only two sperm cells contained within, its companion, vegetative cell. Here, we describe a simple and efficient method to purify SYBR Green-stained sperm and vegetative cell nuclei of Arabidopsis thaliana pollen using fluorescence-activated cell sorting to analyze chromatin and RNA profiles. The method obviates generating transgenic lines expressing cell-type-specific fluorescence reporters and facilitates functional genomic analysis of various mutant lines and accessions. We evaluate the purity and quality of the sorted pollen nuclei and analyze the technique's molecular basis. Our results show that both DNA and RNA contents contribute to SYBR Green-activated nucleus sorting and RNA content differences impact on the separation of sperm and vegetative cell nuclei. We demonstrate the power of the approach by sorting wild-type and polyploid mutant sperm and vegetative cell nuclei from mitotic and meiotic mutants, which is not feasible using cell-type-specific transgenic reporters. Our approach should be applicable to pollen nuclei of crop plants and possibly to cell/nucleus types and cell cycle phases of different species containing substantially different amounts of DNA and/or RNA.

The ARID-HMG DNA-binding protein AtHMGB15 is required for pollen tube growth in Arabidopsis thaliana

In flowering plants, male gametes (sperm cells) develop within male gametophytes (pollen grains) and are delivered to female gametes for double fertilization by pollen tubes. Therefore, pollen tube growth is crucial for reproduction. The mechanisms that control pollen tube growth remain poorly understood. In this study, we demonstrated that the ARID-HMG DNA-binding protein AtHMGB15 plays an important role in pollen tube growth. This protein is preferentially expressed in pollen grains and pollen tubes and is localized in the vegetative nuclei of the tricellular pollen grains and pollen tubes. Knocking down AtHMGB15 expression via a Ds insertion caused retarded pollen tube growth, leading to a significant reduction in the seed set. The athmgb15-1 mutation affected the expression of 1686 genes in mature pollen, including those involved in cell wall formation and modification, cell signaling and cellular transport during pollen tube growth. In addition, it was observed that AtHMGB15 binds to DNA in vitro and interacts with the transcription factors AGL66 and AGL104, which are required for pollen maturation and pollen tube growth. These results suggest that AtHMGB15 functions in pollen tube growth through the regulation of gene expression.

Discovery of novel transcripts and gametophytic functions via RNA-seq analysis of maize gametophytic transcriptomes

BACKGROUND

Plant gametophytes play central roles in sexual reproduction. A hallmark of the plant life cycle is that gene expression is required in the haploid gametophytes. Consequently, many mutant phenotypes are expressed in this phase.

RESULTS

We perform a quantitative RNA-seq analysis of embryo sacs, comparator ovules with the embryo sacs removed, mature pollen, and seedlings to assist the identification of gametophyte functions in maize. Expression levels were determined for annotated genes in both gametophytes, and novel transcripts were identified from de novo assembly of RNA-seq reads. Transposon-related transcripts are present in high levels in both gametophytes, suggesting a connection between gamete production and transposon expression in maize not previously identified in any female gametophytes. Two classes of small signaling proteins and several transcription factor gene families are enriched in gametophyte transcriptomes. Expression patterns of maize genes with duplicates in subgenome 1 and subgenome 2 indicate that pollen-expressed genes in subgenome 2 are retained at a higher rate than subgenome 2 genes with other expression patterns. Analysis of available insertion mutant collections shows a statistically significant deficit in insertions in gametophyte-expressed genes.

CONCLUSIONS

This analysis, the first RNA-seq study to compare both gametophytes in a monocot, identifies maize gametophyte functions, gametophyte expression of transposon-related sequences, and unannotated, novel transcripts. Reduced recovery of mutations in gametophyte-expressed genes is supporting evidence for their function in the gametophytes. Expression patterns of extant, duplicated maize genes reveals that selective pressures based on male gametophytic function have likely had a disproportionate effect on plant genomes.

Phosphorylation of a WRKY transcription factor by MAPKs is required for pollen development and function in Arabidopsis

Plant male gametogenesis involves complex and dynamic changes in gene expression. At present, little is known about the transcription factors involved in this process and how their activities are regulated. Here, we show that a pollen-specific transcription factor, WRKY34, and its close homolog, WRKY2, are required for male gametogenesis in Arabidopsis thaliana. When overexpressed using LAT52, a strong pollen-specific promoter, epitope-tagged WRKY34 is temporally phosphorylated by MPK3 and MPK6, two mitogen-activated protein kinases (MAPKs, or MPKs), at early stages in pollen development. During pollen maturation, WRKY34 is dephosphorylated and degraded. Native promoter-driven WRKY34-YFP fusion also follows the same expression pattern at the protein level. WRKY34 functions redundantly with WRKY2 in pollen development, germination, and pollen tube growth. Loss of MPK3/MPK6 phosphorylation sites in WRKY34 compromises the function of WRKY34 in vivo. Epistasis interaction analysis confirmed that MPK6 belongs to the same genetic pathway of WRKY34 and WRKY2. Our study demonstrates the importance of temporal post-translational regulation of WRKY transcription factors in the control of developmental phase transitions in plants.

Pollen development, or male gametogenesis, is a process by which a haploid uninucleate microspore undergoes cell division and specification to form a mature pollen grain containing two sperm cells. The highly defined cell linage makes pollen development an ideal model to understand the regulation of plant cellular development. Pollen development has multiple phases and involves dynamic changes in gene expression, which highlights the importance of transcription factors and their regulatory pathway(s). In this report, we demonstrate that WRKY34 and WRKY2, two closely related WRKY transcription factors in Arabidopsis, play important roles in pollen development. WRKY34 is phosphorylated by MPK3/MPK6, two functionally redundant mitogen-activated protein kinases (MAPKs or MPKs), at early stages in pollen development. Utilizing a combination of genetic, biochemical, and cytological tools, we determined that this MAPK-WRKY signaling module functions at the early stage of pollen development. Loss of function of this pathway reduces pollen viability, and the surviving pollen has poor germination and reduced pollen tube growth, all of which reduce the transmission rate of the mutant pollen. This study discovers a novel stage-specific signaling pathway in pollen development.

Phytoplasma effector SAP54 hijacks plant reproduction by degrading MADS-box proteins and promotes insect colonization in a RAD23-dependent manner

Pathogens that rely upon multiple hosts to complete their life cycles often modify behavior and development of these hosts to coerce them into improving pathogen fitness. However, few studies describe mechanisms underlying host coercion. In this study, we elucidate the mechanism by which an insect-transmitted pathogen of plants alters floral development to convert flowers into vegetative tissues. We find that phytoplasma produce a novel effector protein (SAP54) that interacts with members of the MADS-domain transcription factor (MTF) family, including key regulators SEPALLATA3 and APETALA1, that occupy central positions in the regulation of floral development. SAP54 mediates degradation of MTFs by interacting with proteins of the RADIATION SENSITIVE23 (RAD23) family, eukaryotic proteins that shuttle substrates to the proteasome. Arabidopsis rad23 mutants do not show conversion of flowers into leaf-like tissues in the presence of SAP54 and during phytoplasma infection, emphasizing the importance of RAD23 to the activity of SAP54. Remarkably, plants with SAP54-induced leaf-like flowers are more attractive for colonization by phytoplasma leafhopper vectors and this colonization preference is dependent on RAD23. An effector that targets and suppresses flowering while simultaneously promoting insect herbivore colonization is unprecedented. Moreover, RAD23 proteins have, to our knowledge, no known roles in flower development, nor plant defence mechanisms against insects. Thus SAP54 generates a short circuit between two key pathways of the host to alter development, resulting in sterile plants, and promotes attractiveness of these plants to leafhopper vectors helping the obligate phytoplasmas reproduce and propagate (zombie plants).

MYB97, MYB101 and MYB120 function as male factors that control pollen tube-synergid interaction in Arabidopsis thaliana fertilization

Pollen tube reception involves a pollen tube-synergid interaction that controls the discharge of sperm cells into the embryo sac during plant fertilization. Despite its importance in the sexual reproduction of plants, little is known about the role of gene regulation in this process. We report here that the pollen-expressed transcription factors MYB97, MYB101 and MYB120 probably control genes whose encoded proteins play important roles in Arabidopsis thaliana pollen tube reception. They share a high amino acid sequence identity and are expressed mainly in mature pollen grains and pollen tubes. None of the single or double mutants of these three genes exhibited any visible defective phenotype. Although the myb97 myb101 myb120 triple mutant was not defective in pollen development, pollen germination, pollen tube growth or tube guidance, the pollen tubes of the triple mutants exhibited uncontrolled growth and failed to discharge their sperm cells after entering the embryo sac. In addition, the myb97 myb101 myb120 triple mutation significantly affected the expression of a group of pollen-expressed genes in mature pollen grains. All these results indicate that MYB97, MYB101 and MYB120 participate in pollen tube reception, possibly by controlling the expression of downstream genes.

Pollen tube reception is an important step of fertilization and is controlled by interactions between the pollen tube and synergid. Components of both the pollen tube and synergid are believed to be involved in the process. Several proteins associated with this process have been identified in synergid cells. However, very little is known about the components contributed by the pollen tube. This work identified a group of Arabidopsis pollen-expressed MYB transcription factors, among which at least three members are involved in pollen tube reception. The myb97 myb101 myb120 triple mutation caused overgrowth of the pollen tube into the embryo sac and disrupted sperm cell discharge, leading to failed fertilization. This study provides novel evidence demonstrating that male factors are involved in pollen tube reception.

A parsimonious model of lineage-specific expansion of MADS-box genes in Physcomitrella patens

The MADS-box gene family expanded in the lineage leading to the moss, Physcomitrella patens , mainly as a result of polyploidisations and/or large-scale segmental duplication events and to a lesser extent by tandem duplications. Plant MADS-box genes comprise a large family best known for the roles of type II MIKC (C) genes in floral organogenesis, but also including type II MIKC* genes, some of which have been implicated in male gametophytic development, and type I genes, a few of which are involved in ontogeny of female gametophytes, seeds and embryos. Genome-wide analyses of the MADS-box family in angiosperms have revealed numeric predominance of type I and MIKC (C) genes and cross-species phylogenetic clustering of the Mα, Mβ and Mγ subtypes of type I genes and of 12 major subgroups of MIKC (C) genes. The genome sequence of Physcomitrella patens has facilitated investigation of its full complement of 26 MADS-box genes, including 6 MIKC (C) genes, 11 MIKC* genes, seven type I genes and two pseudogenes. A much higher degree of similarity in sequence and architecture within the MIKC (C) and MIKC* gene subtypes exists in Physcomitrella than in Arabidopsis. Furthermore, MADS-box and K-box sequence is highly conserved between the MIKC (C) and MIKC* subgroups in Physcomitrella. Nine MIKC* genes and two MIKC (C) genes are located in pairs or triplets on individual DNA scaffolds. Phylogenetic gene clustering, gene architectures and gene linkages (directly determined from examination of the genome sequence) underpin a parsimonious model of two tandem duplications and three segmental duplication events, which can account for lineage-specific expansion of the MADS-box gene family in Physcomitrella from 4 members to 26. Two of these segmental duplication events may be indicative of polyploidisations, one of which has been postulated previously.

Genome-wide expression analysis of soybean MADS genes showing potential function in the seed development

The MADS family is an ancient and best-studied transcription factor and plays fundamental roles in almost every developmental process in plants. In the plant evolutionary history, the whole genome duplication (WGD) events are important not only to the plant species evolution, but to expansion of members of the gene families. Soybean as a model legume crop has experience three rounds of WGD events. Members of some MIKC(C) subfamilies, such as SOC, AGL6, SQUA, SVP, AGL17 and DEF/GLO, were expanded after soybean three rounds of WGD events. And some MIKC(C) subfamilies, MIKC* and type I MADS families had experienced faster birth-and-death evolution and their traces before the Glycine WGD event were not found. Transposed duplication played important roles in tandem arrangements among the members of different subfamilies. According to the expression profiles of type I and MIKC paralog pair genes, the fates of MIKC paralog gene pairs were subfunctionalization, and the fates of type I MADS paralog gene pairs were nonfunctionalization. 137 out of 163 MADS genes were close to 186 loci within 2 Mb genomic regions associated with seed-relative QTLs, among which 115 genes expressed during the seed development. Although MIKC(C) genes kept the important and conserved functions of the flower development, most MIKC(C) genes showed potentially essential roles in the seed development as well as the type I MADS.

Functional conservation of MIKC*-Type MADS box genes in Arabidopsis and rice pollen maturation

There are two groups of MADS intervening keratin-like and C-terminal (MIKC)-type MADS box genes, MIKC(C) type and MIKC* type. In seed plants, the MIKC(C) type shows considerable diversity, but the MIKC* type has only two subgroups, P- and S-clade, which show conserved expression in the gametophyte. To examine the functional conservation of MIKC*-type genes, we characterized all three rice (Oryza sativa) MIKC*-type genes. All three genes are specifically expressed late in pollen development. The single knockdown or knockout lines, respectively, of the S-clade MADS62 and MADS63 did not show a mutant phenotype, but lines in which both S-clade genes were affected showed severe defects in pollen maturation and germination, as did knockdown lines of MADS68, the only P-clade gene in rice. The rice MIKC*-type proteins form strong heterodimeric complexes solely with partners from the other subclade; these complexes specifically bind to N10-type C-A-rich-G-boxes in vitro and regulate downstream gene expression by binding to N10-type promoter motifs. The rice MIKC* genes have a much lower degree of functional redundancy than the Arabidopsis thaliana MIKC* genes. Nevertheless, our data indicate that the function of heterodimeric MIKC*-type protein complexes in pollen development has been conserved since the divergence of monocots and eudicots, roughly 150 million years ago.