Rice expression atlas in reproductive development

Transcriptome-based prediction for polygenic traits in rice using different gene subsets

BACKGROUND

Transcriptome-based prediction of complex phenotypes is a relatively new statistical method that links genetic variation to phenotypic variation. The selection of large-effect genes based on a priori biological knowledge is beneficial for predicting oligogenic traits; however, such a simple gene selection method is not applicable to polygenic traits because causal genes or large-effect loci are often unknown. Here, we used several gene-level features and tested whether it was possible to select a gene subset that resulted in better predictive ability than using all genes for predicting a polygenic trait.

RESULTS

Using the phenotypic values of shoot and root traits and transcript abundances in leaves and roots of 57 rice accessions, we evaluated the predictive abilities of the transcriptome-based prediction models. Leaf transcripts predicted shoot phenotypes, such as plant height, more accurately than root transcripts, whereas root transcripts predicted root phenotypes, such as crown root length, more accurately than leaf transcripts. Furthermore, we used the following three features to train the prediction model: (1) tissue specificity of the transcripts, (2) ontology annotations, and (3) co-expression modules for selecting gene subsets. Although models trained by a gene subset often resulted in lower predictive abilities than the model trained by all genes, some gene subsets showed improved predictive ability. For example, using genes expressed in roots but not in leaves, the predictive ability for crown root diameter was improved by more than 10% (R2 = 0.59 when using all genes; R2 = 0.66, using 1,554 root-specifically expressed genes). Similarly, genes annotated as "gibberellic acid sensitivity" showed higher predictive ability than using all genes for root dry weight.

CONCLUSIONS

Our results highlight both the possibility and difficulty of selecting an appropriate gene subset to predict polygenic traits from transcript abundance, given the current biological knowledge and information. Further integration of multiple sources of information, as well as improvements in gene characterization, may enable the selection of an optimal gene set for the prediction of polygenic phenotypes.

m6A demethylase OsALKBH5 is required for double-strand break formation and repair by affecting mRNA stability in rice meiosis

N6-methyladenosine (m6A) RNA modification is the most prevalent messenger RNA (mRNA) modification in eukaryotes and plays critical roles in the regulation of gene expression. m6A is a reversible RNA modification that is deposited by methyltransferases (writers) and removed by demethylases (erasers). The function of m6A erasers in plants is highly diversified and their roles in cereal crops, especially in reproductive development essential for crop yield, are largely unknown. Here, we demonstrate that rice OsALKBH5 acts as an m6A demethylase required for the normal progression of male meiosis. OsALKBH5 is a nucleo-cytoplasmic protein, highly enriched in rice anthers during meiosis, that associates with P-bodies and exon junction complexes, suggesting that it is involved in regulating mRNA processing and abundance. Mutations of OsALKBH5 cause reduced double-strand break (DSB) formation, severe defects in DSB repair, and delayed meiotic progression, leading to complete male sterility. Transcriptome analysis and m6A profiling indicate that OsALKBH5-mediated m6A demethylation stabilizes the mRNA level of multiple meiotic genes directly or indirectly, including several genes that regulate DSB formation and repair. Our study reveals the indispensable role of m6A metabolism in post-transcriptional regulation of meiotic progression in rice.

Cytological observation and RNA-seq analysis reveal novel miRNAs high expression associated with the pollen fertility of neo-tetraploid rice

BACKGROUND

Neo-tetraploid rice lines exhibit high fertility and strong heterosis and harbor novel specific alleles, which are useful germplasm for polyploid rice breeding. However, the mechanism of the fertility associated with miRNAs remains unknown. In this study, a neo-tetraploid rice line, termed Huaduo21 (H21), was used. Cytological observation and RNA-sequencing were employed to identify the fertility-related miRNAs in neo-tetraploid rice.

RESULTS

H21 showed high pollen fertility (88.08%), a lower percentage of the pollen mother cell (PMC) abnormalities, and lower abnormalities during double fertilization and embryogenesis compared with autotetraploid rice. A total of 166 non-additive miRNAs and 3108 non-additive genes were detected between H21 and its parents. GO and KEGG analysis of non-additive genes revealed significant enrichments in the DNA replication, Chromosome and associated proteins, and Replication and repair pathways. Comprehensive multi-omics analysis identified 32 pairs of miRNA/target that were associated with the fertility in H21. Of these, osa-miR408-3p and osa-miR528-5p displayed high expression patterns, targeted the phytocyanin genes, and were associated with high pollen fertility. Suppression of osa-miR528-5p in Huaduo1 resulted in a low seed set and a decrease in the number of grains. Moreover, transgenic analysis implied that osa-MIR397b-p3, osa-miR5492, and osa-MIR5495-p5 might participate in the fertility of H21.

CONCLUSION

Taken together, the regulation network of fertility-related miRNAs-targets pairs might contribute to the high seed setting in neo-tetraploid rice. These findings enhance our understanding of the regulatory mechanisms of pollen fertility associated with miRNAs in neo-tetraploid rice.

Field and genetic evidence support the photosynthetic performance index (PIABS) as an indicator of rice grain yield

The rice breeding process for grain yield could be effectively enhanced by developing efficient tools that accelerate plant selection through the rapid determination of reliable predictors. In this study, we have described various associations between grain yield and photosynthetic parameters, which can be easily and quickly obtained using a non-invasive technique on the flag leaf during the anthesis stage. Among the analyzed photosynthetic parameters, the photosynthetic performance index (PIABS) stood out due to its strong association with grain yield. A genome-wide association analysis conducted on plants from a rice diversity panel at the tillering stage revealed the presence of a quantitative trait locus on chromosome 9. This locus was characterized by a group of candidate chloroplastic genes that exhibited contrasting haplotypes for PIABS. An analysis of these haplotypes revealed a clear division into two groups. One group consisted of haplotypes linked to high values of PIABS, which were predominantly associated with Japonica spp. subpopulations. The other group consisted of haplotypes linked to low values of PIABS, which were exclusively associated with Indica spp. subpopulations. Japonica spp. genotypes exhibited higher values in the yield component panicle weight compared with the Indica spp. genotypes. The findings of this study indicate that PIABS could serve as an early predictor of yield parameters during the tillering stage in rice breeding processes.

Comparative cytological and transcriptome analyses of ny2 mutant delayed degeneration of tapetal cells and promotes abnormal microspore development in neo-tetraploid rice

We aimed to investigate the genetic defects related to pollen development and infertility in NY2, a novel tetraploid rice germplasm known as Neo-tetraploid rice. This rice variety was created through the crossbreeding and selective breeding of various autotetraploid rice lines and has previously shown high fertility. Our previous research has revealed that the NY2 gene, encoding a eukaryotic translation initiation factor 3 subunit E, regulates pollen fertility. However, the underlying mechanism behind this fertility is yet to be understood. To shed light on this matter, we performed a combined cytological and transcriptome analysis of the NY2 gene. Cytological analysis indicated that ny2 underwent abnormal tapetal cells, microspore, and middle layer development, which led to pollen abortion and ultimately to male sterility. Genetic analysis revealed that the F₁ plants showed normal fertility and an obvious advantage for seed setting compared to ny2. Global gene expression analysis in ny2 revealed a total of 7545 genes were detected at the meiosis stage, and 3925 and 3620 displayed upregulation and downregulation, respectively. The genes were significantly enriched for the gene ontology (GO) term "carbohydrate metabolic process. Moreover, 9 genes related to tapetum or pollen fertility showed down-regulation, such as OsABCG26 (ATP Binding Cassette G26), TMS9-1 (Thermosensitive Male Sterility), EAT1 (Programmed cell death regulatory), KIN14M (Kinesin Motor), OsMT1a (Metallothionein), and OsSTRL2 (Atypical strictosidine synthase), which were validated by qRT-PCR. Further analyses of DEGs identified nine down-regulated transcription factor genes related to pollen development. NY2 is an important regulator of the development of tapetum and microspore. The regulatory gene network described in this study may offer important understandings into the molecular processes that underlie fertility control in tetraploid rice.

Maintenance of Methyl-Esterified Pectin Level in Pollen Mother-Cell Stages Is Required for Microspore Development

Pectin modification and degradation are vital for plant development, although the underlying mechanisms are still not well understood. Furthermore, reports on the function of pectin in early pollen development are limited. We generated OsPME-FOX rice lines with little methyl-esterified pectin even in the early-pollen mother-cell stage due to overexpression of the gene encoding pectin-methylesterase. Overexpression of OsPME1 in rice increased the activity of PME, which decreased the degree of pectin methyl esterification in the cell wall. OsPME1-FOX grew normally and showed abnormal phenotypes in anther and pollen development, especially in terms of the pollen mother-cell stage. In addition, we examined modifications of cell-wall polysaccharides at the cellular level using antibodies against polysaccharides. Immunohistochemical staining using LM19 and LM20 showed that methyl-esterified pectin distribution and the pectin contents in pollen mother-cell wall decreased in OsPME1-FOX compared with the wild type. Thus, the maintenance of methyl-esterified pectin plays a role in degrading and maintaining the pollen mother-cell wall during microspore development.

