Analysis of wheat microspore embryogenesis induction by transcriptome and small RNA sequencing using the highly responsive cultivar "Svilena"

Microinjection of the CRISPR/Cas9 editing system through the germ pore of a wheat microspore induces mutations in the target Ms2 gene

BACKGROUND

Microinjection is a direct procedure for delivering various compounds via micropipette into individual cells. Combined with the CRISPR/Cas9 editing technology, it has been used to produce genetically engineered animal cells. However, genetic micromanipulation of intact plant cells has been a relatively unexplored area of research, partly due to the cytological characteristics of these cells. This study aimed to gain insight into the genetic micromanipulation of wheat microspores using microinjection procedures combined with the CRISPR/Cas9 editing system targeting the Ms2 gene.

METHODS AND RESULTS

Microspores were first reprogrammed by starvation and heat shock treatment to make them structurally suitable for microinjection. The large central vacuole was fragmented and the nucleus with cytoplasm was positioned in the center of the cell. This step and an additional maltose gradient provided an adequate source of intact single cells in the three wheat genotypes. The microcapillary was inserted into the cell through the germ pore to deliver a working solution with a fluorescent marker. This procedure was much more efficient and less harmful to the microspore than inserting the microcapillary through the cell wall. The CRISPR/Cas9 binary vectors injected into reprogrammed microspores induced mutations in the target Ms2 gene with deletions ranging from 1 to 16 bp.

CONCLUSIONS

This is the first report of successful genome editing in an intact microspore/wheat cell using the microinjection technique and the CRISPR/Cas9 editing system. The study presented offers a range of molecular and cellular biology tools that can aid in genetic micromanipulation and single-cell analysis.

Metabolic changes during wheat microspore embryogenesis induction using the highly responsive cultivar Svilena

Androgenetically-derived haploids can be obtained by inducing embryogenesis in microspores. Thus, full homozygosity is achieved in a single generation, oppositely to conventional plant breeding programs. Here, the metabolite profile of embryogenic microspores of Triticum aestivum was acquired and integrated with transcriptomic existing data from the same samples in an effort to identify the key metabolic processes occurring during the early stages of microspore embryogenesis. Primary metabolites and transcription profiles were identified at three time points: prior to and immediately following a low temperature pre-treatment given to uninuclear microspores, and after the first nuclear division. This is the first time an integrative -omics analysis is reported in microspore embryogenesis in T. aestivum. The key findings were that the energy produced during the pre-treatment was obtained from the tricarboxylic acid (TCA) cycle and from starch degradation, while starch storage resumed after the first nuclear division. Intermediates of the TCA cycle were highly demanded from a very active amino acid metabolism. The transcription profiles of genes encoding enzymes involved in amino acid synthesis differed from the metabolite profiles. The abundance of glutamine synthetase was correlated with that of glutamine. Cytosolic glutamine synthetase isoform 1 was found predominantly after the nuclear division. Overall, energy production was shown to represent a major component of the de-differentiation process induced by the pre-treatment, supporting a highly active amino acid metabolism.

The effect of phytosulfokine alpha on haploid embryogenesis and gene expression of Brassica napus microspore cultures

Microspore embryogenesis (ME) is the most powerful tool for creating homozygous lines in plant breeding and molecular biology research. It is still based mainly on the reprogramming of microspores by temperature, osmotic and/or nutrient stress. New compounds are being sought that could increase the efficiency of microspore embryogenesis or even induce the formation of haploid embryos from recalcitrant genotypes. Among these, the mitogenic factor phytosulfokine alpha (PSK-α) is promising due to its broad spectrum of activity in vivo and in vitro. The aim of our study was to investigate the effect of PSK-α on haploid embryogenesis from microspores of oilseed rape (Brassica napus L., DH4079), one of the most important oil crops and a model plant for studying the molecular mechanisms controlling embryo formation. We tested different concentrations (0, 0.01, 0.1 and 1 µM) of the peptide and evaluated its effect on microspore viability and embryo regeneration after four weeks of culture. Our results showed a positive correlation between addition of PSK-α and cultured microspore viability and a positive effect also on the number of developed embryos. The analysis of transcriptomes across three time points (day 0, 2 and 4) with or without PSK-α supplementation (15 RNA libraries in total) unveiled differentially expressed genes pivotal in cell division, microspore embryogenesis, and subsequent regeneration. PCA grouped transcriptomes by RNA sampling time, with the first two principal components explaining 56.8% variability. On day 2 with PSK, 45 genes (15 up- and 30 down-regulated) were differentially expressed when PSK-α was added and their number increased to 304 by day 4 (30 up- and 274 down-regulated). PSK, PSKR, and PSI gene expression analysis revealed dynamic patterns, with PSK2 displaying the highest increase and overall expression during microspore culture at days 2 and 4. Despite some variations, only PSK1 showed significant differential expression upon PSK-α addition. Of 16 ME-related molecular markers, 3 and 15 exhibited significant differential expression in PSK-supplemented cultures at days 2 and 4, respectively. Embryo-specific markers predominantly expressed after 4 days of culture, with higher expression in medium without PSK, while on day 0, numerous sporophyte-specific markers were highly expressed.

