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

Whole-genome landscape of histone H3K4me3 modification during sperm cell lineage development in tomato

BACKGROUND

During male gametogenesis of flowering plants, sperm cell lineage (microspores, generative cells, and sperm cells) differentiated from somatic cells and acquired different cell fates. Trimethylation of histone H3 on lysine 4 (H3K4me3) epigenetically contributes to this process, however, it remained unclear how H3K4me3 influences the gene expression in each cell type. Here, we conducted chromatin immunoprecipitation sequencing (ChIP-seq) to obtain a genome-wide landscape of H3K4me3 during sperm cell lineage development in tomato (Solanum lycopersicum).

RESULTS

We show that H3K4me3 peaks were mainly enriched in the promoter regions, and intergenic H3K4me3 peaks expanded as sperm cell lineage differentiated from somatic cells. H3K4me3 was generally positively associated with transcript abundance and served as a better indicator of gene expression in somatic and vegetative cells, compared to sperm cell lineage. H3K4me3 was mutually exclusive with DNA methylation at 3' proximal of the transcription start sites. The microspore maintained the H3K4me3 features of somatic cells, while generative cells and sperm cells shared an almost identical H3K4me3 pattern which differed from that of the vegetative cell. After microspore division, significant loss of H3K4me3 in genes related to brassinosteroid and cytokinin signaling was observed in generative cells and vegetative cells, respectively.

CONCLUSIONS

Our results suggest the asymmetric division of the microspore significantly reshapes the genome-wide distribution of H3K4me3. Selective loss of H3K4me3 in genes related to hormone signaling may contribute to functional differentiation of sperm cell lineage. This work provides new resource data for the epigenetic studies of gametogenesis in plants.

A repertoire of intronic lariat RNAs reveals tissue-specific regulation and target mimicry potential in plants

Lariat RNA is concomitantly produced by excised intron during RNA splicing, which is usually debranched by DBR1, an RNA debranching enzyme. However, increasing evidence showed that some lariat RNA could escape debranching. Little is known about how and why these lariat RNAs could be retained. By comparing the atlas of lariat RNAs between the non-dividing cell (mature pollen) and three actively dividing tissues (young shoot apex, young seeds, and young roots), we identified hundreds to thousands of lariat RNA naturally retained in each tissue, and the incidence of lariat RNA retention is much less in shoot apex while much more in pollen. Many lariat RNAs derived from the same intron or different lariat RNAs from the same pre-mRNA could be retained in one tissue while degraded in the other tissues. By deciphering lariat RNA sequences, we identified an AG-rich (RAAAAVAAAR) motif and a UC-rich (UCUCUYUCUC) motif for pollen-specific and the other three tissues-retained lariat RNAs, respectively. Reconstitution of the pollen-specific AG-rich motif indeed enhanced lariat RNA retention in plants. Biologically, hundreds of lariat RNAs harbored miRNA binding sites, and dual-luciferase reporter assay showed that these natural lariat RNAs had the potential to protect expression of miRNA target genes. Collectively, our results uncover that selective retention of lariat RNA is an actively regulatory process, and provide new insights into understanding how lariat RNA metabolism may impact miRNA activity.

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.

50 years of sperm cell isolations: from structural to omic studies

The fusion of male and female gametes is a fundamental process in the perpetuation and diversification of species. During the last 50 years, significant efforts have been made to isolate and characterize sperm cells from flowering plants, and to identify how these cells interact with female gametes to achieve double fertilization. The first techniques and analytical approaches not only provided structural and biochemical characterizations of plant sperm cells but also paved the way for in vitro fertilization studies. Further technological advances then led to unique insights into sperm biology at transcriptomic, proteomic and epigenetic level. Starting with a historical overview of sperm cell isolation techniques, we provide examples of how these contributed to create our current knowledge of sperm cell biology, and point out remaining challenges.

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.

Genome-Wide Identification of the Aconitase Gene Family in Tomato (Solanum lycopersicum) and CRISPR-Based Functional Characterization of SlACO2 on Male-Sterility