Cytological, transcriptome and miRNome temporal landscapes decode enhancement of rice grain size

BACKGROUND

Rice grain size (GS) is an essential agronomic trait. Though several genes and miRNA modules influencing GS are known and seed development transcriptomes analyzed, a comprehensive compendium connecting all possible players is lacking. This study utilizes two contrasting GS indica rice genotypes (small-grained SN and large-grained LGR). Rice seed development involves five stages (S1-S5). Comparative transcriptome and miRNome atlases, substantiated with morphological and cytological studies, from S1-S5 stages and flag leaf have been analyzed to identify GS proponents.

RESULTS

Histology shows prolonged endosperm development and cell enlargement in LGR. Stand-alone and comparative RNAseq analyses manifest S3 (5-10 days after pollination) stage as crucial for GS enhancement, coherently with cell cycle, endoreduplication, and programmed cell death participating genes. Seed storage protein and carbohydrate accumulation, cytologically and by RNAseq, is shown to be delayed in LGR. Fourteen transcription factor families influence GS. Pathway genes for four phytohormones display opposite patterns of higher expression. A total of 186 genes generated from the transcriptome analyses are located within GS trait-related QTLs deciphered by a cross between SN and LGR. Fourteen miRNA families express specifically in SN or LGR seeds. Eight miRNA-target modules display contrasting expressions amongst SN and LGR, while 26 (SN) and 43 (LGR) modules are differentially expressed in all stages.

CONCLUSIONS

