A study of the heterochronic sense/antisense RNA representation in florets of sexual and apomictic Paspalum notatum.

BACKGROUND: Apomixis, an asexual mode of plant reproduction, is a genetically heritable trait evolutionarily related to sexuality, which enables the fixation of heterozygous genetic combinations through the development of maternal seeds. Recently, reference floral transcriptomes were generated from sexual and apomictic biotypes of Paspalum notatum, one of the most well-known plant models for the study of apomixis. However, the transcriptome dynamics, the occurrence of apomixis vs. sexual expression heterochronicity across consecutive developmental steps and the orientation of transcription (sense/antisense) remain unexplored. RESULTS: We produced 24 Illumina TruSeq®/ Hiseq 1500 sense/antisense floral transcriptome libraries covering four developmental stages (premeiosis, meiosis, postmeiosis, and anthesis) in biological triplicates, from an obligate apomictic and a full sexual genotype. De novo assemblies with Trinity yielded 103,699 and 100,114 transcripts for the apomictic and sexual samples respectively. A global comparative analysis involving reads from all developmental stages revealed 19,352 differentially expressed sense transcripts, of which 13,205 (68%) and 6147 (32%) were up- and down-regulated in apomictic samples with respect to the sexual ones. Interestingly, 100 differentially expressed antisense transcripts were detected, 55 (55%) of them up- and 45 (45%) down-regulated in apomictic libraries. A stage-by-stage comparative analysis showed a higher number of differentially expressed candidates due to heterochronicity discrimination: the highest number of differential sense transcripts was detected at premeiosis (23,651), followed by meiosis (22,830), postmeiosis (19,100), and anthesis (17,962), while the highest number of differential antisense transcripts were detected at anthesis (495), followed by postmeiosis (164), meiosis (120) and premeiosis (115). Members of the AP2, ARF, MYB and WRKY transcription factor families, as well as the auxin, jasmonate and cytokinin plant hormone families appeared broadly deregulated. Moreover, the chronological expression profile of several well-characterized apomixis controllers was examined in detail. CONCLUSIONS: This work provides a quantitative sense/antisense gene expression catalogue covering several subsequent reproductive developmental stages from premeiosis to anthesis for apomictic and sexual P. notatum, with potential to reveal heterochronic expression between reproductive types and discover sense/antisense mediated regulation. We detected a contrasting transcriptional and hormonal control in apomixis and sexuality as well as specific sense/antisense modulation occurring at the onset of parthenogenesis.

Small RNA-seq reveals novel regulatory components for apomixis in Paspalum notatum.

BACKGROUND: Apomixis is considered an evolutionary deviation of the sexual reproductive pathway leading to the generation of clonal maternal progenies by seeds. Recent evidence from model and non-model species suggested that this trait could be modulated by epigenetic mechanisms involving small RNAs (sRNAs). Here we profiled floral sRNAs originated from apomictic and sexual Paspalum notatum genotypes in order to identify molecular pathways under epigenetic control that might be involved in the transition from sexuality to agamospermy. RESULTS: The mining of genes participating in sRNA-directed pathways from floral Paspalum transcriptomic resources showed these routes are functional during reproductive development, with several members differentially expressed in apomictic and sexual plants. Triplicate floral sRNA libraries derived from apomictic and a sexual genotypes were characterized by using high-throughput sequencing technology. EdgeR was apply to compare the number of sRNA reads between sexual and apomictic libraries that map over all Paspalum floral transcripts. A total of 1525 transcripts showed differential sRNA representation, including genes related to meiosis, plant hormone signaling, biomolecules transport, transcription control and cell cycle. Survey for miRNA precursors on transcriptome and genome references allowed the discovery of 124 entities, including 40 conserved and 8 novel ones. Fifty-six clusters were differentially represented in apomictic and sexual plants. All differentially expressed miRNAs were up-regulated in apomictic libraries but miR2275, which showed different family members with opposed representation. Examination of predicted miRNAs targets detected 374 potential candidates. Considering sRNA, miRNAs and target surveys together, 14 genes previously described as related with auxin metabolism, transport and signaling were detected, including AMINO ACID/AUXIN PERMEASE 15, IAA-AMIDO SYNTHETASE GH3-8, IAA30, miR160, miR167, miR164, miR319, ARF2, ARF8, ARF10, ARF12, AFB2, PROLIFERATING CELL FACTOR 6 and NITRATE TRANSPORTER 1.1. CONCLUSIONS: This work provides a comprehensive survey of the sRNA differential representation in flowers of sexual and apomictic Paspalum notatum plants. An integration of the small RNA profiling data presented here and previous transcriptomic information suggests that sRNA-mediated regulation of auxin pathways is pivotal in promoting apomixis. These results will underlie future functional characterization of the molecular components mediating the switch from sexuality to apomixis.