Transcriptomics integrated with metabolomics reveals the effect of cold stress on rice microspores

BACKGROUND

Microspore culture is one of the important biotechnological tools in plant breeding. The induction of microspore embryogenesis is a critical factor that affects the yield of microspore-derived embryo productions. Cold treatment has been reported to reprogram the gametophytic pathway in various plant species. However, the exact mechanism(s) underlying the effect of cold pre-treatment of floral buds on the efficiency of ME is still not clear.

RESULTS

In this study, the effects of cold stress on the microspore totipotency of rice cultivar Zhonghua 11 were investigated. Our results revealed that a 10-day cold treatment is necessary for microspore embryogenesis initiation. During this period, the survival rate of microspores increased and reached a peak at 7 days post treatment (dpt), before decreasing at 10 dpt. RNA-seq analysis showed that the number of DEGs increased from 3 dpt to 10 dpt, with more downregulated DEGs than upregulated ones at the same time point. GO enrichment analysis showed a shift from 'Response to abiotic stimulus' at 3 dpt to 'Metabolic process' at 7 and 10 dpt, with the most significant category in the cellular component being 'cell wall'. KEGG analysis of the pathways revealed changes during cold treatment. Mass spectrometry was used to evaluate the variations in metabolites at 10 dpt compared to 0 dpt, with more downregulated DEMs being determined in both GC-MS and LC-MS modes. These DEMs were classified into 11 categories, Most of the DEMs belonged to 'lipids and lipid-like molecules'. KEGG analysis of DEMs indicates pathways related to amino acid and nucleotide metabolism being upregulated and those related to carbohydrate metabolism being downregulated. An integration analysis of transcriptomics and metabolomics showed that most pathways belonged to 'Amino acid metabolism' and 'Carbohydrate metabolism'. Four DEMs were identified in the interaction network, with stearidonic acid involving in the most correlations, suggesting the potential role in microspore totipotency.

CONCLUSIONS

Our findings exhibited the molecular events occurring during stress-induced rice microspore. Pathways related to 'Amino acid metabolism' and 'Carbohydrate metabolism' may play important roles in rice microspore totipotency. Stearidonic acid was identified, which may participate in the initiation of microspore embryogenesis.

Microspore embryogenesis induction by mannitol and TSA results in a complex regulation of epigenetic dynamics and gene expression in bread wheat

Reprogramming of microspores development towards embryogenesis mediated by stress treatment constitutes the basis of doubled haploid production. Recently, compounds that alter histone post-translational modifications (PTMs) have been reported to enhance microspore embryogenesis (ME), by altering histones acetylation or methylation. However, epigenetic mechanisms underlying ME induction efficiency are poorly understood. In this study, the epigenetic dynamics and the expression of genes associated with histone PTMs and ME induction were studied in two bread wheat cultivars with different ME response. Microspores isolated at 0, 3 and 5 days, treated with 0.7M mannitol (MAN) and 0.7M mannitol plus 0.4µM trichostatin A (TSA), which induced ME more efficiently, were analyzed. An additional control of gametophytic development was included. Microspores epigenetic state at the onset of ME induction was distinctive between cultivars by the ratio of H3 variants and their acetylated forms, the localization and percentage of labeled microspores with H3K9ac, H4K5ac, H4K16ac, H3K9me2 and H3K27me3, and the expression of genes related to pollen development. These results indicated that microspores of the high responding cultivar could be at a less advanced stage in pollen development. MAN and TSA resulted in a hyperacetylation of H3.2, with a greater effect of TSA. Histone PTMs were differentially affected by both treatments, with acetylation being most concerned. The effect of TSA was observed in the H4K5ac localization pattern at 3dT in the mid-low responding cultivar. Three gene networks linked to ME response were identified. TaHDT1, TaHAG2, TaYAO, TaNFD6-A, TabZIPF1 and TaAGO802-B, associated with pollen development, were down-regulated. TaHDA15, TaHAG3, TaHAM, TaYUC11D, Ta-2B-LBD16 TaMS1 and TaDRM3 constituted a network implicated in morphological changes by auxin signaling and cell wall modification up-regulated at 3dT. The last network included TaHDA18, TaHAC1, TaHAC4, TaABI5, TaATG18fD, TaSDG1a-7A and was related to ABA and ethylene hormone signaling pathways, DNA methylation and autophagy processes, reaching the highest expression at 5dT. The results indicated that TSA mainly modified the regulation of genes related to pollen and auxin signaling. This study represents a breakthrough in identifying the epigenetic dynamics and the molecular mechanisms governing ME induction efficiency, with relevance to recalcitrant wheat genotypes and other crops.

Phospholipase C is a novel regulator at the early stages of microspore embryogenesis in Nicotiana tabacum

Microspore transfers the developmental fate into embryogenesis in vitro regulated by determinant factors of stress-induced. However, the key regulators of microspore embryogenesis (ME) are still largely undiscovered to reveal the mechanism of cell fate transition. Here, we report that Phospholipase C (PLC) is involved at the early stages of ME in Nicotiana tabacum. NtPLC2/3/4 are expressed at the initial stages of ME. The expression levels of NtPLC2/3 are transient activated after 3 days in culture, while the expression level of NtPLC4 maintains relatively stable. Inhibition of PLCs induces the decrease in NtPLC2/3/4 expression level and decline of ME yield. We confirm that lipids in microspore are degraded and then re-accumulate at first embryonic division stage. Inhibition of PLCs suppresses the lipids metabolism at the early stages of ME. Thus, we propose that PLCs-mediated lipid metabolism is a novel regulator at the early stages of ME.

