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

Spatiotemporal profiles of gene activity in stamen delineate nucleo-cytoplasmic interaction in a male-sterile somatic cybrid citrus

Cytoplasmic male sterility (CMS) has long been used to produce seedless fruits in perennial woody crops like citrus. A male-sterile somatic cybrid citrus (G1 + HBP) was generated by protoplast fusion between a CMS callus parent 'Guoqing No. 1' Satsuma mandarin (Citrus unshiu, G1) and a fertile mesophyll parent Hirado Buntan pummelo (Citrus grandis, HBP). To uncover the male-sterile mechanism of G1 + HBP, we compared the transcriptome profiles of stamen organ and cell types at five stages between G1 + HBP and HBP, including the initial stamen primordia, enlarged stamen primordia, pollen mother cells, tetrads, and microspores captured by laser microdissection. The stamen organ and cell types showed distinct gene expression profiles. A majority of genes involved in stamen development were differentially expressed, especially CgAP3.2, which was downregulated in enlarged stamen primordia and upregulated in tetrads of G1 + HBP compared with HBP. Jasmonic acid- and auxin-related biological processes were enriched among the differentially expressed genes of stamen primordia, and the content of jasmonic acid biosynthesis metabolites was higher in flower buds and anthers of G1 + HBP. In contrast, the content of auxin biosynthesis metabolites was lower in G1 + HBP. The mitochondrial tricarboxylic acid cycle and oxidative phosphorylation processes were enriched among the differentially expressed genes in stamen primordia, meiocytes, and microspores, indicating the dysfunction of mitochondria in stamen organ and cell types of G1 + HBP. Taken together, the results indicate that malfunction of mitochondria-nuclear interaction might cause disorder in stamen development, and thus lead to male sterility in the citrus cybrid.

Characterization and analysis of the promoter region of monodehydroascorbate reductase 4 (CpMDAR4) in papaya

Differential spatial and temporal expression patterns due to regulatory cis-elements and two different isoforms are detected among CpMDAR4 alleles in papaya. The aim of this research was to study the effects of cis-element differences between the X, Y and Y^h alleles on the expression of CpMDAR4, a potential candidate gene for sex differentiation in papaya, using a transcriptional reporter system in a model species Arabidopsis thaliana. Possible effects of a retrotransposon insertion in the Y and Y^h alleles on the transcription and expression of CpMDAR4 alleles in papaya flowers were also examined. When comparing promoters and cis-regulatory elements among genes in the non-recombining region of the sex chromosomes, paired genes exhibited differences. Our results showed that differences in the promoter sequences of the CpMDAR4 alleles drove the expression of a reporter gene to different flower tissues in Arabidopsis. β-glucuronidase staining analysis of T₂ and T₃ lines for constructs containing 5' deletions of native Y and Y^h allele promoters showed the loss of specific expression of the reporter gene in the anthers, confirming the existence and location of cis-regulatory element POLLEN1LELAT52. The expression analysis of CpMDAR4 alleles in papaya flowers also showed that all alleles are actively expressed in different flower tissues, with the existence of a shorter truncated isoform, with unknown function, for the Y and Y^h alleles due to an LTR-RT insertion in the Y and Y^h chromosomes. The observed expression patterns in Arabidopsis thaliana flowers and the expression patterns of CpMDAR4 alleles in papaya flowers suggest that MDAR4 might have a role on development of reproductive organs in papaya, and that it constitutes an important candidate for sex differentiation.

Rapid Identification of Pollen- and Anther-Specific Genes in Response to High-Temperature Stress Based on Transcriptome Profiling Analysis in Cotton