A reference floral transcriptome of sexual and apomictic Paspalum notatum.

BACKGROUND: Paspalum notatum Flügge is a subtropical grass native to South America, which includes sexual diploid and apomictic polyploid biotypes. In the past decade, a number of apomixis-associated genes were discovered in this species through genetic mapping and differential expression surveys. However, the scarce information on Paspalum sequences available in public databanks limited annotations and functional predictions for these candidates. RESULTS: We used a long-read 454/Roche FLX+ sequencing strategy to produce robust reference transcriptome datasets from florets of sexual and apomictic Paspalum notatum genotypes and delivered a list of transcripts showing differential representation in both reproductive types. Raw data originated from floral samples collected from premeiosis to anthesis was assembled in three libraries: i) sexual (SEX), ii) apomictic (APO) and iii) global (SEX + APO). A group of physically-supported Paspalum mRNA and EST sequences matched with high level of confidence to both sexual and apomictic libraries. A preliminary trial allowed discovery of the whole set of putative alleles/paralogs corresponding to 23 previously identified apomixis-associated candidate genes. Moreover, a list of 3,732 transcripts and several co-expression and protein -protein interaction networks associated with apomixis were identified. CONCLUSIONS: The use of the 454/Roche FLX+ transcriptome database will allow the detailed characterization of floral alleles/paralogs of apomixis candidate genes identified in prior and future work. Moreover, it was used to reveal additional candidate genes differentially represented in apomictic and sexual flowers. Gene ontology (GO) analyses of this set of transcripts indicated that the main molecular pathways altered in the apomictic genotype correspond to specific biological processes, like biotic and abiotic stress responses, growth, development, cell death and senescence. This data collection will be of interest to the plant reproduction research community and, particularly, to Paspalum breeding projects.

First insight into divergence, representation and chromosome distribution of reverse transcriptase fragments from L1 retrotransposons in peanut and wild relative species.

Peanut is an allotetraploid (2n = 2x = 40, AABB) of recent origin. Arachis duranensis and A. ipaënsis, the most probable diploid ancestors of the cultigen, and several other wild diploid species with different genomes (A, B, D, F and K) are used in peanut breeding programs. However, the genomic relationships and the evolutionary pathways of genome differentiation of these species are poorly understood. We performed a sequence-based phylogenetic analysis of the L1 reverse transcriptase and estimated its representation and chromosome distribution in species of five genomes and three karyotype groups with the aim of contributing to the knowledge of the genomic structure and evolution of peanut and wild diploid relatives. All the isolated rt fragments were found to belong to plant L1 lineage and were named ALI. The best supported phylogenetic groups were not concordant with the genomes or karyotype groups. The copy number of ALI sequences was higher than the expected one for plants and directly related to genome size. FISH experiments revealed that ALI is mainly located on the euchromatin of interstitial and distal regions of most chromosome arms. Divergence of ALI sequences would have occurred before the differentiation of the genomes and karyotype groups of Arachis. The representation and chromosome distribution of ALI in peanut was almost additive of those of the parental species suggesting that the spontaneous hybridization of the two parental species of peanut followed by chromosome doubling would not have induced a significant burst of ALI transposition.

Characterization and expression analysis of SOMATIC EMBRYOGENESIS RECEPTOR KINASE (SERK) genes in sexual and apomictic Paspalum notatum.

The SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE (SERK) gene plays a fundamental role in somatic embryogenesis of angiosperms, and is associated with apomixis in Poa pratensis. The objective of this work was to isolate, characterize and analyze the expression patterns of SERK genes in apomictic and sexual genotypes of Paspalum notatum. A conserved 200-bp gene fragment was amplified from genomic DNA with heterologous primers, and used to initiate a chromosomal walking strategy for cloning the complete sequence. This procedure allowed the isolation of two members of the P. notatum SERK family; PnSERK1, which is similar to PpSERK1, and PnSERK2, which is similar to ZmSERK2 and AtSERK1. Phylogenetic analyses indicated that PnSERK1 and PnSERK2 represent paralogous sequences. Southern-blot hybridization indicated the presence of at least three copies of SERK genes in the species. qRT-PCR analyses revealed that PnSERK2 was expressed at significantly higher levels than PnSERK1 in roots, leaves, reproductive tissues and embryogenic calli. Moreover, in situ hybridization experiments revealed that PnSERK2 displayed a spatially and chronologically altered expression pattern in reproductive organs of the apomictic genotype with respect to the sexual one. PnSERK2 is expressed in nucellar cells of the apomictic genotype at meiosis, but only in the megaspore mother cell in the sexual genotype. Therefore, apomixis onset in P. notatum seems to be correlated with the expression of PnSERK2 in nucellar tissue.

A methylation status analysis of the apomixis-specific region in Paspalum spp. suggests an epigenetic control of parthenogenesis.

Apomixis, a clonal plant reproduction by seeds, is controlled in Paspalum spp. by a single locus which is blocked in terms of recombination. Partial sequence analysis of the apomixis locus revealed structural features of heterochromatin, namely the presence of repetitive elements, gene degeneration, and de-regulation. To test the epigenetic control of apomixis, a study on the distribution of cytosine methylation at the apomixis locus and the effect of artificial DNA demethylation on the mode of reproduction was undertaken in two apomictic Paspalum species. The 5-methylcytosine distribution in the apomixis-controlling genomic region was studied in P. simplex by methylation-sensitive restriction fragment length polymorphism (RFLP) analysis and in P. notatum by fluorescene in situ hybridization (FISH). The effect of DNA demethylation was studied on the mode of reproduction of P. simplex by progeny test analysis of apomictic plants treated with the demethylating agent 5'-azacytidine. A high level of cytosine methylation was detected at the apomixis-controlling genomic region in both species. By analysing a total of 374 open pollination progeny, it was found that artificial demethylation had little or no effect on apospory, whereas it induced a significant depression of parthenogenesis. The results suggested that factors controlling repression of parthenogenesis might be inactivated in apomictic Paspalum by DNA methylation.

Sequence characterization, in silico mapping and cytosine methylation analysis of markers linked to apospory in Paspalum notatum.

In previous studies we reported the identification of several AFLP, RAPD and RFLP molecular markers linked to apospory in Paspalum notatum. The objective of this work was to sequence these markers, obtain their flanking regions by chromosome walking and perform an in silico mapping analysis in rice and maize. The methylation status of two apospory-related sequences was also assessed using methylation-sensitive RFLP experiments. Fourteen molecular markers were analyzed and several protein-coding sequences were identified. Copy number estimates and RFLP linkage analysis showed that the sequence PnMAI3 displayed 2-4 copies per genome and linkage to apospory. Extension of this marker by chromosome walking revealed an additional protein-coding sequence mapping in silico in the apospory-syntenic regions of rice and maize. Approximately 5 kb corresponding to different markers were characterized through the global sequencing procedure. A more refined analysis based on sequence information indicated synteny with segments of chromosomes 2 and 12 of rice and chromosomes 3 and 5 of maize. Two loci associated with apomixis locus were tested in methylation-sensitive RFLP experiments using genomic DNA extracted from leaves. Although both target sequences were methylated no methylation polymorphisms associated with the mode of reproduction were detected.

The Auxin-Response Repressor IAA30 Is Down-Regulated in Reproductive Tissues of Apomictic Paspalum notatum.