Tomato (Solanum lycopersicum) is one of the most cultivated vegetables in the world due to its consumption in a large variety of raw, cooked, or processed foods. Tomato breeding and productivity highly depend on the use of hybrid seeds and their higher yield, environmental adaption, and disease tolerance. However, the emasculation procedure during hybridization raises tomato seed production costs and labor expenses. Using male sterility is an effective way to reduce the cost of hybrid seeds and ensure cultivar purity. Recent developments in CRISPR genome editing technology enabled tomato breeders to investigate the male sterility genes and to develop male-sterile tomato lines. In the current study, the tomato Acotinase (SlACO) gene family was investigated via in silico tools and functionally characterized with CRISPR/Cas9-mediated gene disruption. Genome-wide blast and HMM search represented two SlACO genes located on different tomato chromosomes. Both genes were estimated to have a segmental duplication in the tomato genome due to their identical motif and domain structure. One of these genes, SlACO2, showed a high expression profile in all generative cells of tomato. Therefore, the SlACO2 gene was targeted with two different gRNA/Cas9 constructs to identify their functional role in tomatoes. The gene was mutated in a total of six genome-edited tomato lines, two of which were homozygous. Surprisingly, pollen viability was found to be extremely low in mutant plants compared to their wild-type (WT) counterparts. Likewise, the number of seeds per fruit also sharply decreased more than fivefold in mutant lines (10-12 seeds) compared to that in WT (67 seeds). The pollen shape, anther structures, and flower colors/shapes were not significantly varied between the mutant and WT tomatoes. The mutated lines were also subjected to salt and mannitol-mediated drought stress to test the effect of SlACO2 on abiotic stress tolerance. The results of the study indicated that mutant tomatoes have higher tolerance with significantly lower MDA content under stress conditions. This is the first CRISPR-mediated characterization of ACO genes on pollen viability, seed formation, and abiotic stress tolerance in tomatoes.

From birth to function: Male gametophyte development in flowering plants

Male germline development in flowering plants involves two distinct and successive phases, microsporogenesis and microgametogenesis, which involve one meiosis followed by two rounds of mitosis. Many aspects of distinctions after mitosis between the vegetative cell and the male germ cells are seen, from morphology to structure, and the differential functions of the two cell types in the male gametophyte are differentially needed and required for double fertilization. The two sperm cells, carriers of the hereditary substances, depend on the vegetative cell/pollen tube to be delivered to the female gametophyte for double fertilization. Thus, the intercellular communication and coordinated activity within the male gametophyte probably represent the most subtle regulation in flowering plants to guarantee the success of reproduction. This review will focus on what we have known about the differentiation process and the functional diversification of the vegetative cell and the male germ cell, the most crucial cell types for plant fertility and crop production.

The rice VCS1 is identified as a molecular tool to mark and visualize the vegetative cell of pollen

Cell-type-specific markers are valuable tools to reveal developmental processes underlying cell differentiation during plant reproduction. Here we report the pollen vegetative cell marker gene VCS1 (Vegetative Cell Specific 1) of rice (Oryza sativa japonica). VCS1 was expressed specifically in late pollen and was predicted to encode a small FAF domain-containing protein of 205 amino acid residues (aa). The expression of reporter fusion proteins showed that VCS1 was exclusively targeted to the vegetative nucleus of pollen. Upon pollen germination, VCS1 lost vegetative nucleus localization, and appeared diffused in the vegetative cytoplasm of pollen grain but not in the pollen tube. T-DNA insertional mutation which disrupted the carboxyl-terminal 21 aa of VCS1 did not affect plant vegetative growth and pollen development, while destruction of VCS1 by CRISPR/Cas9 only moderately affect pollen viability. VCS1 is evolutionally conserved in monocots but appeared absent in dicotyledons. This study reveals a molecular tool for visualizing the vegetative cell of rice and possible other monocots, which has potential values in the genetic engineering of male-sterile lines.

Comparative transcriptomic analysis reveals conserved programmes underpinning organogenesis and reproduction in land plants