Anther indehiscence and pollen sterility caused by high temperature (HT) stress have become a major problem that decreases the yield of cotton. Pollen- and anther-specific genes play a critical role in the process of male reproduction and the response to HT stress. In order to identify pollen-specific genes that respond to HT stress, a comparative transcriptome profiling analysis was performed in the pollen and anthers of Gossypium hirsutum HT-sensitive Line H05 against other tissue types under normal temperature (NT) conditions, and the analysis of a differentially expressed gene was conducted in the pollen of H05 under NT and HT conditions. In total, we identified 1111 pollen-specific genes (PSGs), 1066 anther-specific genes (ASGs), and 833 pollen differentially expressed genes (DEGs). Moreover, we found that the late stage of anther included more anther- and pollen-specific genes (APSGs). Stress-related cis-regulatory elements (CREs) and hormone-responsive CREs are enriched in the promoters of APSGs, suggesting that APSGs may respond to HT stress. However, 833 pollen DEGs had only 10 common genes with 1111 PSGs, indicating that PSGs are mainly involved in the processes of pollen development and do not respond to HT stress. Promoters of these 10 common genes are enriched for stress-related CREs and MeJA-responsive CREs, suggesting that these 10 common genes are involved in the process of pollen development while responding to HT stress. This study provides a pathway for rapidly identifying cotton pollen-specific genes that respond to HT stress.

Overexpression of the Panax ginseng CYP703 Alters Cutin Composition of Reproductive Tissues in Arabidopsis

Cytochrome P450 (CYP) catalyzes a wide variety of monooxygenation reactions in plant primary and secondary metabolisms. Land plants contain CYP703, belonging to the CYP71 clan, which catalyzes the biochemical pathway of fatty acid hydroxylation, especially in male reproductive tissues. Korean/Asian ginseng (Panax ginseng Meyer) has been regarded as one of important medicinal plant for a long time, however the molecular mechanism is less known on its development. In this study, we identified and characterized a CYP703A gene in P. ginseng (PgCYP703A4), regarding reproductive development. PgCYP703A4 shared a high-sequence identity (81–83%) with predicted amino acid as CYP703 in Dancus carota, Pistacia vera, and Camellia sinensis as well as 76% of amino acid sequence identity with reported CYP703 in Arabidopsis thaliana and 75% with Oryza sativa. Amino acid alignment and phylogenetic comparison of P. ginseng with higher plants and known A. thaliana members clearly distinguish the CYP703 members, each containing the AATDTS oxygen binding motif and PERH as a clade signature. The expression of PgCYP704B1 was only detected in P. ginseng flower buds, particularly in meiotic cells and the tapetum layer of developing anther, indicating the conserved role on male reproduction with At- and Os- CYP703. To acquire the clue of function, we transformed the PgCYP703A4 in A. thaliana. Independent overexpressing lines (PgCYP703A4ox) increased silique size and seed number, and altered the contents of fatty acids composition of cutin monomer in the siliques. Our results indicate that PgCYP703A4 is involved in fatty acid hydroxylation which affects cutin production and fruit size.

SDR7-6, a short-chain alcohol dehydrogenase/reductase family protein, regulates light-dependent cell death and defence responses in rice

Lesion mimic mutants resembling the hypersensitive response without pathogen attack are an ideal material to understand programmed cell death, the defence response, and the cross-talk between defence response and development in plants. In this study, mic, a lesion mimic mutant from cultivar Yunyin treated with ethyl methanesulphonate (EMS), was screened. By map-based cloning, a short-chain alcohol dehydrogenase/reductase with an atypical active site HxxxK was isolated and designated as SDR7-6. It functions as a homomultimer in rice and is localized at the endoplasmic reticulum. The lesion mimic phenotype of the mutant is light-dependent. The mutant displayed an increased resistance response to bacterial blight, but reduced resistance to rice blast disease. The mutant and knockout lines showed increased reactive oxygen species, jasmonic acid content, antioxidant enzyme activity, and expression of pathogenicity-related genes, while chlorophyll content was significantly reduced. The knockout lines showed significant reduction in grain size, seed setting rate, 1000-grain weight, grain weight per plant, panicle length, and plant height. SDR7-6 is a new lesion mimic gene that encodes a short-chain alcohol dehydrogenase with atypical catalytic site. Disruption of SDR7-6 led to cell death and had adverse effects on multiple agricultural characters. SDR7-6 may act at the interface of the two defence pathways of bacterial blight and rice blast disease in rice.

A Rapid Pipeline for Pollen- and Anther-Specific Gene Discovery Based on Transcriptome Profiling Analysis of Maize Tissues