Comprehensive, integrative genomic analysis of microRNA expression profiles in different tissues of two wheat cultivars with different traits

Wheat is one of the most important food sources on Earth. MicroRNAs (miRNA) play important roles in wheat productivity. To identify wheat miRNAs, we constructed and sequenced sRNA libraries from leaves and roots of two wheat cultivars (RAC875 and Kukri) with many different traits. Given that available miRNA wheat complement in the plant-specific database PmiREN ( https://pmiren.com ) does not include root tissues and root-associated miRNAs might thus be missing, we performed first the prediction of novel miRNAs using the sRNAbench tool. We found a total of 150 putatively novel miRNA genes with expression of both arms from 289 unique mature sequences and nearly 30% of all miRNA reads in roots corresponded to novel miRNAs. In contrast, this figure in leaves dropped to under 3%, confirming the undersampling of roots in the complement of known miRNAs. By using 120 publicly available wheat datasets, 598 Zea mays small RNA libraries, 64 plant species genomes, wheat degradome library, and functional enrichment analysis, a subset of novel miRNAs were confirmed as bona-fide miRNAs. Of the total 605 miRNAs identified in this study inclusive of 316 known miRNAs, 528 miRNAs were shared by both cultivars, 429 miRNAs were shared by both root tissues and 329 miRNAs were shared by both leaf tissues. In addition, 32 miRNAs were specific to Kukri while 45 miRNAs were specific to RAC875. These miRNAs had diverse functions, such as regulation of gene transcription, protein translation, energy metabolism, and cell cycle progression. Our data provide a genome-wide miRNA expression profile in these two wheat cultivars and help functional studies of wheat genomics.

Androgenesis-Based Doubled Haploidy: Past, Present, and Future Perspectives

Androgenesis, which entails cell fate redirection within the microgametophyte, is employed widely for genetic gain in plant breeding programs. Moreover, androgenesis-responsive species provide tractable systems for studying cell cycle regulation, meiotic recombination, and apozygotic embryogenesis within plant cells. Past research on androgenesis has focused on protocol development with emphasis on temperature pretreatments of donor plants or floral buds, and tissue culture optimization because androgenesis has different nutritional requirements than somatic embryogenesis. Protocol development for new species and genotypes within responsive species continues to the present day, but slowly. There is more focus presently on understanding how protocols work in order to extend them to additional genotypes and species. Transcriptomic and epigenetic analyses of induced microspores have revealed some of the cellular and molecular responses required for or associated with androgenesis. For example, microRNAs appear to regulate early microspore responses to external stimuli; trichostatin-A, a histone deacetylase inhibitor, acts as an epigenetic additive; ά-phytosulfokine, a five amino acid sulfated peptide, promotes androgenesis in some species. Additionally, present work on gene transfer and genome editing in microspores suggest that future endeavors will likely incorporate greater precision with the genetic composition of microspores used in doubled haploid breeding, thus likely to realize a greater impact on crop improvement. In this review, we evaluate basic breeding applications of androgenesis, explore the utility of genomics and gene editing technologies for protocol development, and provide considerations to overcome genotype specificity and morphogenic recalcitrance in non-model plant systems.

Occurrence of albinism during wheat androgenesis is correlated with repression of the key genes required for proper chloroplast biogenesis

The phenomenon of albinism in wheat androgenesis is linked to the transcriptional repression of specific genes involved in chloroplast biogenesis during the first weeks of in vitro culture. Isolated microspore culture is widely used to accelerate breeding programs and produce new cultivars. However, in cereals and particularly in wheat, the use of this technique is limited due to the high proportion of regenerated albino plantlets. The causes and mechanisms leading to the formation of albino plantlets in wheat remain largely unknown and, to date, no concrete solution has been found to make it possible to overcome this barrier. We performed a molecular study of proplastid-to-chloroplast differentiation within wheat microspore cultures by analyzing the expression of 20 genes specifically involved in chloroplast biogenesis. Their expression levels were compared between two wheat genotypes that exhibit differential capacities to regenerate green plantlets, i.e., Pavon and Paledor, which produce high and low rates of green plants, respectively. We observed that chloroplast biogenesis within wheat microspores was affected as of the very early stages of the androgenesis process. A successful transition from a NEP- to a PEP-dependent transcription during early plastid development was found to be strongly correlated with the formation of green plantlets, while failure of this transition was strongly correlated with the regeneration of albino plantlets. The very low expression of plastid-encoded 16S and 23S rRNAs within plastids of the recalcitrant genotype Paledor suggests a low translation activity in albino plastids. Furthermore, a delay in the activation of the transcription of nuclear encoded key genes like GLK1 related to chloroplast biogenesis was observed in multicellular structures and pro-embryos of the genotype Paledor. These data help to understand the phenomenon of albinism in wheat androgenesis, which appears to be linked to the transcriptional activation of specific genes involved in the initial steps of chloroplast biogenesis that occurs between days 7 and 21 of in vitro culture.

Comparative proteomic analysis provides new insights into regulation of microspore embryogenesis induction in winter triticale (× Triticosecale Wittm.) after 5-azacytidine treatment

Effective microspore embryogenesis (ME) requires substantial modifications in gene expression pattern, followed by changes in the cell proteome and its metabolism. Recent studies have awakened also interest in the role of epigenetic factors in microspore de-differentiation and reprogramming. Therefore, demethylating agent (2.5-10 μM 5-azacytidine, AC) together with low temperature (3 weeks at 4 °C) were used as ME-inducing tiller treatment in two doubled haploid (DH) lines of triticale and its effect was analyzed in respect of anther protein profiles, expression of selected genes (TAPETUM DETERMINANT1 (TaTPD1-like), SOMATIC EMBRYOGENESIS RECEPTOR KINASE 2 (SERK2) and GLUTATHIONE S-TRANSFERASE (GSTF2)) and ME efficiency. Tiller treatment with 5.0 µM AC was the most effective in ME induction; it was associated with (1) suppression of intensive anabolic processes-mainly photosynthesis and light-dependent reactions, (2) transition to effective catabolism and mobilization of carbohydrate reserve to meet the high energy demand of cells during microspore reprograming and (3) effective defense against stress-inducing treatment, i.e. protection of proper folding during protein biosynthesis and effective degradation of dysfunctional or damaged proteins. Additionally, 5.0 µM AC enhanced the expression of all genes previously identified as being associated with embryogenic potential of microspores (TaTPD1-like, SERK and GSTF2).

Proteins, Small Peptides and Other Signaling Molecules Identified as Inconspicuous but Possibly Important Players in Microspores Reprogramming Toward Embryogenesis

In this review, we describe and integrate the latest knowledge on the signaling role of proteins and peptides in the stress-induced microspore embryogenesis (ME) in some crop plants with agricultural importance (i.e., oilseed rape, tobacco, barley, wheat, rice, triticale, rye). Based on the results received from the most advanced omix analyses, we have selected some inconspicuous but possibly important players in microspores reprogramming toward embryogenic development. We provide an overview of the roles and downstream effect of stress-related proteins (e.g., β-1,3-glucanases, chitinases) and small signaling peptides, especially cysteine—(e.g., glutathione, γ-thionins, rapid alkalinization factor, lipid transfer, phytosulfokine) and glycine-rich peptides and other proteins (e.g., fasciclin-like arabinogalactan protein) on acclimation ability of microspores and the cell wall reconstruction in a context of ME induction and haploids/doubled haploids (DHs) production. Application of these molecules, stimulating the induction and proper development of embryo-like structures and green plant regeneration, brings significant improvement of the effectiveness of DHs procedures and could result in its wider incorporation on a commercial scale. Recent advances in the design and construction of synthetic peptides–mainly cysteine-rich peptides and their derivatives–have accelerated the development of new DNA-free genome-editing techniques. These new systems are evolving incredibly fast and soon will find application in many areas of plant science and breeding.

Reinvestigation of THOUSAND-GRAIN WEIGHT 6 grain weight genes in wheat and rice indicates a role in pollen development rather than regulation of auxin content in grains

Phylogenetic and expression analyses of grain weight genes TaTGW6 and OsTGW6 and investigation of substrate availability indicate TGW6 does not regulate auxin content of grains but may affect pollen development. The THOUSAND-GRAIN WEIGHT 6 genes (TaTGW6 and OsTGW6) are reported to result in larger grains of wheat and rice by reducing production of indole-3-acetic acid (IAA) in developing grains. However, a critical comparison of data on TaTGW6 and OsTGW6 with other reports on IAA synthesis in cereal grains requires that this hypothesis be reinvestigated. Here, we show that TaTGW6 and OsTGW6 are members of a large gene family that has undergone major, lineage-specific gene expansion. Wheat has nine genes, and rice three genes encoding proteins with more than 80% amino acid identity with TGW6, making it difficult to envisage how a single inactive allele could have a major effect on IAA levels in grains. In our study, we show that neither TaTGW6 nor OsTGW6 is expressed in developing grains. Instead, both genes and their close homologues are exclusively expressed in pre-emergent inflorescences; TaTGW6 is expressed particularly in microspores prior to mitosis. This evidence, combined with our observation that developing wheat grains have undetectable levels of ester IAA in comparison to free IAA and do not express an IAA-glucose synthase suggests that TaTGW6 and OsTGW6 do not regulate grain size via the hydrolysis of IAA-glucose. Instead, their similarity to rice strictosidine synthase-like (OsSTRL2) suggests they play a key role in pollen development.

MicroRNA Omics Analysis of Camellia sinesis Pollen Tubes in Response to Low-Temperature and Nitric Oxide

Nitric oxide (NO) as a momentous signal molecule participates in plant reproductive development and responds to various abiotic stresses. Here, the inhibitory effects of the NO-dominated signal network on the pollen tube growth of Camellia sinensis under low temperature (LT) were studied by microRNA (miRNA) omics analysis. The results showed that 77 and 71 differentially expressed miRNAs (DEMs) were induced by LT and NO treatment, respectively. Gene ontology (GO) analysis showed that DEM target genes related to microtubules and actin were enriched uniquely under LT treatment, while DEM target genes related to redox process were enriched uniquely under NO treatment. In addition, the target genes of miRNA co-regulated by LT and NO are only located on the cell membrane and cell wall, and most of them are enriched in metal ion binding and/or transport and cell wall organization. Furthermore, DEM and its target genes related to metal ion binding/transport, redox process, actin, cell wall organization and carbohydrate metabolism were identified and quantified by functional analysis and qRT-PCR. In conclusion, miRNA omics analysis provides a complex signal network regulated by NO-mediated miRNA, which changes cell structure and component distribution by adjusting Ca²⁺ gradient, thus affecting the polar growth of the C. sinensis pollen tube tip under LT.

Enzymatic characterization of ancestral/group-IV clade xyloglucan endotransglycosylase/hydrolase enzymes reveals broad substrate specificities

Xyloglucan endotransglycosylase/hydrolase (XTH) enzymes play important roles in cell wall remodelling. Although previous studies have shown a pathway of evolution for XTH genes from bacterial licheninases, through plant endoglucanases (EG16), the order of development within the phylogenetic clades of true XTHs is yet to be elucidated. In addition, recent studies have revealed interesting and potentially useful patterns of transglycosylation beyond the standard xyloglucan-xyloglucan donor/acceptor substrate activities. To study evolutionary relationships and to search for enzymes with useful broad substrate specificities, genes from the 'ancestral' XTH clade of two monocots, Brachypodium distachyon and Triticum aestivum, and two eudicots, Arabidopsis thaliana and Populus tremula, were investigated. Specific activities of the heterologously produced enzymes showed remarkably broad substrate specificities. All the enzymes studied had high activity with the cellulose analogue HEC (hydroxyethyl cellulose) as well as with mixed-link β-glucan as donor substrates, when compared with the standard xyloglucan. Even more surprising was the wide range of acceptor substrates that these enzymes were able to catalyse reactions with, opening a broad range of possible roles for these enzymes, both within plants and in industrial, pharmaceutical and medical fields. Genome screening and expression analyses unexpectedly revealed that genes from this clade were found only in angiosperm genomes and were predominantly or solely expressed in reproductive tissues. We therefore posit that this phylogenetic group is significantly different and should be renamed as the group-IV clade.

β-1,3-Glucanases and chitinases participate in the stress-related defence mechanisms that are possibly connected with modulation of arabinogalactan proteins (AGP) required for the androgenesis initiation in rye (Secale cereale L.)

This work presents the biochemical, cytochemical and molecular studies on two groups of PR proteins, β-1,3-glucanases and chitinases, and the arabinogalactan proteins (AGP) during the early stages of androgenesis induction in two breeding lines of rye (Secale cereale L.) with different androgenic potential. The process of androgenesis was initiated by tillers pre-treatments with low temperature, mannitol and/or reduced glutathione and resulted in microspores reprogramming and formation of androgenic structures what was associated with high activity of β-1,3-glucanases and chitinases. Some isoforms of β-1,3-glucanases, namely several acidic isoforms of about 26 kDa; appeared to be anther specific. Chitinases were well represented but were less variable. RT-qPCR revealed that the cold-responsive chitinase genes Chit1 and Chit2 were expressed at a lower level in the microspores and whole anthers while the cold-responsive Glu2 and Glu3 were not active. The stress pre-treatments modifications promoted the AGP accumulation. An apparent dominance of some AGP epitopes (LM2, JIM4 and JIM14) was detected in the androgenesis-responsive rye line. An abundant JIM13 epitopes in the vesicles and inner cell walls of the microspores and in the cell walls of the anther cell layers appeared to be the most specific for embryogenesis.

The commitment of barley microspores into embryogenesis correlates with miRNA-directed regulation of members of the SPL, GRF and HD-ZIPIII transcription factor families

Microspore embryogenesis is a model for developmental plasticity and cell fate decisions. To investigate the role of miRNAs in this development, we sequenced sRNAs and the degradome of barley microspores collected prior to (day 0) and after (days 2 and 5) the application of a stress treatment known to induce embryogenesis. Microspores isolated at these timepoints were uniform in both appearance and in their complements of sRNAs. We detected 68 miRNAs in microspores. The abundance of 51 of these miRNAs differed significantly during microspore development. One group of miRNAs was induced when the stress treatment was applied, prior to being repressed when microspores transitioned to embryogenesis. Another group of miRNAs were up-regulated in day-2 microspores and their abundance remained stable or increased in day-5 microspores, a timepoint at which the first clear indications of the transition toward embryogenesis were visible. Collectively, these miRNAs might play a role in the modulation of the stress response, the repression of gametic development, and/or the gain of embryogenic potential. A degradome analysis allowed us to validate the role of miRNAs in regulating 41 specific transcripts. We showed that the transition of microspores toward the embryogenesis pathway involves miRNA-directed regulation of members of the ARF, SPL, GRF, and HD-ZIPIII transcription factor families. We noted that 41.5% of these targets were shared between day-2 and day-5 microspores while 26.8% were unique to day-5 microspores. The former set may act to disrupt transcripts involved in pollen development while the latter set may drive the commitment to embryogenesis.

Differential Expression Profiling Reveals Stress-Induced Cell Fate Divergence in Soybean Microspores

Stress-induced microspore embryogenesis is a widely employed method to achieve homozygosity in plant breeding programs. However, the molecular mechanisms that govern gametophyte de- and redifferentiation are understood poorly. In this study, RNA-Seq was used to evaluate global changes across the microspore transcriptome of soybean (Glycine max [L.] Merrill) as a consequence of pretreatment low-temperature stress. Expression analysis revealed more than 20,000 differentially expressed genes between treated and control microspore populations. Functional enrichment illustrated that many of these genes (e.g., those encoding heat shock proteins and cytochrome P450s) were upregulated to maintain cellular homeostasis through the mitigation of oxidative damage. Moreover, transcripts corresponding to saccharide metabolism, vacuolar transport, and other pollen-related developmental processes were drastically downregulated among treated microspores. Temperature stress also triggered cell wall modification and cell proliferation-characteristics that implied putative commitment to an embryonic pathway. These findings collectively demonstrate that pretreatment cold stress induces soybean microspore reprogramming through suppression of the gametophytic program while concomitantly driving sporophytic development.

Trichostatin A Affects Developmental Reprogramming of Bread Wheat Microspores towards an Embryogenic Route

Microspores can be developmentally reprogrammed by the application of different stress treatments to initiate an embryogenic pathway leading to the production of doubled haploid (DH) plants. Epigenetic modifications are involved in cell reprogramming and totipotency in response to stress. To increase microspore embryogenesis (ME) efficiency in bread wheat, the effect of the histone deacetylase inhibitor trichostatin A (TSA) has been examined in two cultivars of wheat with different microspore embryogenesis response. Diverse strategies were assayed using 0-0.4 µM TSA as a single induction treatment and after or simultaneously with cold or mannitol stresses. The highest efficiency was achieved when 0.4 µM TSA was applied to anthers for 5 days simultaneously with a 0.7 M mannitol treatment, producing a four times greater number of green DH plants than mannitol. Ultrastructural studies by transmission electron microscopy indicated that mannitol with TSA and mannitol treatments induced similar morphological changes in early stages of microspore reprogramming, although TSA increased the number of microspores with 'star-like' morphology and symmetric divisions. The effect of TSA on the transcript level of four ME marker genes indicated that the early signaling pathways in ME, involving the TaTDP1 and TAA1b genes, may be mediated by changes in acetylation patterns of histones and/or other proteins.

Research Tools for the Functional Genomics of Plant miRNAs During Zygotic and Somatic Embryogenesis

During early plant embryogenesis, some of the most fundamental decisions on fate and identity are taken making it a fascinating process to study. It is no surprise that higher plant embryogenesis was intensively analysed during the last century, while somatic embryogenesis is probably the most studied regeneration model. Encoded by the MIRNA, short, single-stranded, non-coding miRNAs, are commonly present in all Eukaryotic genomes and are involved in the regulation of the gene expression during the essential developmental processes such as plant morphogenesis, hormone signaling, and developmental phase transition. During the last few years dedicated to miRNAs, analytical methods and tools have been developed, which have afforded new opportunities in functional analyses of plant miRNAs, including (i) databases for in silico analysis; (ii) miRNAs detection and expression approaches; (iii) reporter and sensor lines for a spatio-temporal analysis of the miRNA-target interactions; (iv) in situ hybridisation protocols; (v) artificial miRNAs; (vi) MIM and STTM lines to inhibit miRNA activity, and (vii) the target genes resistant to miRNA. Here, we attempted to summarise the toolbox for functional analysis of miRNAs during plant embryogenesis. In addition to characterising the described tools/methods, examples of the applications have been presented.

Effect of Zearalenone and Hormone Regulators on Microspore Embryogenesis in Anther Culture of Wheat

The purpose of this work was to assess the impact of zearalenone (ZEN) and selected hormone regulators on the effectiveness of microspore embryogenesis in anther culture of wheat. The plant material comprised F₁ hybrids of winter and spring wheat. Six combinations of media inducing microspore proliferation and formation of embryogenic structures were investigated: two combinations of growth regulators (D - 2,4-D + dicamba, K - 2,4-D + kinetin), each with three ZEN concentrations (0 mL/L, 0.1 mL/L, 0.2 mL/L). A significant increase in microspore embryogenesis effectiveness on media with the addition of ZEN was observed both at the stages of its induction and the formation of green plants in some genotypes. In case of both combinations of growth regulators, an increased concentration of ZEN resulted in more effective induction of microspore embryogenesis. The most effective induction medium was the D medium supplemented with 0.2 mL/L ZEN. As a result of the use of zearalenone together with two combinations of growth regulators, all genotypes tested produced androgenic structures, which indicates the breakdown of genotypic recalcitrant in the analysed hybrids. In addition, green plants were obtained from 18 out of 19 tested hybrids. The addition of ZEN to the medium did not affect the number of regenerated albino plants nor the number of spontaneous genome doublings proportion.

Chemically-induced DNA de-methylation alters the effectiveness of microspore embryogenesis in triticale

Microspores exposed to some stress factors may display cell totipotency and could be reprogrammed towards embryogenic development. Plant breeding and genetic engineering widely use haploids/doubled haploids (DHs) derived from in vitro-cultured microspores, but the mechanism of this process remains poorly understood. Recently published data suggest that microspore embryogenesis (ME) is accompanied by changes in DNA methylation and chromatin reorganization. Here, we used two triticale DH lines (DH19 and DH28), significantly different with respect to embryogenic potential. To change DNA methylation levels, we applied two cytosine-analogs: 5-azacytidine (AC) and 2'-deoxy-5-azacytidine (DAC) treatments. We found that chemically-induced DNA demethylation caused chromatin relaxation and dysregulation of marker genes (TaTPD1-like, GSTF2, GSTA2, CHI3, Tad1, TaNF-YA7, SERK2, TaME1) related to ME. Both drugs showed significant cytotoxicity in a dose-dependent manner. We noticed that lines varied in terms of overall DNA methylation levels and responded in a different way to hypomethylation caused by the drugs. DH19 (low embryogenic) after inhibitors treatment, showed higher microspore viability, but its recalcitrancy was not overcome. For highly embryogenic DH28, we noted significantly higher effectiveness of embryo-like structure production and plant regeneration. In summary, our study provides new insight into the role of DNA methylation in ME initiation. They suggest potential benefits resulting from the utilization of epigenetic inhibitors to improve the process of DHs production.

De novo assembly of wheat root transcriptomes and transcriptional signature of longitudinal differentiation

Hidden underground, root systems constitute an important part of the plant for its development, nourishment and sensing the soil environment around it, but we know very little about its genetic regulation in crop plants like wheat. In the present study, we de novo assembled the root transcriptomes in reference cultivar Chinese Spring from RNA-seq reads generated by the 454-GS-FLX and HiSeq platforms. The FLX reads were assembled into 24,986 transcripts with completeness of 54.84%, and the HiSeq reads were assembled into 91,543 high-confidence protein-coding transcripts, 2,404 low-confidence protein-coding transcripts, and 13,181 non-coding transcripts with the completeness of >90%. Combining the FLX and HiSeq assemblies, we assembled a root transcriptome of 92,335 ORF-containing transcripts. Approximately 7% of the coding transcripts and ~2% non-coding transcripts are not present in the current wheat genome assembly. Functional annotation of both assemblies showed similar gene ontology patterns and that ~7% coding and >5% non-coding transcripts are root-specific. Transcription quantification identified 1,728 differentially expressed transcripts between root tips and maturation zone, and functional annotation of these transcripts captured a transcriptional signature of longitudinal development of wheat root. With the transcriptomic resources developed, this study provided the first view of wheat root transcriptome under different developmental zones and laid a foundation for molecular studies of wheat root development and growth using a reverse genetic approach.

Transcriptomics analyses reveal the molecular roadmap and long non-coding RNA landscape of sperm cell lineage development

Sperm cell (SC) lineage development from the haploid microspore to SCs represents a unique biological process in which the microspore generates a larger vegetative cell (VC) and a smaller generative cell (GC) enclosed in the VC, then the GC further develops to functionally specified SCs in the VC for double fertilization. Understanding the mechanisms of SC lineage development remains a critical goal in plant biology. We isolated individual cells of the three cell types, and characterized the genome-wide atlas of long non-coding (lnc) RNAs and mRNAs of haploid SC lineage cells. Sperm cell lineage development involves global repression of genes for pluripotency, somatic development and metabolism following asymmetric microspore division and coordinated upregulation of GC/SC preferential genes. This process is accompanied by progressive loss of the active marks H3K4me3 and H3K9ac, and accumulation of the repressive methylation mark H3K9. The SC lineage has a higher ratio of lncRNAs to mRNAs and preferentially expresses a larger percentage of lncRNAs than does the non-SC lineage. A co-expression network showed that the largest set of lncRNAs in these nodes, with more than 100 links, are GC-preferential, and a small proportion of lncRNAs co-express with their neighboring genes. Single molecular fluorescence in situ hybridization showed that several candidate genes may be markers distinguishing the three cell types of the SC lineage. Our findings reveal the molecular programming and potential roles of lncRNAs in SC lineage development.

Insights into the cotton anther development through association analysis of transcriptomic and small RNA sequencing

BACKGROUND

Plant anther development is a systematic and complex process precisely controlled by genes. Regulation genes and their regulatory mechanisms for this process remain elusive. In contrast to numerous researches on anther development with respect to mRNAs or miRNAs in many crops, the association analysis combining both omics has not been reported on cotton anther.

RESULTS

In this study, the molecular mechanism of cotton anther development was investigated with the employment of association analysis of transcriptome and small RNA sequencing during the predefined four stages of cotton anther development, sporogenuous cell proliferation (SCP), meiotic phase (MP), microspore release period (MRP) and pollen maturity (PM). Analysis revealed that the differentially expressed genes are increasingly recruited along with the developmental progress. Expression of functional genes differed significantly among developmental stages. The genes related with cell cycle, progesterone-mediated oocyte maturation, and meiosis are predominantly expressed at the early stage of anther development (SCP and MP), and the expression of genes involved in energy metabolism, flavonoid biosynthesis, axon guidance and phospholipase D signaling pathways is mainly enriched at the late stage of anther development (MRP and PM). Analysis of expression patterns revealed that there was the largest number of differentially expressed genes in the MP and the expression profiles of differentially expressed genes were significantly increased, which implied the importance of MP in the entire anther development cycle. In addition, prediction and analysis of miRNA targeted genes suggested that miRNAs play important roles in anther development. The miRNAs ghr-miR393, Dt_chr12_6065 and At_chr9_3080 participated in cell cycle, carbohydrate metabolism and auxin anabolism through the target genes, respectively, to achieve the regulation of anther development.

CONCLUSIONS

Through the association analysis of mRNA and miRNA, our work gives a better understanding of the preferentially expressed genes and regulation in different developmental stages of cotton anther and the importance of meiotic phase, and also the involvement of miRNAs in precise regulation for this process, which would be valuable for clarifying the mechanism of plant anther development in response to internal and external environments.

Genetic Loci Governing Androgenic Capacity in Perennial Ryegrass (Lolium perenne L.)

Immature pollen can be induced to switch developmental pathways from gametogenesis to embryogenesis and subsequently regenerate into homozygous, diploid plants. Such androgenic production of doubled haploids is particularly useful for species where inbreeding is hampered by effective self-incompatibility systems. Therefore, increasing the generally low androgenic capacity of perennial ryegrass (Lolium perenne L.) germplasm would enable the efficient production of homozygous plant material, so that a more effective exploitation of heterosis through hybrid breeding schemes can be realized. Here, we present the results of a genome-wide association study in a heterozygous, multiparental population of perennial ryegrass (n = 391) segregating for androgenic capacity. Genotyping-by-sequencing was used to interrogate gene- dense genomic regions and revealed over 1,100 polymorphic sites. Between one and 10 quantitative trait loci (QTL) were identified for anther response, embryo and total plant production, green and albino plant production and regeneration. Most traits were under polygenic control, although a major QTL on linkage group 5 was associated with green plant regeneration. Distinct genetic factors seem to affect green and albino plant recovery. Two intriguing candidate genes, encoding chromatin binding domains of the developmental phase transition regulator, Polycomb Repressive Complex 2, were identified. Our results shed the first light on the molecular mechanisms behind perennial ryegrass microspore embryogenesis and enable marker-assisted introgression of androgenic capacity into recalcitrant germplasm of this forage crop of global significance.

Differential Expression Profiling of Microspores During the Early Stages of Isolated Microspore Culture Using the Responsive Barley Cultivar Gobernadora

In barley, it is possible to induce embryogenesis in the haploid and uninucleate microspore to obtain a diploid plant that is perfectly homozygous. To change developmental fates in this fashion, microspores need to engage in cellular de-differentiation, interrupting the pollen formation, and restore totipotency prior to engaging in embryogenesis. In this work, we used the barley cultivar Gobernadora to characterize the transcriptome of microspores prior to (day 0) and immediately after (days 2 and 5) the application of a stress pretreatment. A deep RNA-seq analysis revealed that microspores at these three time points exhibit a transcriptome of ∼14k genes, ∼90% of which were shared. An expression analysis identified a total of 3,382 differentially expressed genes (DEGs); of these, 2,155 and 2,281 DEGs were respectively identified when contrasting expression at days 0 and 2 and at days 2 and 5. These define 8 expression profiles in which DEGs share a common up- or down-regulation at these time points. Up-regulation of numerous glutathione S-transferase and heat shock protein genes as well as down-regulation of ribosomal subunit protein genes was observed between days 0 and 2. The transition from microspores to developing embryos (days 2 vs. 5) was marked by the induction of transcription factor genes known to play important roles in early embryogenesis, numerous genes involved in hormone biosynthesis and plant hormonal signal transduction in addition to genes involved in secondary metabolism. This work sheds light on transcriptional changes accompanying an important developmental shift and provides candidate biomarkers for embryogenesis in barley.

Dynamics of post-translationally modified histones during barley pollen embryogenesis in the presence or absence of the epi-drug trichostatin A

Improving pollen embryogenesis. Despite the agro-economic importance of pollen embryogenesis, the mechanisms underlying this process are still poorly understood. We describe the dynamics of chromatin modifications (histones H3K4me2, H3K9ac, H3K9me2, and H3K27me3) and chromatin marks (RNA polymerase II CDC phospho-Ser5, and CENH3) during barley pollen embryogenesis. Immunolabeling results show that, in reaction to stress, immature pollen rapidly starts reorganizing several important chromatin modifications indicative of a change in cell fate. This new chromatin modification pattern was accomplished within 24 h from whereon it remained unaltered during subsequent mitotic activity. This indicates that cell fate transition, the central element of pollen embryogenesis, is completed early on during the induction process. Application of the histone deacetylase inhibitor trichostatin A stimulated pollen embryogenesis when used on pollen with a gametophytic style chromatin pattern. However, when this drug was administered to embryogenic pollen, the chromatin markers reversed toward a gametophytic profile, embryogenesis was halted and all pollen invariably died.