The capacity for apomixis in Paspalum notatum is controlled by a single-dominant genomic region, which shows strong synteny to a portion of rice chromosome 12 long arm. The locus LOC_Os12g40890, encoding the Auxin/Indole-3-Acetic Acid (Aux/IAA) family member OsIAA30, is located in this rice genomic segment. The objectives of this work were to identify transcripts coding for Aux/IAA proteins expressed in reproductive tissues of P. notatum, detect the OsIAA30 putative ortholog and analyze its temporal and spatial expression pattern in reproductive organs of sexual and apomictic plants. Thirty-three transcripts coding for AUX/IAA proteins were identified. Predicted protein alignment and phylogenetic analysis detected a highly similar sequence to OsIAA30 (named as PnIAA30) present in both sexual and apomictic samples. The expression assays of PnIAA30 showed a significant down-regulation in apomictic spikelets compared to sexual ones at the stages of anthesis and post-anthesis, representation levels negatively correlated with apospory expressivity and different localizations in sexual and apomictic ovules. Several PnIAA30 predicted interactors also appeared differentially regulated in the sexual and apomictic floral transcriptomes. Our results showed that an auxin-response repressor similar to OsIAA30 is down-regulated in apomictic spikelets of P. notatum and suggests a contrasting regulation of auxin signaling during sexual and asexual seed formation.

Spotting the Targets of the Apospory Controller TGS1 in Paspalum notatum.

Sexuality and apomixis are interconnected plant reproductive routes possibly behaving as polyphenic traits under the influence of the environment. In the subtropical grass Paspalum notatum, one of the controllers of apospory, a main component of gametophytic apomixis reproduction, is TRIMETHYLGUANOSINE SYNTHASE 1 (TGS1), a multifunctional gene previously associated with RNA cleavage regulation (including mRNA splicing as well as rRNA and miRNA processing), transcriptional modulation and the establishment of heterochromatin. In particular, the downregulation of TGS1 induces a sexuality decline and the emergence of aposporous-like embryo sacs. The present work was aimed at identifying TGS1 target RNAs expressed during reproductive development of Paspalum notatum. First, we mined available RNA databases originated from spikelets of sexual and apomictic plants, which naturally display a contrasting TGS1 representation, to identify differentially expressed mRNA splice variants and miRNAs. Then, the role of TGS1 in the generation of these particular molecules was investigated in antisense tgs1 sexual lines. We found that CHLOROPHYLL A-B BINDING PROTEIN 1B-21 (LHC Ib-21, a component of the chloroplast light harvesting complex), QUI-GON JINN (QGJ, encoding a MAP3K previously associated with apomixis) and miR2275 (a meiotic 24-nt phasi-RNAs producer) are directly or indirectly targeted by TGS1. Our results point to a coordinated control exercised by signal transduction and siRNA machineries to induce the transition from sexuality to apomixis.

Expression of lorelei-like genes in aposporous and sexual Paspalum notatum plants.

Gametophytic apomictic plants form non-reduced embryo sacs that generate clonal embryos by parthenogenesis, in the absence of both meiosis and egg-cell fertilization. Here we report the sequence and expression analysis of a lorelei-like Paspalum notatum gene, n20gap-1, which encodes a GPI-anchored protein previously associated with apomixis in this species. Phylogeny trees showed that n20gap-1 was evolutionary related to the Arabidopsis thaliana lorelei genes At4g26466 and At5g56170. The lorelei At4g26466 disruption was shown to be detrimental to sperm cell release in arabidopsis. RFLP (Restriction Fragment Length Polymorphism) analysis revealed the occurrence of several homologous sequences in the Paspalum notatum genome, exhibiting polymorphisms genetically linked to apomixis. Real-time PCR showed that lorelei-family genes present a minor activity peak at pre-meiosis and a major one at anthesis. The apomictic genotype analyzed showed a significantly increased activity at pre-meiosis, post-meiosis and anthesis with respect to a sexual genotype. In situ hybridization assays revealed expression in integuments, nucellus and the egg-cell apparatus. Several n20gap-1 alleles differing mainly at the 3' UTR sequence were identified. Allele-specific real-time PCR experiments showed that allele 28 was significantly induced in reproductive tissues of the apomictic genotype with respect to the sexual genotype at anthesis. Our results indicate that P. notatum lorelei-like genes are differentially expressed in representative sexual (Q4188) and apomictic (Q4117) genotypes, and might play a role in the final stages of the apomixis developmental cascade. However, the association of n20gap-1 expression with the trait should be confirmed in significant number of sexual and apomictic genotypes.