Recently, crop breeders have widely adopted a new biotechnology-based process, termed Seed Production Technology (SPT), to produce hybrid varieties. The SPT does not produce nuclear male-sterile lines, and instead utilizes transgenic SPT maintainer lines to pollinate male-sterile plants for propagation of nuclear-recessive male-sterile lines. A late-stage pollen-specific promoter is an essential component of the pollen-inactivating cassette used by the SPT maintainers. While a number of plant pollen-specific promoters have been reported so far, their usefulness in SPT has remained limited. To increase the repertoire of pollen-specific promoters for the maize community, we conducted a comprehensive comparative analysis of transcriptome profiles of mature pollen and mature anthers against other tissue types. We found that maize pollen has much less expressed genes (>1 FPKM) than other tissue types, but the pollen grain has a large set of distinct genes, called pollen-specific genes, which are exclusively or much higher (100 folds) expressed in pollen than other tissue types. Utilizing transcript abundance and correlation coefficient analysis, 1215 mature pollen-specific (MPS) genes and 1009 mature anther-specific (MAS) genes were identified in B73 transcriptome. These two gene sets had similar GO term and KEGG pathway enrichment patterns, indicating that their members share similar functions in the maize reproductive process. Of the genes, 623 were shared between the two sets, called mature anther- and pollen-specific (MAPS) genes, which represent the late-stage pollen-specific genes of the maize genome. Functional annotation analysis of MAPS showed that 447 MAPS genes (71.7% of MAPS) belonged to genes encoding pollen allergen protein. Their 2-kb promoters were analyzed for cis-element enrichment and six well-known pollen-specific cis-elements (AGAAA, TCCACCA, TGTGGTT, [TA]AAAG, AAATGA, and TTTCT) were found highly enriched in the promoters of MAPS. Interestingly, JA-responsive cis-element GCC box (GCCGCC) and ABA-responsive cis-element-coupling element1 (ABRE-CE1, CCACC) were also found enriched in the MAPS promoters, indicating that JA and ABA signaling likely regulate pollen-specific MAPS expression. This study describes a robust and straightforward pipeline to discover pollen-specific promotes from publicly available data while providing maize breeders and the maize industry a number of late-stage (mature) pollen-specific promoters for use in SPT for hybrid breeding and seed production.

Diverse Roles of MAX1 Homologues in Rice

Cytochrome P450 enzymes encoded by MORE AXILLARY GROWTH1 (MAX1)-like genes produce most of the structural diversity of strigolactones during the final steps of strigolactone biosynthesis. The diverse copies of MAX1 in Oryza sativa provide a resource to investigate why plants produce such a wide range of strigolactones. Here we performed in silico analyses of transcription factors and microRNAs that may regulate each rice MAX1, and compared the results with available data about MAX1 expression profiles and genes co-expressed with MAX1 genes. Data suggest that distinct mechanisms regulate the expression of each MAX1. Moreover, there may be novel functions for MAX1 homologues, such as the regulation of flower development or responses to heavy metals. In addition, individual MAX1s could be involved in specific functions, such as the regulation of seed development or wax synthesis in rice. Our analysis reveals potential new avenues of strigolactone research that may otherwise not be obvious.

Ubiquitinome Profiling Reveals the Landscape of Ubiquitination Regulation in Rice Young Panicles

Ubiquitination, an essential post-transcriptional modification (PTM), plays a vital role in nearly every biological process, including development and growth. Despite its functions in plant reproductive development, its targets in rice panicles remain unclear. In this study, we used proteome-wide profiling of lysine ubiquitination in rice (O. sativa ssp. indica) young panicles. We created the largest ubiquitinome dataset in rice to date, identifying 1638 lysine ubiquitination sites on 916 unique proteins. We detected three conserved ubiquitination motifs, noting that acidic glutamic acid (E) and aspartic acid (D) were most frequently present around ubiquitinated lysine. Enrichment analysis of Gene Ontology (GO) annotations and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of these ubiquitinated proteins revealed that ubiquitination plays an important role in fundamental cellular processes in rice young panicles. Interestingly, enrichment analysis of protein domains indicated that ubiquitination was enriched on a variety of receptor-like kinases and cytoplasmic tyrosine and serine-threonine kinases. Furthermore, we analyzed the crosstalk between ubiquitination, acetylation, and succinylation, and constructed a potential protein interaction network within our rice ubiquitinome. Moreover, we identified ubiquitinated proteins related to pollen and grain development, indicating that ubiquitination may play a critical role in the physiological functions in young panicles. Taken together, we reported the most comprehensive lysine ubiquitinome in rice so far, and used it to reveal the functional role of lysine ubiquitination in rice young panicles.

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.