Integration of all analyses concludes in a "Domino effect" model for GS regulation highlighting chronology and fruition of each event. This study delineates the essence of GS regulation, providing scope for future exploits. The rice grain development database (RGDD) ( www.nipgr.ac.in/RGDD/index.php ; https://doi.org/10.5281/zenodo.7762870 ) has been developed for easy access of data generated in this paper.

An uncharacterized protein NY1 targets EAT1 to regulate anther tapetum development in polyploid rice

BACKGROUND

Autotetraploid rice is a useful germplasm for the breeding of polyploid rice; however, low fertility is a major hindrance for its utilization. Neo-tetraploid rice with high fertility was developed from the crossing of different autotetraploid rice lines. Our previous research showed that the mutant (ny1) of LOC_Os07g32406 (NY1), which was generated by CRISPR/Cas9 knock-out in neo-tetraploid rice, showed low pollen fertility, low seed set, and defective chromosome behavior during meiosis. However, the molecular genetic mechanism underlying the fertility remains largely unknown.

RESULTS

Here, cytological observations of the NY1 mutant (ny1) indicated that ny1 exhibited abnormal tapetum and middle layer development. RNA-seq analysis displayed a total of 5606 differentially expressed genes (DEGs) in ny1 compared to wild type (H1) during meiosis, of which 2977 were up-regulated and 2629 were down-regulated. Among the down-regulated genes, 16 important genes associated with tapetal development were detected, including EAT1, CYP703A3, CYP704B2, DPW, PTC1, OsABCG26, OsAGO2, SAW1, OsPKS1, OsPKS2, and OsTKPR1. The mutant of EAT1 was generated by CRISPR/Cas9 that showed abnormal tapetum and pollen wall formation, which was similar to ny1. Moreover, 478 meiosis-related genes displayed down-regulation at same stage, including 9 important meiosis-related genes, such as OsREC8, OsSHOC1, SMC1, SMC6a and DCM1, and their expression levels were validated by qRT-PCR.

CONCLUSIONS

Taken together, these results will aid in identifying the key genes associated with pollen fertility, which offered insights into the molecular mechanism underlying pollen development in tetraploid rice.

Interaction of ONION2 ketoacyl CoA synthase with ketoacyl CoA reductase of rice

BACKGROUND

Fatty acid elongases (FAEs), which catalyse elongation reactions of a carbon chain of very-long-chain fatty acids, play an important role in shoot development in rice. The elongation reactions consist of four sequential reactions catalysed by distinct enzymes, which are assumed to form an elongation complex. However, no interacting proteins of ONION1 (ONI1) and ONI2, which are ketoacyl CoA synthase catalyzing the first step and are required for shoot development in rice, are reported.

METHODS AND RESULTS

In this study ketoacyl CoA reductase (KCR) that interacts with ONI1 and ONI2 was searched. A database search identified 10 KCR genes in the rice genome. Among the genes, the expression pattern of KCR1 was similar to that of ONI2. Yeast two-hybrid analysis showed interaction of ONI2 with KCR1, which was confirmed by GST pull-down assay. No interacting partner of ONI1 was identified.

CONCLUSIONS

Our results suggest that ONI2 and KCR1 form an FAE complex that may play a role in biosynthesizing VLCFAs during shoot development.

Association mapping of drought tolerance and agronomic traits in rice (Oryza sativa L.) landraces

BACKGROUND

Asian cultivars were predominantly represented in global rice panel selected for sequencing and to identify novel alleles for drought tolerance. Diverse genetic resources adapted to Indian subcontinent were not represented much in spite harboring useful alleles that could improve agronomic traits, stress resilience and productivity. These rice accessions are valuable genetic resource in developing rice varieties suited to different rice ecosystem that experiences varying drought stress level, and at different crop stages. A core collection of rice germplasm adapted to Southwestern Indian peninsular genotyped using SSR markers and characterized by contrasting water regimes to associate genomic regions for physiological, root traits and yield related traits. Genotyping-By-Sequencing of selected accessions within the diverse panel revealed haplotype variation in genic content within genomic regions mapped for physiological, morphological and root traits.

RESULTS

Diverse rice panel (99 accessions) were evaluated in field and measurements on plant physiological, root traits and yield related traits were made over five different seasons experiencing varying drought stress intensity at different crop stages. Traits like chlorophyll stability index, leaf rolling, days to 50% flowering, chlorophyll content, root volume and root biomass were identified as best predictors of grain yield under stress. Association mapping revealed genetic variation among accessions and revealed 14 genomic targets associated with different physiological, root and plant production traits. Certain accessions were found to have beneficial allele to improve traits, plant height, root length and spikelet fertility, that contribute to the grain yield under stress. Genomic characterization of eleven accessions revealed haplotype variation within key genomic targets on chromosomes 1, 4, 6 and 11 for potential use as molecular markers to combine drought avoidance and tolerance traits. Genes mined within the genomic QTL intervals identified were prioritized based on tissue specific expression level in publicly available rice transcriptome data.

CONCLUSION

The genetic and genomic resources identified will enable combining traits with agronomic value to optimize yield under stress and hasten trait introgression into elite cultivars. Alleles associated with plant height, specific leaf area, root length from PTB8 and spikelet fertility and grain weight from PTB26 can be harnessed in future rice breeding program.

Genetic control of flowering in greater yam (Dioscorea alata L.)

BACKGROUND

Greater yam (Dioscorea alata L.) is a major tropical and subtropical staple crop cultivated for its starchy tubers. Breeding of this dioecious species is hampered by its erratic flowering, yet little is currently known on the genetic determinism of its sexual reproduction.

RESULT

Here we used a genome-wide association approach and identified a major genetic barrier to reproduction in yam on chromosome 1, as represented by two candidate genes. A deleterious effect on male fitness could be hypothesized considering the involvement of these two genes in male reproduction and the low frequency of this non-flowering dominant allele within the male genepool. We also extended the hypothesis of a XX/XY sex-determination system located on chromosome 6 in D. alata to encompass most of the species diversity. Moreover, a kompetitive allele-specific PCR (KASPar) marker was designed and validated that enables accurate cultivar sex estimation. The reconstruction of chromosome 6 associated with the detection of highly putative structural variations confirmed the possible involvement of a major part of the chromosome.

CONCLUSION

The findings of this study, combined with proper estimation of accession ploidy levels to avoid endosperm incompatibility issues, could facilitate the design of future promising parental combinations in D. alata breeding programs. Moreover, the discovery of this genetic barrier to reproduction opens new avenues for gaining insight into yam reproductive biology and diversification.

Cytological Observations and Bulked-Segregant Analysis Coupled Global Genome Sequencing Reveal Two Genes Associated with Pollen Fertility in Tetraploid Rice

Neo-tetraploid rice with high fertility is a useful germplasm for polyploid rice breeding, which was developed from the crossing of different autotetraploid rice lines. However, little information is available on the molecular mechanism underlying the fertility of neo-tetraploid rice. Here, two contrasting populations of tetraploid rice, including one with high fertility (hereafter referred to as JG) and another with low fertility (hereafter referred to as JD), were generated by crossing Huaduo 3 (H3), a high fertility neo-tetraploid rice that was developed by crossing Jackson-4x with 96025-4x, and Huajingxian74-4x (T452), a low fertility autotetraploid rice parent. Cytological, global genome sequencing-based bulked-segregant (BSA-seq) and CRISPR/Cas9 technology were employed to study the genes associated with pollen fertility in neo-tetraploid rice. The embryo sacs of JG and JD lines were normal; however, pollen fertility was low in JD, which led to scarce fertilization and low seed setting. Cytological observations displayed low pollen fertility (25.1%) and approximately 31.3 and 27.2% chromosome lagging at metaphase I and II, and 28.8 and 24.8% chromosome straggling at anaphase I and II in JD, respectively. BSA-seq of F_2–3 generations and RNA-seq of F₄ generation detected a common fragment, i.e., 18,915,234–19,500,000, at chromosome 7, which was comprised of 78 genes associated with fertility. Among 78 genes, 9 genes had been known to be involved in meiosis and pollen development. Two mutants ny1 (LOC_Os07g32406) and ny2 (LOC_Os07g32040) were generated by CRISPR/Cas9 knockout in neo-tetraploid rice, and which exhibited low pollen fertility and abnormal chromosome behavior. Our study revealed that two unknown genes, LOC_Os07g32406 (NY1) and LOC_Os07g32040 (NY2) play an important role in pollen development of neo-tetraploid rice and provides a new perspective about the genetic mechanisms of fertility in polyploid rice.

Barley Anther and Meiocyte Transcriptome Dynamics in Meiotic Prophase I

In flowering plants, successful germinal cell development and meiotic recombination depend upon a combination of environmental and genetic factors. To gain insights into this specialized reproductive development program we used short- and long-read RNA-sequencing (RNA-seq) to study the temporal dynamics of transcript abundance in immuno-cytologically staged barley (Hordeum vulgare) anthers and meiocytes. We show that the most significant transcriptional changes in anthers occur at the transition from pre-meiosis to leptotene-zygotene, which is followed by increasingly stable transcript abundance throughout prophase I into metaphase I-tetrad. Our analysis reveals that the pre-meiotic anthers are enriched in long non-coding RNAs (lncRNAs) and that entry to meiosis is characterized by their robust and significant down regulation. Intriguingly, only 24% of a collection of putative meiotic gene orthologs showed differential transcript abundance in at least one stage or tissue comparison. Argonautes, E3 ubiquitin ligases, and lys48 specific de-ubiquitinating enzymes were enriched in prophase I meiocyte samples. These developmental, time-resolved transcriptomes demonstrate remarkable stability in transcript abundance in meiocytes throughout prophase I after the initial and substantial reprogramming at meiosis entry and the complexity of the regulatory networks involved in early meiotic processes.

Cell wall Glycine-rich Protein2 is involved in tapetal differentiation and pollen maturation

The tapetum plays important roles in anther development by providing materials for pollen-wall formation and nutrients for pollen development. Here, we report the characterization of a male-sterile mutant of glycine-rich protein 2 (OsGRP2), which exhibits irregular cell division and dysfunction of the tapetum. GRP is a cellwall structural protein present in the cell walls of diverse plant species, but its function is unclear in pollen development. We found that few GRP genes are expressed in rice and thus focused on one highly expressed gene, OsGRP2. The tapetal cell walls of an OsGRP2 mutant did not thicken at the pollen mothercell stage, as a result, pollen maturation and fertility rate decreased. High OsGRP2 expression was detected in male-floral organs, and OsGRP2 was distributed in the tapetum. OsGRP2 participated in establishment of the cellwall network during early tapetum development. In conclusion, our results indicate that OsGRP2 plays important roles in the differentiation and function of the tapetum.

Transcriptome and Gene Editing Analyses Reveal MOF1a Defect Alters the Expression of Genes Associated with Tapetum Development and Chromosome Behavior at Meiosis Stage Resulting in Low Pollen Fertility of Tetraploid Rice

Autotetraploid rice is a useful rice germplasm for polyploid rice breeding. However, low fertility limits its commercial production. A neo-tetraploid rice with high fertility was developed from the progenies of crossing between autotetraploid lines by our research group. Our previous study showed that a myeloblastosis (MYB) transcription factor, MOF1, might be associated with the pollen development in tetraploid rice. However, little information is available about its role in pollen development in tetraploid rice. Here, we identified a new haplotype of MOF1 from neo-tetraploid rice and marked it as MOF1a. Transcriptome and qRT-PCR analysis demonstrated that MOF1a highly expressed in anthers, and displayed differential expression in neo-tetraploid rice compared to tetraploid rice line with low pollen fertility. The mutant (mof1a) of MOF1a, which was generated by CRISPR/Cas9 knockout in neo-tetraploid rice, showed low pollen fertility, and also exhibited abnormal tapetum and middle layer development, and defective chromosome behaviors during meiosis. A total of 13 tapetal related genes were found to be up-regulated in meiotic anthers of MOF1a compared with wild type plants by RNA-seq analysis, including CYP703A3, PTC1, and OsABCG26, which had been demonstrated to affect tapetal development. Moreover, 335 meiosis-related genes displayed differential expression patterns at same stage, including nine important meiosis-related genes, such as metallothionein OsMT1a. These results demonstrated that MOF1a plays an important role in pollen development and provides a foundation for understanding the molecular mechanism underlying MOF1a in reproduction of tetraploid rice.

Comparative study on cytogenetics and transcriptome between diploid and autotetraploid rice hybrids harboring double neutral genes

Double pollen fertility neutral genes, San and Sbn, can control pollen sterility in intersubspecific (indica × japonica) rice hybrids, which has excellent potential to increase rice yield. Previous studies showed that polyploidy could increase the interaction of three pollen sterility loci, i.e. Sa, Sb and Sc, which cause pollen sterility in autotetraploid rice hybrids, and hybrid fertility could be improved by double neutral genes, San and Sbn, in autotetraploid rice hybrids. We compared cytological and transcriptome data between autotetraploid and diploid rice hybrid during meiosis and single microspore stages to understand the molecular mechanism of neutral genes for overcoming pollen sterility in autotetraploid rice hybrids, which harbored double neutral genes. Cytological results revealed that the double neutral genes resulted in higher pollen fertility (76.74%) and lower chromosomal abnormalities in autotetraploid hybrid than in parents during metaphase I, metaphase II, anaphase I and anaphase II. Moreover, autotetraploid rice hybrid displayed stronger heterosis than a diploid hybrid. Compared with diploid rice hybrid, a total of 904 and 68 differently expressed genes (DEGs) were identified explicitly in autotetraploid hybrid at meiosis and single microspore stages, respectively. Of these, 133 and 41 genes were detected in higher-parent dominance and transgressive up-regulation dominance, respectively, which were considered autotetraploid potential heterosis genes, including a meiosis-related gene (Os01g0917500, MSP1) and two meiosis specific-genes (Os07g0624900 and Os04g0208600). Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes pathway (KEGG) analysis revealed that DEGs significantly enriched in amino acid metabolism and photosynthesis metabolism. These results indicated that meiosis-specific and meiosis-related genes, and amino acids and photosynthesis metabolism-related genes contribute to higher yield and pollen fertility in autotetraploid rice hybrid. This study provides a theoretical basis for molecular mechanisms of heterosis in autotetraploid rice harboring double neutral genes for pollen fertility.

Comparative Cytological and Transcriptome Analysis Revealed the Normal Pollen Development Process and Up-Regulation of Fertility-Related Genes in Newly Developed Tetraploid Rice

Autotetraploid rice is a useful germplasm for polyploid rice breeding; however, low seed setting is a major hindrance for its utilization. Here, we reported the development of a new tetraploid rice, Huoduo1 (H1), which has the characteristic of high fertility, from crossing generations of autotetraploid rice. Cytological observations displayed the high fertility of the pollen (95.62%) in H1, a lower percentage of pollen mother cell (PMC) abnormalities, and stable chromosome configurations during the pollen development process compared with its parents. Using RNA-seq analysis, we detected 440 differentially expressed genes (DEGs) in H1 compared with its parents. Of these DEGs, 193 were annotated as pollen fertility-related genes, and 129 (~66.8%) exhibited significant up-regulation in H1 compared with the parents, including three environmentally sensitive genic male sterility genes (TMS9-1, TMS5, and CSA), one meiosis gene (RAD51D), and three tapetal-related genes (MIL2, OsAP25, and OsAP37), which were validated by qRT-PCR in this study. Two genes, TMS9-1 and TMS5, were knocked out using CRISPR/Cas9 technology, and their mutants displayed low fertility and the abnormal development of pollen. Our findings provide evidence for the regulatory mechanisms of fertility in tetraploid rice and indicated that the up-regulation of pollen fertility-related genes may contribute to the high fertility in new tetraploid rice.

Production and roles of IAA and ABA during development of superior and inferior rice grains

Current understanding of the role of plant hormones during cereal grain filling is confounded by contradictory reports on hormone production that is based on poor methodology. We report here on the accurate measurement of indole-3-acetic acid (IAA) and abscisic acid (ABA) by combined liquid chromatography-tandem mass spectrometry in multiple reaction-monitoring mode with heavy isotope labelled internal standards. ABA and IAA contents of superior versus inferior rice grains (ABA maxima 159 ng g-1 FW and 109 ng g-1 FW, IAA maxima 2 µg g-1 FW and 1.7 µg g-1 FW respectively) correlated with the expression of biosynthetic genes and with grain fill. Results confirm that grain ABA is produced primarily by OsNCED2(5), but suggest that ABA import and metabolism also play important roles in ABA regulation. The IAA content of grains is primarily influenced by OsYUC9 and OsYUC11. However, the distinct expression profile of OsYUC12 suggests a specific role for IAA produced by this enzyme. Co-expression of OsYUC12 with OsIAA29 indicates their involvement in a common signalling pathway. Co-expression and cis-element analysis identified several aleurone-specific transcriptional regulators as well as glutelin as strong candidates for detailed investigation for direct regulation by the auxin-signalling pathway.

Global identification and analysis revealed differentially expressed lncRNAs associated with meiosis and low fertility in autotetraploid rice

BACKGROUND

Autotetraploid rice is a useful germplasm for polyploid rice breeding. Our previous research showed that non-coding RNAs might be associated with low fertility in autotetraploid rice. However, little information is available on long non-coding RNAs (lncRNAs) involved in the low fertility of autotetraploid rice. In the present study, RNA-seq was employed to detect the differentially expressed meiosis-related lncRNAs in autotetraploid rice, and gene overexpression and knock out experiments were used to validate the potential function of candidate lncRNA.

RESULTS

A total of 444 differentially expressed lncRNAs (DEL) were detected during anther and ovary meiosis in autotetraploid rice. Of these, 328 DEL were associated with the transposable elements, which displayed low expression levels during meiosis in autotetraploid rice. We used rapid amplification of cDNA ends (RACE) assay to validate 10 DEL and found that the lncRNAs were not assembly artifacts, and six of them were conserved in tetraploid rice. Moreover, 237 and 20 lncRNAs were associated with pollen mother cell (PMC) and embryo sac mother cell (EMC) meiosis in autotetraploid rice, respectively. The differential expressions of some meiosis-related targets and its DEL regulator, including MEL1 regulated by TCONS_00068868, LOC_Os12g41350 (meiotic asynaptic mutant 1) by TCONS_00057811 in PMC, and LOC_Os12g39420 by TCONS_00144592 in EMC, were confirmed by qRT-PCR. TCONS_00057811, TCONS_00055980 and TCONS_00130461 showed anther specific expression patterns and were found to be highly expressed during meiosis. CRISPR/Cas9 editing of lncRNA57811 displayed similar morphology compared to wild type. The overexpression of lncRNA57811 resulted in low pollen fertility (29.70%) and seed setting (33%) in rice.

CONCLUSION

The differential expression levels of lncRNAs, associated with transposable elements and meiosis-regulated targets, might be endogenous noncoding regulators of pollen/embryo sac development that cause low fertility in autotetraploid rice. The results enhance our understanding about rice lncRNAs, and facilitate functional research in autotetraploid rice.

Transcriptional trajectories of anther development provide candidates for engineering male fertility in sorghum

Sorghum is a self-pollinated crop with multiple economic uses as cereal, forage, and biofuel feedstock. Hybrid breeding is a cornerstone for sorghum improvement strategies that currently relies on cytoplasmic male sterile lines. To engineer genic male sterility, it is imperative to examine the genetic components regulating anther/pollen development in sorghum. To this end, we have performed transcriptomic analysis from three temporal stages of developing anthers that correspond to meiotic, microspore and mature pollen stages. A total of 5286 genes were differentially regulated among the three anther stages with 890 of them exhibiting anther-preferential expression. Differentially expressed genes could be clubbed into seven distinct developmental trajectories using K-means clustering. Pathway mapping revealed that genes involved in cell cycle, DNA repair, regulation of transcription, brassinosteroid and auxin biosynthesis/signalling exhibit peak expression in meiotic anthers, while those regulating abiotic stress, carbohydrate metabolism, and transport were enriched in microspore stage. Conversely, genes associated with protein degradation, post-translational modifications, cell wall biosynthesis/modifications, abscisic acid, ethylene, cytokinin and jasmonic acid biosynthesis/signalling were highly expressed in mature pollen stage. High concurrence in transcriptional dynamics and cis-regulatory elements of differentially expressed genes in rice and sorghum confirmed conserved developmental pathways regulating anther development across species. Comprehensive literature survey in conjunction with orthology analysis and anther-preferential accumulation enabled shortlisting of 21 prospective candidates for in-depth characterization and engineering male fertility in sorghum.

Cytological and transcriptome analyses reveal OsPUB73 defect affects the gene expression associated with tapetum or pollen exine abnormality in rice

BACKGROUND

As one of the main crops in the world, sterility of rice (Oryza sativa L.) significantly affects the production and leads to yield decrease. Our previous research showed that OsPUB73, which encodes U-box domain-containing protein 73, may be associated with male sterility. However, little information is available on this gene that is required for anther development. In the present study, we knocked out OsPUB73 by using the CRISPR/Cas9 system and studied the cytological and transcriptome of the gene-defect associated with pollen development and sterility in the rice variety (Taichung 65).

RESULTS

The sequence analysis indicated that OsPUB73 was comprised of 3 exons and 2 introns, of which CDS encoded 586 amino acids including a U-box domain. The expression pattern of OsPUB73 showed that it was highly expressed in the anther during meiosis stage. The ospub73 displayed low pollen fertility (19.45%), which was significantly lower than wild type (WT) (85.37%). Cytological observation showed tapetum vacuolated at the meiosis stage and pollen exine was abnormal at the bi-cellular pollen stage of ospub73. RNA-seq analysis detected 2240 down and 571 up-regulated genes in anther of ospub73 compared with WT during meiosis stage. Among of 2240 down-regulated genes, seven known genes were associated with tapetal cell death or pollen exine development, including CYP703A3 (Cytochrome P450 Hydroxylase703A3), CYP704B2 (Cytochrome P450 Hydroxylase704B2), DPW (Defective Pollen Wall), PTC1 (Persistant Tapetal Cell1), UDT1 (Undeveloped Tapetum1), OsAP37 (Aspartic protease37) and OsABCG15 (ATP binding cassette G15), which were validated by quantitative real-time polymerase chain reaction (qRT-PCR). These results suggested OsPUB73 may play an important role in tapetal or pollen exine development and resulted in pollen partial sterility.

CONCLUSION

Our results revealed that OsPUB73 plays an important role in rice male reproductive development, which provides valuable information about the molecular mechanisms of the U-box in rice male reproductive development.

Generation of a Comprehensive Transcriptome Atlas and Transcriptome Dynamics in Medicinal Cannabis

Cannabinoids are the main medicinal compounds of interest in the plant Cannabis sativa, that are primarily synthesised in the glandular trichomes; found on female floral buds. The content, composition and yield of secondary metabolites (cannabinoids and terpenoids) is influenced by the plant's genetics and environment. Some initial gene expression experiments have been performed from strains of this plant species that contrasted in cannabinoid production, however the present knowledge about detailed trichome transcriptomics in this species is limited. An extensive transcriptome atlas was generated by RNA sequencing using root, shoot, flower and trichome tissues from a female plant strain (Cannbio-2) and was enhanced with the addition of vegetative and reproductive tissues from a male cannabis plant. Differential gene expression analysis identified genes preferentially expressed in different tissues. Detailed trichomics was performed from extractions specifically from glandular trichomes as well as female floral tissues at varying developmental stages, to identify stage-specific differentially expressed genes. Candidate genes involved in terpene and cannabinoid synthesis were identified and the majority were found to have an abundant expression in trichomes. The comprehensive transcriptome is a significant resource in cannabis for further research of functional genomics to improve the yield of specialised metabolites with high pharmacological value.

Cytological and transcriptome analysis reveal that interaction at Sb pollen sterility locus cause down-regulation of important meiosis-related genes associated with high pollen sterility in autotetraploid rice hybrids

Polyploidy could increase the interactions of pollen sterility loci and Sb locus interaction cause higher pollen abortion than other loci. Therefore, we focused on the interaction at Sb pollen sterility locus in autotetraploid rice compared to diploid rice hybrid using the near-isogenic lines in the present study. Cytological observations indicated that interaction at Sb locus cause high pollen sterility (69.9%) and abnormal chromosome behavior (37.02%) at Metaphase II in autotetraploid rice hybrid. A total of 139 meiosis-related or meiosis stage-specific genes were detected in the autotetraploid rice hybrid harboring interaction at Sb locus and 27 of these meiosis-related or specific genes displayed significant down-regulation, including four pollen fertility related genes (Rad51, XRI1, PSS1 and MIL1). These results revealed a stronger interaction at Sb pollen sterility locus than other loci, which cause down-regulation of many important meiosis-related genes that were associated with higher pollen sterility in autotetraploid rice hybrids.

Transcriptome analysis reveals potential roles of a barley ASR gene that confers stress tolerance in transgenic rice

Control of gene expression and induction of cellular protection mechanisms are two important processes that plants employ to protect themselves against abiotic stresses. ABA-, stress, and ripening-induced (ASR) proteins have been identified to participate in such responses. Previous studies have proposed that these proteins can act as transcription factors and as molecular chaperones protecting transgenic plants against stresses; however a gene network regulated by ASRs has not been explored. To expand our knowledge on the function of these proteins in cereals, we present the functional characterization of a barley ASR gene. Expression of HvASR5 was almost ubiquitous in different organs and responded to ABA and to different stress treatments. When expressed ectopically, HvASR5 was able to confer drought and salt stress tolerance to Arabidopsis thaliana and to improve growth performance of rice plants under stress conditions. A transcriptomic analysis with two transgenic rice lines overexpressing HvASR5 helped to identify potential downstream targets and understand ASR-regulated cellular processes. HvASR5 up-regulated the expression of a distinct set of genes associated with stress responses, metabolic processes (particularly carbohydrate metabolism), as well as reproduction and development. These data, together with the confirmed nuclear and cytoplasmic localization of HvASR5, further support the hypothesis that HvASR5 can also carry out roles as molecular protector and transcriptional regulator.

Carbohydrate metabolism and fertility related genes high expression levels promote heterosis in autotetraploid rice harboring double neutral genes

BACKGROUND

Autotetraploid rice hybrids have great potential to increase the production, but hybrid sterility is a major hindrance in the utilization of hybrid vigor in polyploid rice, which is mainly caused by pollen abortion. Our previous study showed that double pollen fertility neutral genes, Sa-n and Sb-n, can overcome hybrid sterility in autotetraploid rice. Here, we used an autotetraploid rice line harboring double neutral genes to develop hybrids by crossing with auto- and neo-tetraploid rice, and evaluated heterosis and its underlying molecular mechanism.

RESULTS

All autotetraploid rice hybrids, which harbored double pollen fertility neutral genes, Sa-n and Sb-n, displayed high seed setting and significant positive heterosis for yield and yield-related traits. Cytological observations revealed normal chromosome behaviors and higher frequency of bivalents in the hybrid than parents during meiosis. Transcriptome analysis revealed significantly higher expressions of important saccharides metabolism and starch synthase related genes, such as OsBEIIb and OsSSIIIa, in the grains of hybrid than parents. Furthermore, many meiosis-related and specific genes, including DPW and CYP703A3, displayed up-regulation in the hybrid compared to a parent with low seed setting. Many non-additive genes were detected in the hybrid, and GO term of carbohydrate metabolic process was significantly enriched in all the transcriptome tissues except flag leaf (three days after flowering). Moreover, many differentially expressed genes (DEGs) were identified in the yield-related quantitative trait loci (QTLs) regions as possible candidate genes.

CONCLUSION

Our results revealed that increase in the number of bivalents improved the seed setting of hybrid harboring double pollen fertility neutral genes. Many important genes, including meiosis-related and meiosis-specific genes and saccharides metabolism and starch synthase related genes, exhibited heterosis specific expression patterns in polyploid rice during different development stages. The functional analysis of important genes will provide valuable information for molecular mechanisms of heterosis in polyploid rice.

Genome-wide analysis of DNA polymorphisms, the methylome and transcriptome revealed that multiple factors are associated with low pollen fertility in autotetraploid rice

Autotetraploid rice is a useful germplasm with high biomass production; however, low fertility is the main barrier in commercial utilization. In our previous study, differential expression of meiosis-related miRNAs was found to be involved in the pollen sterility of autotetraploid rice. However, genome-wide DNA variations and methylomes associated with low fertility of autotetraploid rice are still poorly understood. Here, we measured both global DNA variations and the methylome and compared them with the transcriptome during pollen development in autotetraploid rice by high-throughput sequencing. A total of 34416 SNPs, 6993 InDels, 1003 SVs and 25 CNVs were detected, and 11367 and 41117 differentially methylated regions showed hypermethylation and hypomethylation in 02428-4x. In total, 1110 genes displayed differentially expression in 02428-4x during meiosis, of these six harbored CNVs, including four upregulated genes with gain CNVs, such as LOC_Os11g38620. We identified 122 genes by comparing with the previous data that might be associated with low fertility during pollen development in 02428-4x. Of the 122 gens, 98 were displayed methylation and differential expression, including OsMADS98, CYP703A3 and OsABCG26. The downregulation of these three genes were confirmed by qPCR during meiosis of 02428-4x, which played pivotal roles in pollen fertility. These results indicate that the low fertility of autotetraploid rice is not only caused by the differential expression of genes involved in pollen development, but also by sequence variation and differential methylation, suggesting that the reason for pollen sterility in autotetraploid rice is complex and might be affected by multiple factors.

Cytological and transcriptome analyses reveal abrupt gene expression for meiosis and saccharide metabolisms that associated with pollen abortion in autotetraploid rice

Autotetraploid rice is a useful germplasm that has four chromosome sets and strong biological advantages; however, low fertility limits its commercial utilization. Little information is available about the DNA variation and differential gene expressions associated with low fertility in autotetraploid rice. In the present study, 81 SNPs and 182 InDels were identified in T449 (an autotetraploid rice line with low fertility) compared to E249 (diploid counterpart) by whole-genome re-sequencing. We detected only three non-synonymous SNPs and six large-effect InDels, which were associated with three and six genes, respectively. A total of 75 meiosis-related differentially expressed genes were detected during the meiosis stage by transcriptome analysis, including OsMTOPVIB, which is essential for meiotic DSB formation, and OsMOF, which takes part in homologous chromosome pairing and synapsis. Approximately 20.69% lagging chromosome at metaphase I and 4.65% abnormal tetrad were observed in T449. Moreover, transcriptome analysis revealed down-regulation of a sucrose transporter (OsSUT5) and two monosaccharide transporters (OsMST1 and OsMST8) in T449 at the single microspore stage, and their expression levels were verified by qRT-PCR. Cytological observation of saccharide distribution showed abnormal accumulation of saccharides in T449 and the contents of fructose and glucose were markedly higher in T449 than E249 at the single microspore stage. Our results suggested that polyploidy not only induces abrupt expression changes in the meiosis-related genes that lead to abnormal chromosome behavior, but also causes changes in the saccharide distribution and expression patterns of saccharide-related genes, which jointly causes sterility in the autotetraploid rice.

Effects of threshold on the topology of gene co-expression networks

Several developments regarding the analysis of gene co-expression profiles using complex network theory have been reported recently. Such approaches usually start with the construction of an unweighted gene co-expression network, therefore requiring the selection of a suitable threshold defining which pairs of vertices will be connected. We aimed at addressing such an important problem by suggesting and comparing five different approaches for threshold selection. Each of the methods considers a respective biologically-motivated criterion for electing a potentially suitable threshold. A set of 21 microarray experiments from different biological groups was used to investigate the effect of applying the five proposed criteria to several biological situations. For each experiment, we used the Pearson correlation coefficient to measure the relationship between each gene pair, and the resulting weight matrices were thresholded considering several values, generating respective adjacency matrices (co-expression networks). Each of the five proposed criteria was then applied in order to select the respective threshold value. The effects of these thresholding approaches on the topology of the resulting networks were compared by using several measurements, and we verified that, depending on the database, the impact on the topological properties can be large. However, a group of databases was verified to be similarly affected by most of the considered criteria. Based on such results, it can be suggested that when the generated networks present similar measurements, the thresholding method can be chosen with greater freedom. If the generated networks are markedly different, the thresholding method that better suits the interests of each specific research study represents a reasonable choice.

Genome-wide analyses of late pollen-preferred genes conserved in various rice cultivars and functional identification of a gene involved in the key processes of late pollen development

BACKGROUND

Understanding late pollen development, including the maturation and pollination process, is a key component in maintaining crop yields. Transcriptome data obtained through microarray or RNA-seq technologies can provide useful insight into those developmental processes. Six series of microarray data from a public transcriptome database, the Gene Expression Omnibus of the National Center for Biotechnology Information, are related to anther and pollen development.

RESULTS

We performed a systematic and functional study across the rice genome of genes that are preferentially expressed in the late stages of pollen development, including maturation and germination. By comparing the transcriptomes of sporophytes and male gametes over time, we identified 627 late pollen-preferred genes that are conserved among japonica and indica rice cultivars. Functional classification analysis with a MapMan tool kit revealed a significant association between cell wall organization/metabolism and mature pollen grains. Comparative analysis of rice and Arabidopsis demonstrated that genes involved in cell wall modifications and the metabolism of major carbohydrates are unique to rice. We used the GUS reporter system to monitor the expression of eight of those genes. In addition, we evaluated the significance of our candidate genes, using T-DNA insertional mutant population and the CRISPR/Cas9 system. Mutants from T-DNA insertion and CRISPR/Cas9 systems of a rice gene encoding glycerophosphoryl diester phosphodiesterase are defective in their male gamete transfer.

CONCLUSION

Through the global analyses of the late pollen-preferred genes from rice, we found several biological features of these genes. First, biological process related to cell wall organization and modification is over-represented in these genes to support rapid tube growth. Second, comparative analysis of late pollen preferred genes between rice and Arabidopsis provide a significant insight on the evolutional disparateness in cell wall biogenesis and storage reserves of pollen. In addition, these candidates might be useful targets for future examinations of late pollen development, and will be a valuable resource for accelerating the understanding of molecular mechanisms for pollen maturation and germination processes in rice.

Pervasive interactions of Sa and Sb loci cause high pollen sterility and abrupt changes in gene expression during meiosis that could be overcome by double neutral genes in autotetraploid rice

BACKGROUND

Intersubspecific autotetraploid rice hybrids possess high hybrid vigor; however, low pollen fertility is a critical hindrance in its commercial utilization. Our previous study demonstrated that polyploidy could increase the multi-loci interaction and cause high pollen abortion in autotetraploid rice hybrids. However, there is little known about the critical role of pollen sterility locus or loci in the intersubspecific hybrids. We developed autotetraploid rice hybrids harboring heterozygous genotypes (S i S i S j S j ) at different pollen sterility loci by using the near isogenic lines of Taichung65-4×. Moreover, autotetraploid lines carrying double neutral genes, Sa n and Sb n , were used to assess their effect on fertility restoration.

RESULTS

Cytological studies showed that the deleterious genetic interactions at Sa and Sb pollen sterility loci resulted in higher pollen sterility (76.83%) and abnormal chromosome behavior (24.59%) at metaphase I of meiosis in autotetraploid rice hybrids. Transcriptome analysis revealed 1092 differentially expressed genes (DEG) in a hybrid with the pervasive interactions at Sa and Sb pollen sterility loci, and most of the genes (about 83%) exhibited down regulation. Of the DEG, 60 were associated with transcription regulation and 18 genes were annotated as meiosis-related genes. Analysis on the hybrids developed by using autotetraploid rice harboring double neutral genes, Sa n and Sb n , revealed normal pollen fertility, and transcriptome analysis showed non-significant difference in number of DEG among different hybrids.

CONCLUSIONS

Our finding revealed that pervasive interactions at Sa and Sb pollen sterility loci cause high sterility in the autotetraploid hybrids that lead to the down-regulation of important meiosis-related genes and transcription regulation factors. Moreover, we also found that the hybrids sterility could be overcome by double neutral genes, Sa n and Sb n , in autotetraploid rice hybrids. The present study provided a strong evidence for the utilization of heterosis in autotetraploid rice hybrids.

RiceAntherNet: a gene co-expression network for identifying anther and pollen development genes

In plants, normal anther and pollen development involves many important biological events and complex molecular regulatory coordination. Understanding gene regulatory relationships during male reproductive development is essential for fundamental biology and crop breeding. In this work, we developed a rice gene co-expression network for anther development (RiceAntherNet) that allows prediction of gene regulatory relationships during pollen development. RiceAntherNet was generated from 57 rice anther tissue microarrays across all developmental stages. The microarray datasets from nine rice male sterile mutants, including msp1-4, ostdl1a, gamyb-2, tip2, udt1-1, tdr, eat1-1, ptc1 and mads3-4, were used to explore and test the network. Among the changed genes, three clades showing differential expression patterns were constructed to identify genes associated with pollen formation. Many of these have known roles in pollen development, for example, seven genes in Clade 1 (OsABCG15, OsLAP5, OsLAP6, DPW, CYP703A3, OsNP1 and OsCP1) are involved in rice pollen wall formation. Furthermore, Clade 1 contained 12 genes whose predicted orthologs in Arabidopsis have been reported as key during pollen development and may play similar roles in rice. Genes in Clade 2 are expressed earlier than Clade 1 (anther stages 2-9), while genes in Clade 3 are expressed later (stages 10-12). RiceAntherNet serves as a valuable tool for identifying novel genes during plant anther and pollen development. A website is provided (https://www.cpib.ac.uk/anther/riceindex.html) to present the expression profiles for gene characterization. This will assist in determining the key relationships between genes, thus enabling characterization of critical genes associated with anther and pollen regulatory networks.

miR3954 is a trigger of phasiRNAs that affects flowering time in citrus

In plant, a few 22-nt miRNAs direct cleavages of their targets and trigger the biogenesis of phased small interfering RNAs (phasiRNAs) in plant. In this study, we characterized a miRNA triggering phasiRNAs generation, miR3954, and explored its downstream target genes and potential function. Our results demonstrated that miR3954 showed specific expression in the flowers of citrus species, and it targeted a NAC transcription factor (Cs7 g22460) and two non-coding RNA transcripts (lncRNAs, Cs1 g09600 and Cs1 g09635). The production of phasiRNAs was detected from transcripts targeted by miR3954, and was further verified in both sequencing data and transient expression experiments. PhasiRNAs derived from the two lncRNAs targeted not only miR3954-targeted NAC gene but also additional NAC homologous genes. No homologous genes of these two lncRNAs were found in plants other than citrus species, implying that this miR3954-lncRNAs-phasiRNAs-NAC pathway is likely citrus-specific. Transgenic analysis indicated that the miR3954-overexpressing lines showed decreased transcripts of lncRNA, elevated abundance of phasiRNAs and reduced expression of NAC genes. Interestingly, the overexpression of miR3954 leads to early flowering in citrus plants. In summary, our results illustrated a model of the regulatory network of miR3954-lncRNA-phasiRNAs-NAC, which may be functionally involved in flowering in citrus.

Tobacco rattle virus (TRV)-Mediated Silencing of Nicotiana benthamiana ARGONAUTES (NbAGOs) Reveals New Antiviral Candidates and Dominant Effects of TRV-NbAGO1

The objective of this study was to determine the contribution of different ARGONAUTE proteins in Nicotiana benthamiana (NbAGOs) to the defense against silencing sensitive GFP-expressing viral constructs based on Tomato bushy stunt virus (TBSV) (Tombusvirus), Sunn-hemp mosaic virus (Tobamovirus), and Foxtail mosaic virus (Potexvirus). Upon Tobacco rattle virus (TRV)-mediated down-regulation of NbAGO1, 4, 5, or 6, no effects were noted on susceptibility to any virus construct, whereas knockdown of NbAGO2 specifically prevented silencing of P19-defective TBSV (TGdP19). Down-regulation of a new gene referred to as NbAGO5L showed some reduced silencing for TGdP19 but not for the other two virus constructs, whereas silencing of NbAGO7 gave rise to a subtle increase in susceptibility to all three viruses. Co-infiltrating different TRV-NbAGO constructs simultaneously did not enhance virus susceptibility. However, an unexpected finding was that whenever the TRV-NbAGO1 construct was present, this compromised silencing of genes targeted by co-infiltrated constructs, as shown upon co-infiltration of TRV-NbAGO1 with either TRV-NbAGO2 or TRV-Sul (targeting Magnesium chelatase I). Only after a prolonged period (approximately 2 months) did TRV-Sul-mediated systemic bleaching occur in these co-infected plants, suggesting that TRV-NbAGO1 hinders the silencing ability of other TRV-NbAGO constructs. In conclusion, this study revealed new antiviral NbAGOs and dominant effects of silencing NbAGO1.

Anther-preferential expressing gene PMR is essential for the mitosis of pollen development in rice

Phenotype identification, expression examination, and function prediction declared that the anther-preferential expressing gene PMR may participate in regulation of male gametophyte development in rice. Male germline development in flowering plants produces the pair of sperm cells for double fertilization and the pollen mitosis is a key process of it. Although the structural features of male gametophyte have been defined, the molecular mechanisms regulating the mitotic cell cycle are not well elucidated in rice. Here, we reported an anther-preferential expressing gene in rice, PMR (Pollen Mitosis Relative), playing an essential role in male gametogenesis. When PMR gene was suppressed via RNAi, the mitosis of microspore was severely damaged, and the plants formed unmatured pollens containing only one or two nucleuses at the anthesis, ultimately leading to serious reduction of pollen fertility and seed-setting. The CRISPR mutants, pmr-1 and pmr-2, both showed the similar defects as the PMR-RNAi lines. Further analysis revealed that PMR together with its co-expressing genes were liable to participate in the regulation of DNA metabolism in the nucleus, and affected the activities of some enzymes related to the cell cycle. We finally discussed that unknown protein PMR contained the PHD, SWIB and Plus-3 domains and they might have coordinating functions in regulation pathway of the pollen mitosis in rice.

A long noncoding RNA involved in rice reproductive development by negatively regulating osa-miR160

Long noncoding RNAs (lncRNAs) participate in the regulation of multiple biological processes via diverse manners, one of which is functioning as endogenous target mimics (eTMs) to modulate microRNAs (miRNAs) by competing for their targets. Previously, we have predicted one lncRNA (osa-eTM160) as an endogenous repressor of osa-miR160 and validated the target mimicry ability of osa-eTM160 for ath-miR160 in Arabidopsis thaliana, yet the functions of osa-eTM160 in rice remain obscure. Here, we demonstrated that osa-eTM160 attenuated the repression of osa-miR160 on osa-ARF18 mRNAs during early anther developmental stages through the target mimicry manner, therefore to regulate rice seed setting and seed size. These findings revealed the roles of osa-eTM160 in rice, and indicated that lncRNAs with eTM functions may serve as temporal regulators to modulate the effects of miRNAs at specific developmental stages.

Analysis of small RNAs revealed differential expressions during pollen and embryo sac development in autotetraploid rice

Background

Partial pollen and embryo sac sterilities are the two main reasons for low fertility in autotetraploid rice. Our previous study revealed that small RNAs changes may associate with pollen fertility in autotetraploid rice. However, knowledge on comparative analysis between the development of pollen and embryo sac by small RNAs in autotetraploid rice is still unknown. In the present study, WE-CLSM (whole-mount eosin B-staining confocal laser scanning microscopy) and high-throughput sequencing technology was employed to examine the cytological variations and to analyze small RNAs changes during pollen and embryo sac development in autotetraploid rice compared with its diploid counterpart.

Results

A total of 321 and 368 differentially expressed miRNAs (DEM) were detected during pollen and embryo sac development in autotetraploid rice, respectively. Gene Ontology enrichment analysis on the targets of DEM associated with embryo sac and pollen development revealed 30 prominent functional gene classes, such as cell differentiation and signal transduction during embryo sac development, while only 7 prominent functional gene classes, such as flower development and transcription factor activity, were detected during pollen development in autotetraploid rice. The expression levels of 39 DEM, which revealed interaction with meiosis-related genes, showed opposite expression patterns during pollen and embryo sac development. Of these DEM, osa-miR1436_L + 3_1ss5CT and osa-miR167h-3p were associated with the female meiosis, while osa-miR159a.1 and osa-MIR159a-p5 were related with the male meiosis. 21 nt-phasiRNAs were detected during both pollen and embryo sac development, while 24 nt-phasiRNAs were found only in pollen development, which displayed down-regulation in autotetraploid compared to diploid rice and their spatial-temporal expression patterns were similar to osa-miR2275d. 24 nt TEs-siRNAs were found to be up-regulated in embryo sac but down-regulated in pollen development.

Conclusion

The above results not only provide the small RNAs changes during four landmark stages of pollen and embryo sac development in autotetraploid rice but also have identified specifically expressed miRNAs, especially meiosis-related miRNAs, pollen-specific-24 nt-phasiRNAs and TEs-siRNAs in autotetraploid rice. Together, these findings provide a foundation for understanding the effect of polyploidy on small RNAs expression patterns during pollen and embryo sac development that may lead to different abnormalities in autotetraploid rice.

Electronic supplementary material

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

Genome-wide association analysis for heat tolerance at flowering detected a large set of genes involved in adaptation to thermal and other stresses

Fertilization sensitivity to heat in rice is a major issue within climate change scenarios in the tropics. A panel of 167 indica landraces and improved varieties was phenotyped for spikelet sterility (SPKST) under 38°C during anthesis and for several secondary traits potentially affecting panicle micro-climate and thus the fertilization process. The panel was genotyped with an average density of one marker per 29 kb using genotyping by sequencing. Genome-wide association analyses (GWAS) were conducted using three methods based on single marker regression, haplotype regression and simultaneous fitting of all markers, respectively. Fourteen loci significantly associated with SPKST under at least two GWAS methods were detected. A large number of associations was also detected for the secondary traits. Analysis of co-localization of SPKST associated loci with QTLs detected in progenies of bi-parental crosses reported in the literature allowed to narrow -down the position of eight of those QTLs, including the most documented one, qHTSF4.1. Gene families underlying loci associated with SPKST corresponded to functions ranging from sensing abiotic stresses and regulating plant response, such as wall-associated kinases and heat shock proteins, to cell division and gametophyte development. Analysis of diversity at the vicinity of loci associated with SPKST within the rice three thousand genomes, revealed widespread distribution of the favourable alleles across O. sativa genetic groups. However, few accessions assembled the favourable alleles at all loci. Effective donors included the heat tolerant variety N22 and some Indian and Taiwanese varieties. These results provide a basis for breeding for heat tolerance during anthesis and for functional validation of major loci governing this trait.

Transcriptome analysis of neo-tetraploid rice reveals specific differential gene expressions associated with fertility and heterosis

Polyploid rice hybrids have a powerful biological and yield potential that may become a new way for rice breeding; however, low fertility is major hindrance in commercial utilization. Here, we developed a neo-tetraploid rice that could overcome the sterility of autotetraploid rice and produce high heterosis. Transcriptome analysis of F₁ hybrid developed by crossing neo-tetraploid with autotetraploid rice displayed 807, 663 and 866 differentially expressed genes that uniquely associated with F₁ and specific to (DEGFu-sp) anther, ovary and leaf, respectively. Of the DEGFu-sp, 1224 genes displayed nonadditive expression; 44 and 10 genes were annotated as TFs and methyltransferase or hydroxymethyltransferase, respectively. Gene ontology enrichment and co-expression analysis revealed specific differential gene expressions in the DEGFu-sp to leaf, anther and ovary, such as genes related to photosynthesis, metabolic process and transport, and co-expression network including fertility, resistance and epigenetic elements. Of the DEGFu-sp to anther, 42 meiosis stage-specific genes, eight meiosis-related genes, such as RAD51 and SMC2, were identified. We identified 38 miRNAs from DEGFu-sp to anther, and their targets were associated with pollen fertility and retrotransposon protein. Our study provides new germplasm for polyploid rice breeding, and revealed complex regulatory mechanisms that might be associated with heterosis and fertility.

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

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

Application of rice microspore-preferred promoters to manipulate early pollen development in Arabidopsis: a heterologous system

Rice microspore-promoters. Based on microarray data analyzed for developing anthers and pollen grains, we identified nine rice microspore-preferred (RMP) genes, designated RMP1 through RMP9. To extend their biotechnological applicability, we then investigated the activity of RMP promoters originating from monocotyledonous rice in a heterologous system of dicotyledonous Arabidopsis. Expression of GUS was significantly induced in transgenic plants from the microspore to the mature pollen stages and was driven by the RMP1, RMP3, RMP4, RMP5, and RMP9 promoters. We found it interesting that, whereas RMP2 and RMP6 directed GUS expression in microspore at the early unicellular and bicellular stages, RMP7 and RMP8 seemed to be expressed at the late tricellular and mature pollen stages. Moreover, GUS was expressed in seven promoters, RMP3 through RMP9, during the seedling stage, in immature leaves, cotyledons, and roots. To confirm microspore-specific expression, we used complementation analysis with an Arabidopsis male-specific gametophytic mutant, sidecar pollen-2 (scp-2), to verify the activity of three promoters. That mutant shows defects in microspore development prior to pollen mitosis I. These results provide strong evidence that the SIDECAR POLLEN gene, driven by RMP promoters, successfully complements the scp-2 mutation, and they strongly suggest that these promoters can potentially be applied for manipulating the expression of target genes at the microspore stage in various species.

Identification and molecular characterization of the nicotianamine synthase gene family in bread wheat

Nicotianamine (NA) is a non-protein amino acid involved in fundamental aspects of metal uptake, transport and homeostasis in all plants and constitutes the biosynthetic precursor of mugineic acid family phytosiderophores (MAs) in graminaceous plant species. Nicotianamine synthase (NAS) genes, which encode enzymes that synthesize NA from S-adenosyl-L-methionine (SAM), are differentially regulated by iron (Fe) status in most plant species and plant genomes have been found to contain anywhere from 1 to 9 NAS genes. This study describes the identification of 21 NAS genes in the hexaploid bread wheat (Triticum aestivum L.) genome and their phylogenetic classification into two distinct clades. The TaNAS genes are highly expressed during germination, seedling growth and reproductive development. Fourteen of the clade I NAS genes were up-regulated in root tissues under conditions of Fe deficiency. Protein sequence analyses revealed the presence of endocytosis motifs in all of the wheat NAS proteins as well as chloroplast, mitochondrial and secretory transit peptide signals in four proteins. These results greatly expand our knowledge of NAS gene families in graminaceous plant species as well as the genetics underlying Fe nutrition in bread wheat.

Rice caryopsis development I: Dynamic changes in different cell layers

Rice caryopsis as one of the most important food sources for humans has a complex structure that is composed of maternal tissues including the pericarp and testa and filial tissues including the endosperm and embryo. Although rice caryopsis studies have been conducted previously, a systematic characterization throughout the entire developmental process is still lacking. In this study, detailed morphological examinations of caryopses were made during the entire 30-day developmental process. We observed some rapid changes in cell differentiation events and cataloged how cellular degeneration processes occurred in maternal tissues. The differentiations of tube cells and cross cells were achieved by 9 days after pollination (DAP). In the testa, the outer integument was degenerated by 3 DAP, while the outer layer of the inner integument degenerated by 7 DAP. In the nucellus, all tissues with the exception of the nucellar projection and the nucellar epidermis degenerated in the first 5 DAP. By 21 DAP, all maternal tissues, including vascular bundles, the nucellar projection and the nucellar epidermal cells were degenerated. In summary, this study provides a complete atlas of the dynamic changes in cell differentiation and degeneration for individual maternal cell layers of rice caryopsis.

Genome-wide identification and analysis of rice genes preferentially expressed in pollen at an early developmental stage

Microspore production using endogenous developmental programs has not been well studied. The main limitation is the difficulty in identifying genes preferentially expressed in pollen grains at early stages. To overcome this limitation, we collected transcriptome data from anthers and microspore/pollen and performed meta-expression analysis. Subsequently, we identified 410 genes showing preferential expression patterns in early developing pollen samples of both japonica and indica cultivars. The expression patterns of these genes are distinguishable from genes showing pollen mother cell or tapetum-preferred expression patterns. Gene Ontology enrichment and MapMan analyses indicated that microspores in rice are closely linked with protein degradation, nucleotide metabolism, and DNA biosynthesis and regulation, while the pollen mother cell or tapetum are strongly associated with cell wall metabolism, lipid metabolism, secondary metabolism, and RNA biosynthesis and regulation. We also generated transgenic lines under the control of the promoters of eight microspore-preferred genes and confirmed the preferred expression patterns in plants using the GUS reporting system. Furthermore, cis-regulatory element analysis revealed that pollen specific elements such as POLLEN1LELAT52, and 5659BOXLELAT5659 were commonly identified in the promoter regions of eight rice genes with more frequency than estimation. Our study will provide new sights on early pollen development in rice, a model crop plant.

Transcriptomes of Plant Gametophytes Have a Higher Proportion of Rapidly Evolving and Young Genes than Sporophytes

Reproductive traits in plants tend to evolve rapidly due to various causes that include plant-pollinator coevolution and pollen competition, but the genomic basis of reproductive trait evolution is still largely unknown. To characterize evolutionary patterns of genome wide gene expression in reproductive tissues in the gametophyte and to compare them to developmental stages of the sporophyte, we analyzed evolutionary conservation and genetic diversity of protein-coding genes using microarray-based transcriptome data from three plant species, Arabidopsis thaliana, rice (Oryza sativa), and soybean (Glycine max). In all three species a significant shift in gene expression occurs during gametogenesis in which genes of younger evolutionary age and higher genetic diversity contribute significantly more to the transcriptome than in other stages. We refer to this phenomenon as "evolutionary bulge" during plant reproductive development because it differentiates the gametophyte from the sporophyte. We show that multiple, not mutually exclusive, causes may explain the bulge pattern, most prominently reduced tissue complexity of the gametophyte, a varying extent of selection on reproductive traits during gametogenesis as well as differences between male and female tissues. This highlights the importance of plant reproduction for understanding evolutionary forces determining the relationship of genomic and phenotypic variation in plants.

Ascribing Functions to Genes: Journey Towards Genetic Improvement of Rice Via Functional Genomics

Rice, one of the most important cereal crops for mankind, feeds more than half the world population. Rice has been heralded as a model cereal owing to its small genome size, amenability to easy transformation, high synteny to other cereal crops and availability of complete genome sequence. Moreover, sequence wealth in rice is getting more refined and precise due to resequencing efforts. This humungous resource of sequence data has confronted research fraternity with a herculean challenge as well as an excellent opportunity to functionally validate expressed as well as regulatory portions of the genome. This will not only help us in understanding the genetic basis of plant architecture and physiology but would also steer us towards developing improved cultivars. No single technique can achieve such a mammoth task. Functional genomics through its diverse tools viz. loss and gain of function mutants, multifarious omics strategies like transcriptomics, proteomics, metabolomics and phenomics provide us with the necessary handle. A paradigm shift in technological advances in functional genomics strategies has been instrumental in generating considerable amount of information w.r.t functionality of rice genome. We now have several databases and online resources for functionally validated genes but despite that we are far from reaching the desired milestone of functionally characterizing each and every rice gene. There is an urgent need for a common platform, for information already available in rice, and collaborative efforts between researchers in a concerted manner as well as healthy public-private partnership, for genetic improvement of rice crop better able to handle the pressures of climate change and exponentially increasing population.

Pollen Killer Gene S35 Function Requires Interaction with an Activator That Maps Close to S24, Another Pollen Killer Gene in Rice

Pollen killer genes disable noncarrier pollens, and are responsible for male sterility and segregation distortion in hybrid populations of distantly related plant species. The genetic networks and the molecular mechanisms underlying the pollen killer system remain largely unknown. Two pollen killer genes, S24 and S35, have been found in an intersubspecific cross of Oryza sativa ssp. indica and japonica The effect of S24 is counteracted by an unlinked locus EFS Additionally, S35 has been proposed to interact with S24 to induce pollen sterility. These genetic interactions are suggestive of a single S24-centric genetic pathway (EFS-S24-S35) for the pollen killer system. To examine this hypothetical genetic pathway, the S35 and the S24 regions were further characterized and genetically dissected in this study. Our results indicated that S35 causes pollen sterility independently of both the EFS and S24 genes, but is dependent on a novel gene close to the S24 locus, named incentive for killing pollen (INK). We confirmed the phenotypic effect of the INK gene separately from the S24 gene, and identified the INK locus within an interval of less than 0.6 Mb on rice chromosome 5. This study characterized the genetic effect of the two independent genetic pathways of INK-S35 and EFS-S24 in indica-japonica hybrid progeny. Our results provide clear evidence that hybrid male sterility in rice is caused by several pollen killer networks with multiple factors positively and negatively regulating pollen killer genes.

Comparative Small RNA Analysis of Pollen Development in Autotetraploid and Diploid Rice

MicroRNAs (miRNAs) play key roles in plant reproduction. However, knowledge on microRNAome analysis in autotetraploid rice is rather limited. Here, high-throughput sequencing technology was employed to analyze miRNAomes during pollen development in diploid and polyploid rice. A total of 172 differentially expressed miRNAs (DEM) were detected in autotetraploid rice compared to its diploid counterpart, and 57 miRNAs were specifically expressed in autotetraploid rice. Of the 172 DEM, 115 and 61 miRNAs exhibited up- and down-regulation, respectively. Gene Ontology analysis on the targets of up-regulated DEM showed that they were enriched in transport and membrane in pre-meiotic interphase, reproduction in meiosis, and nucleotide binding in single microspore stage. osa-miR5788 and osa-miR1432-5p_R+1 were up-regulated in meiosis and their targets revealed interaction with the meiosis-related genes, suggesting that they may involve in the genes regulation associated with the chromosome behavior. Abundant 24 nt siRNAs associated with transposable elements were found in autotetraploid rice during pollen development; however, they significantly declined in diploid rice, suggesting that 24 nt siRNAs may play a role in pollen development. These findings provide a foundation for understanding the effect of polyploidy on small RNA expression patterns during pollen development that cause pollen sterility in autotetraploid rice.

Two Tightly Linked Genes at the hsa1 Locus Cause Both F1 and F2 Hybrid Sterility in Rice

Molecular mechanisms of hybrid breakdown associated with sterility (F2 sterility) are poorly understood as compared with those of F1 hybrid sterility. Previously, we characterized three unlinked epistatic loci, hybrid sterility-a1 (hsa1), hsa2, and hsa3, responsible for the F2 sterility in a cross between Oryza sativa ssp. indica and japonica. In this study, we identified that the hsa1 locus contains two interacting genes, HSA1a and HSA1b, within a 30-kb region. HSA1a-j (japonica allele) encodes a highly conserved plant-specific domain of unknown function protein (DUF1618), whereas the indica allele (HSA1a-i(s)) has two deletion mutations that cause disruption of domain structure. The second gene, HSA1b-i(s), encodes an uncharacterized protein with some similarity to a nucleotide-binding protein. Homozygous introgression of indica HSA1a-i(s)-HSA1b-i(s) alleles into japonica showed female gamete abortion at an early mitotic stage. The fact that the recombinant haplotype HSA1a-j-HSA1b-i(s) caused semi-sterility in the heterozygous state with the HSA1a-i(s)-HSA1b-i(s) haplotype suggests that variation in the hsa1 locus is a possible cause of the wide-spectrum sterility barriers seen in F1 hybrids and successive generations in rice. We propose a simple genetic model to explain how a single causal mechanism can drive both F1 and F2 hybrid sterility.

OryzaGenome: Genome Diversity Database of Wild Oryza Species

The species in the genus Oryza, encompassing nine genome types and 23 species, are a rich genetic resource and may have applications in deeper genomic analyses aiming to understand the evolution of plant genomes. With the advancement of next-generation sequencing (NGS) technology, a flood of Oryza species reference genomes and genomic variation information has become available in recent years. This genomic information, combined with the comprehensive phenotypic information that we are accumulating in our Oryzabase, can serve as an excellent genotype-phenotype association resource for analyzing rice functional and structural evolution, and the associated diversity of the Oryza genus. Here we integrate our previous and future phenotypic/habitat information and newly determined genotype information into a united repository, named OryzaGenome, providing the variant information with hyperlinks to Oryzabase. The current version of OryzaGenome includes genotype information of 446 O. rufipogon accessions derived by imputation and of 17 accessions derived by imputation-free deep sequencing. Two variant viewers are implemented: SNP Viewer as a conventional genome browser interface and Variant Table as a text-based browser for precise inspection of each variant one by one. Portable VCF (variant call format) file or tab-delimited file download is also available. Following these SNP (single nucleotide polymorphism) data, reference pseudomolecules/scaffolds/contigs and genome-wide variation information for almost all of the closely and distantly related wild Oryza species from the NIG Wild Rice Collection will be available in future releases. All of the resources can be accessed through http://viewer.shigen.info/oryzagenome/.