The appearance of plant organs mediated the explosive radiation of land plants, which shaped the biosphere and allowed the establishment of terrestrial animal life. The evolution of organs and immobile gametes required the coordinated acquisition of novel gene functions, the co-option of existing genes and the development of novel regulatory programmes. However, no large-scale analyses of genomic and transcriptomic data have been performed for land plants. To remedy this, we generated gene expression atlases for various organs and gametes of ten plant species comprising bryophytes, vascular plants, gymnosperms and flowering plants. A comparative analysis of the atlases identified hundreds of organ- and gamete-specific orthogroups and revealed that most of the specific transcriptomes are significantly conserved. Interestingly, our results suggest that co-option of existing genes is the main mechanism for evolving new organs. In contrast to female gametes, male gametes showed a high number and conservation of specific genes, which indicates that male reproduction is highly specialized. The expression atlas capturing pollen development revealed numerous transcription factors and kinases essential for pollen biogenesis and function.

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.

Male gametophyte development in flowering plants: A story of quarantine and sacrifice

Over 160 years ago, scientists made the first microscopic observations of angiosperm pollen. Unlike in animals, male meiosis in angiosperms produces a haploid microspore that undergoes one asymmetric division to form a vegetative cell and a generative cell. These two cells have distinct fates: the vegetative cell exits the cell cycle and elongates to form a tip-growing pollen tube; the generative cell divides once more in the pollen grain or within the growing pollen tube to form a pair of sperm cells. The concept that male germ cells are less active than the vegetative cell came from biochemical analyses and pollen structure anatomy early in the last century and is supported by the pollen transcriptome data of the last decade. However, the mechanism of how and when the transcriptional repression in male germ cells occurs is still not fully understood. In this review, we provide a brief account of the cytological and metabolic differentiation between the vegetative cell and male germ cells, with emphasis on the role of temporary callose walls, dynamic nuclear pore density, transcription repression, and histone variants. We further discuss the intercellular movement of small interfering RNA (siRNA) derived from transposable elements (TEs) and reexamine the function of TE expression in male germ cells.

DNA methylation dynamics of sperm cell lineage development in tomato

During the sexual reproduction of higher plants, DNA methylation and transcription are broadly changed to reshape a microspore into two sperm cells (SCs) and a vegetative cell (VC). However, when and how the DNA methylation of SCs is established remains not fully understood. Here we investigate the DNA methylation (5 mC) dynamics of SC lineage and the VC in tomato using whole-genome bisulfite sequencing. We find the asymmetric division of the microspore gives its two daughter cells differential methylome. Compared with the generative cell (GC), the VC is hypomethylated at CG sites while hypermethylated at CHG and CHH sites, with the majority of differentially methylation regions targeted to transposable elements (TEs). SCs have a nearly identical DNA methylome to the GC, suggesting that the methylation landscape in SCs may be pre-established following the asymmetric division or inherited from the GC. The random forest classifier for predicting gene and TE expression shows that methylation within the gene body is a more powerful predictor for gene expression. Among all tested samples, gene and TE expression in the microspore may be more predictable by DNA methylation. Our results depict an intact DNA methylome landscape of SC lineage in higher plants, and reveal that the impact of DNA methylation on transcription is variant in different cell types.

Transcriptomics of Arabidopsis sperm cells at single-cell resolution

We present a detailed protocol for isolation of single sperm cells and transcriptome analysis to study variation in gene expression between sperm cells. Male gametophyte development in flowering plants begins with a microspore mother cell, which upon two consecutive cell divisions forms a mature pollen grain containing a vegetative nucleus and two sperm cells. Pollen development is a highly dynamic process, involving changes at both the transcriptome and epigenome levels of vegetative nuclei and the pair of sperm cells that have their own cytoplasm and nucleus. While the overall transcriptome of Arabidopsis pollen development is well documented, studies at single-cell level, in particular of sperm cells, are still lacking. Such studies would be essential to understand whether and how the two sperm cells are transcriptionally different, in particular once the pollen tube grows through the transmitting tissue of the pistil. Here we describe a detailed protocol for isolation of single sperm cells from growing pollen tubes and analysis of their transcriptome.