Heavy Metal Stress-Associated Proteins in Rice and Arabidopsis: Genome-Wide Identification, Phylogenetics, Duplication, and Expression Profiles Analysis

Heavy metal exposure is a serious environmental stress in plants. However, plants have evolved several strategies to improve their heavy metal tolerance. Heavy metal-associated proteins (HMPs) participate in heavy metal detoxification. Here, we identified 46 and 55 HMPs in rice and Arabidopsis, respectively, and named them OsHMP 1–46 and AtHMP 1–55 according to their chromosomal locations. The HMPs from both plants were divided into six clades based on the characteristics of their heavy metal-associated domains (HMA). The HMP gene structures and motifs varied greatly among the different classifications. The HMPs had high collinearity and were segmentally duplicated. A cis-element analysis revealed that the HMPs may be regulated by different transcription factors. An expression profile analysis disclosed that only eight OsHMPs were constitutive in rice tissues. Of these, the expression of OsHMP37 was far higher than that of the other seven genes while OsHMP28 was expressed exclusively in the roots. For Arabidopsis, nine AtHMPs presented with very high transcript levels in all organs. Most of the selected OsHMPs were differentially expressed in various tissues under different heavy metal stresses. Only OsHMP09, OsHMP18, and OsHMP22 showed higher expression levels in all tissues under different heavy metal stresses. In contrast, most of the selected AtHMPs had nearly constant expression levels in different tissues under various heavy metal stresses. The AtHMP20, AtHMP23, AtHMP25, AtHMP31, AtHMP35, AtHMP46 expression levels under different heavy metal stresses were higher in the leaves and roots. The foregoing discoveries elucidated HMP evolution in monocotyledonous and dicotyledonous plants and may helpful functionally characterize HMPs in the future.

In Silico Identification of QTL-Based Polymorphic Genes as Salt-Responsive Potential Candidates through Mapping with Two Reference Genomes in Rice

Recent advances in next generation sequencing have created opportunities to directly identify genetic loci and candidate genes for abiotic stress responses in plants. With the objective of identifying candidate genes within the previously identified QTL-hotspots, the whole genomes of two divergent cultivars for salt responses, namely At 354 and Bg 352, were re-sequenced using Illumina Hiseq 2500 100PE platform and mapped to Nipponbare and R498 genomes. The sequencing results revealed approximately 2.4 million SNPs and 0.2 million InDels with reference to Nipponbare while 1.3 million and 0.07 million with reference to R498 in two parents. In total, 32,914 genes were reported across all rice chromosomes of this study. Gene mining within QTL hotspots revealed 1236 genes, out of which 106 genes were related to abiotic stress. In addition, 27 abiotic stress-related genes were identified in non-QTL regions. Altogether, 32 genes were identified as potential genes containing polymorphic non-synonymous SNPs or InDels between two parents. Out of 10 genes detected with InDels, tolerant haplotypes of Os01g0581400, Os10g0107000, Os11g0655900, Os12g0622500, and Os12g0624200 were found in the known salinity tolerant donor varieties. Our findings on different haplotypes would be useful in developing resilient rice varieties for abiotic stress by haplotype-based breeding studies.

Deficiency of rice hexokinase HXK5 impairs synthesis and utilization of starch in pollen grains and causes male sterility