Characterization of retrotransposon sequences expressed in inflorescences of apomictic and sexual Paspalum notatum plants.

Apomixis, an asexual mode of reproduction through seeds, holds much promise for agricultural advances. However, the molecular mechanisms underlying this trait are still poorly understood. We previously isolated several transcripts representing novel sequences differentially expressed in reproductive tissues of sexual and apomictic plants. Here, we report the characterization of two of these unknown RNA transcripts (experimental codes N17 and N22). Since original fragments showed no significant homologies to sequences at databases, preliminary genomic PCR experiments were carried out to discard possible contaminations. RACE extension on flanking regions provided longer sequences for the candidates and additional related transcripts, which revealed similarity to LTR retrotransposons carrying short transduplicated segments of protein-coding genes. Interestingly, some transduplicated segments corresponded to genes previously associated with apomictic development. Gene copy number estimations revealed a moderate representation of the elements in the genome, with significantly increased numbers in a sexual genotype with respect to an apomictic one. Genetic mapping of N17 showed that a copy of this particular element was located onto Paspalum notatum linkage group F3c, at a central non-recombinant region resembling a centromere. Expression analysis showed an increased activity of N17 and N22 sense strands in ovules of the sexual genotypes. A retrotransposon-specific differential display analysis aimed at detecting related sequences allowed the identification of a complex family, with the majority of its members represented in the sexual genotype. Our results suggest that these elements could be participating in regulatory pathways related to apomixis and sexuality.

Single-cell transcriptome profiling of buffelgrass (Cenchrus ciliaris) eggs unveils apomictic parthenogenesis signatures.

Apomixis, a type of asexual reproduction in angiosperms, results in progenies that are genetically identical to the mother plant. It is a highly desirable trait in agriculture due to its potential to preserve heterosis of F1 hybrids through subsequent generations. However, no major crops are apomictic. Deciphering mechanisms underlying apomixis becomes one of the alternatives to engineer self-reproducing capability into major crops. Parthenogenesis, a major component of apomixis, commonly described as the ability to initiate embryo formation from the egg cell without fertilization, also can be valuable in plant breeding for doubled haploid production. A deeper understanding of transcriptional differences between parthenogenetic and sexual or non-parthenogenetic eggs can assist with pathway engineering. By conducting laser capture microdissection-based RNA-seq on sexual and parthenogenetic egg cells on the day of anthesis, a de novo transcriptome for the Cenchrus ciliaris egg cells was created, transcriptional profiles that distinguish the parthenogenetic egg from its sexual counterpart were identified, and functional roles for a few transcription factors in promoting natural parthenogenesis were suggested. These transcriptome data expand upon previous gene expression studies and will be a resource for future research on the transcriptome of egg cells in parthenogenetic and sexual genotypes.

Differential Epigenetic Marks Are Associated with Apospory Expressivity in Diploid Hybrids of Paspalum rufum.

Apomixis seems to emerge from the deregulation of preexisting genes involved in sexuality by genetic and/or epigenetic mechanisms. The trait is associated with polyploidy, but diploid individuals of Paspalum rufum can form aposporous embryo sacs and develop clonal seeds. Moreover, diploid hybrid families presented a wide apospory expressivity variation. To locate methylation changes associated with apomixis expressivity, we compare relative DNA methylation levels, at CG, CHG, and CHH contexts, between full-sib P. rufum diploid genotypes presenting differential apospory expressivity. The survey was performed using a methylation content-sensitive enzyme ddRAD (MCSeEd) strategy on samples at premeiosis/meiosis and postmeiosis stages. Based on the relative methylation level, principal component analysis and heatmaps, clearly discriminate samples with contrasting apospory expressivity. Differential methylated contigs (DMCs) showed 14% of homology to known transcripts of Paspalum notatum reproductive transcriptome, and almost half of them were also differentially expressed between apomictic and sexual samples. DMCs showed homologies to genes involved in flower growth, development, and apomixis. Moreover, a high proportion of DMCs aligned on genomic regions associated with apomixis in Setaria italica. Several stage-specific differential methylated sequences were identified as associated with apospory expressivity, which could guide future functional gene characterization in relation to apomixis success at diploid and tetraploid levels.