There is little known about the function of rice hexokinases (HXKs) in planta. We characterized hxk5-1, a Tos17 mutant of OsHXK5 that is up-regulated in maturing pollen, a stage when starch accumulates. Progeny analysis of self-pollinated heterozygotes of hxk5-1 and reciprocal crosses between the wild-type and heterozygotes revealed that loss of HXK5 causes male sterility. Homozygous hxk5-1, produced via anther culture, and additional homozygous hxk5-2, hxk5-3 and hxk5-4 lines created by CRISPR/Cas9 confirmed the male-sterile phenotype. In vitro pollen germination ability and in vivo pollen tube growth rate were significantly reduced in the hxk5 mutant pollen. Biochemical analysis of anthers with the mutant pollen revealed significantly reduced hexokinase activity and starch content, although they were sufficient to produce some viable seed. However, the mutant pollen was unable to compete successfully against wild-type pollen. Expression of the catalytically inactive OsHXK5-G113D did not rescue the hxk5 male-sterile phenotype, indicating that its catalytic function was responsible for pollen fertility, rather than its role in sugar sensing and signaling. Our results demonstrate that OsHXK5 contributes to a large portion of the hexokinase activity necessary for the starch utilization pathway during pollen germination and tube growth, as well as for starch biosynthesis during pollen maturation.

Genome-Wide analysis of aluminum-activated malate transporter family genes in six rosaceae species, and expression analysis and functional characterization on malate accumulation in Chinese white pear

Aluminum-activated malate transporters (ALMTs) exhibit a variety of physiological roles in plants to regulate fruit quality, but the evolutionary history of the ALMT family in the Rosaceae species remains unknown. In this study, a total of 113 ALMT homologous genes were identified from six Rosaceae species (Pyrus bretschneideri, Malus × domestica, Prunus persica, Fragaria vesca, Prunus mume, and Pyrus communis), and 27 of these sequences came from Chinese white pear, designated PbrALMT. Based on the phylogenetic analysis, we divided these ALMT genes into three main clusters (A-C). Conserved domain analysis indicated that all PbrALMT proteins contained the ALMT domain and the FUSC_2 domain, and fewer proteins included the FUSC domain. The results of subcellular localization experiments showed that parts of PbrALMT proteins containing the FUSC domain were located in the membrane. Collinearity analysis revealed that segmental and dispersed duplications were the primary forces underlying ALMT gene family expansion in the Rosaceae. Calculation of Ka/Ks between the paralogous pairs indicated that all of the genes in the PbrALMT family have evolved under negative selection. Combining the changes of malate content and transcriptome data analysis, five genes belonging to Cluster B were chosen for qRT-PCR, and the results revealed that Pbr020270.1, as a candidate gene, may play important roles in malate accumulation during pear fruit development. Further transgenic assay confirmed the above conclusion. The present study provides a foundation to better understand the molecular evolution of ALMT genes in pear and the functional characterization of these genes in the future.

Crosstalk between diurnal rhythm and water stress reveals an altered primary carbon flux into soluble sugars in drought-treated rice leaves

Plants retain rhythmic physiological responses when adapting to environmental challenges. However, possible integrations between drought conditions and those responses have not received much focus, especially regarding crop plants, and the relationship between abiotic stress and the diurnal cycle is generally not considered. Therefore, we conducted a genome-wide analysis to identify genes showing both diurnal regulation and water-deficiency response in rice (Oryza sativa). Among the 712 drought-responsive genes primary identified, 56.6% are diurnally expressed while 47.6% of the 761 that are down-regulated by drought are also diurnal. Using the β-glucuronidase reporter system and qRT-PCR analyses, we validated expression patterns of two candidate genes, thereby supporting the reliability of our transcriptome data. MapMan analysis indicated that diurnal genes up-regulated by drought are closely associated with the starch-sucrose pathway while those that are down-regulated are involved in photosynthesis. We then confirmed that starch-sucrose contents and chlorophyll fluorescence are altered in a diurnal manner under drought stress, suggesting these metabolic diurnal alterations as a novel indicator to evaluate the drought response in rice leaves. We constructed a functional gene network associated with the starch-sucrose KEGG metabolic pathway for further functional studies, and also developed a regulatory pathway model that includes OsbZIP23 transcription factor.

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.