The MAP3K-Coding QUI-GON JINN (QGJ) Gene Is Essential to the Formation of Unreduced Embryo Sacs in Paspalum.

Apomixis is a clonal mode of reproduction via seeds, which results from the failure of meiosis and fertilization in the sexual female reproductive pathway. In previous transcriptomic surveys, we identified a mitogen-activated protein kinase kinase kinase (N46) displaying differential representation in florets of sexual and apomictic Paspalum notatum genotypes. Here, we retrieved and characterized the N46 full cDNA sequence from sexual and apomictic floral transcriptomes. Phylogenetic analyses showed that N46 was a member of the YODA family, which was re-named QUI-GON JINN (QGJ). Differential expression in florets of sexual and apomictic plants was confirmed by qPCR. In situ hybridization experiments revealed expression in the nucellus of aposporous plants' ovules, which was absent in sexual plants. RNAi inhibition of QGJ expression in two apomictic genotypes resulted in significantly reduced rates of aposporous embryo sac formation, with respect to the level detected in wild type aposporous plants and transformation controls. The QGJ locus segregated independently of apospory. However, a probe derived from a related long non-coding RNA sequence (PN_LNC_QGJ) revealed RFLP bands cosegregating with the Paspalum apospory-controlling region (ACR). PN_LNC_QGJ is expressed in florets of apomictic plants only. Our results indicate that the activity of QGJ in the nucellus of apomictic plants is necessary to form non-reduced embryo sacs and that a long non-coding sequence with regulatory potential is similar to sequences located within the ACR.

Haploid embryo production in rice and maize induced by PsASGR-BBML transgenes.

KEY MESSAGE: The PsASGR - BBML transgene, derived from a wild apomictic grass species, can induce parthenogenesis, embryo formation without fertilization, in rice and maize, leading to the formation of haploid plants. The ability to engineer apomictic crop plants using genes identified from naturally occurring apomicts will depend on the ability of those genes to function in crop plants. The PsASGR-BBML transgene, derived from the apomictic species Pennisetum squamulatum, promotes parthenogenesis in sexual pearl millet, a member of the same genus, leading to the formation of haploid embryos. This study determined that the PsASGR-BBML transgene can induce haploid embryo development in two major monocot crops, maize and rice. Transgene variations tested included two different promoters and the use of both genomic and cDNA PsASGR-BBML-derived sequences. Haploid plants were recovered from mature caryopses (seed) of rice and maize lines at variable rates. The PsASGR-BBML transgenes failed to induce measurable haploid seed development in the model genetic plant system Arabidopsis thaliana. Complexity of embryo development, as documented in transgenic rice lines, identifies the need for further characterization of the PsASGR-BBML gene.

Gene expression analysis at the onset of aposporous apomixis in Paspalum notatum.

Apomixis is a route of asexual reproduction through seeds, that progresses in the absence of meiosis and fertilization to generate maternal clonal progenies. Gametophytic apomicts are usually polyploid and probably arose from sexual ancestors through a limited number of mutations in the female reproductive pathway. A differential display analysis was carried out on immature inflorescences of sexual and apomictic tetraploid genotypes of Paspalum notatum, in order to identify genes associated with the emergence of apospory. Analysis of approximately 10,000 transcripts led to the identification of 94 high-quality differentially expressed sequences. Assembling analysis, plus validation, rendered 65 candidate unigenes, organized as 14 contigs and 51 singletons. Thirty-four unigenes were isolated from apomictic plants and 31 from sexual ones. A total of 45 (69.2%) unigenes were functionally categorized. While several of the differentially expressed sequences appeared to be components of an extracellular receptor kinase (ERK) signal transduction cascade, others seemed to participate in a variety of central cellular processes like cell-cycle control, protein turnover, intercellular signalling, transposon activity, transcriptional regulation and endoplasmic reticulum-mediated biosynthesis. In silico mapping revealed that a particular group of five genes silenced in apomictic plants clustered in a rice genomic area syntenic with the region governing apospory in Paspalum notatum and Brachiaria brizantha. Two of these genes mapped within the set of apo-homologues in P. notatum. Four genes previously reported to be controlled by ploidy were identified among those expressed differentially between apomictic and sexual plants. In situ hybridization experiments were performed for selected clones.