Comprehensive network analysis of anther-expressed genes in rice by the combination of 33 laser microdissection and 143 spatiotemporal microarrays

Volatile metabolome and floral transcriptome analyses reveal the volatile components of strongly fragrant progeny of Malus × robusta

Floral fragrance is an important trait that contributes to the ornamental properties and pollination of crabapple. However, research on the physiological and molecular biology of the floral volatile compounds of crabapple is rarely reported. In this study, metabolomic and transcriptomic analyses of the floral volatile compounds of standard Malus robusta flowers (Mr), and progeny with strongly and weakly fragrant flowers (SF and WF, respectively), were conducted. Fifty-six floral volatile compounds were detected in the plant materials, mainly comprising phenylpropane/benzene ring-type compounds, fatty acid derivatives, and terpene compounds. The volatile contents were significantly increased before the early flowering stage (ES), and the contents of SF flowers were twice those of WF and Mr flowers. Odor activity values were determined for known fragrant volatiles and 10-11 key fragrant volatiles were identified at the ES. The predominant fragrant volatiles were methyl benzoate, linalool, leaf acetate, and methyl anthranilate. In the petals, stamens, pistil, and calyx of SF flowers, 26 volatiles were detected at the ES, among which phenylpropane/benzene ring-type compounds were the main components accounting for more than 75% of the total volatile content. Functional analysis of transcriptome data revealed that the phenylpropanoid biosynthesis pathway was significantly enriched in SF flowers. By conducting combined analyses between volatiles and differentially expressed genes, transcripts of six floral scent-related genes were identified and were associated with the contents of the key fragrant volatiles, and other 23 genes were potentially correlated with the key volatile compounds. The results reveal possible mechanisms for the emission of strong fragrance by SF flowers, and provide a foundation for improvement of the floral fragrance and development of new crabapple cultivars.

Identification of major candidate genes for multiple abiotic stress tolerance at seedling stage by network analysis and their validation by expression profiling in rice (Oryza sativa L.)

A wealth of microarray and RNA-seq data for studying abiotic stress tolerance in rice exists but only limited studies have been carried out on multiple stress-tolerance responses and mechanisms. In this study, we identified 6657 abiotic stress-responsive genes pertaining to drought, salinity and heat stresses from the seedling stage microarray data of 83 samples and used them to perform unweighted network analysis and to identify key hub genes or master regulators for multiple abiotic stress tolerance. Of the total 55 modules identified from the analysis, the top 10 modules with 8-61 nodes comprised 239 genes. From these 10 modules, 10 genes common to all the three stresses were selected. Further, based on the centrality properties and highly dense interactions, we identified 7 intra-modular hub genes leading to a total of 17 potential candidate genes. Out of these 17 genes, 15 were validated by expression analysis using a panel of 4 test genotypes and a pair of standard check genotypes for each abiotic stress response. Interestingly, all the 15 genes showed upregulation under all stresses and in all the genotypes, suggesting that they could be representing some of the core abiotic stress-responsive genes. More pertinently, eight of the genes were found to be co-localized with the stress-tolerance QTL regions. Thus, in conclusion, our study not only provided an effective approach for studying abiotic stress tolerance in rice, but also identified major candidate genes which could be further validated by functional genomics for abiotic stress tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03182-7.

Potential biomarkers of acute myocardial infarction based on co‑expression network analysis

Acute myocardial infarction (AMI) is a common cause of death in numerous countries. Understanding the molecular mechanisms of the disease and analyzing potential biomarkers of AMI is crucial. However, specific diagnostic biomarkers have thus far not been fully established and candidate regulatory targets for AMI remain to be determined. In the present study, the AMI gene chip dataset GSE48060 comprising blood samples from control subjects with normal cardiac function (n=21) and patients with AMI (n=26) was downloaded from Gene Expression Omnibus. The differentially expressed genes (DEGs) between the AMI and control groups were identified with the online tool GEO2R. The co-expression network of DEGs was analyzed by calculating the Pearson correlation coefficient of all gene pairs, mutual rank screening and cutoff threshold screening. Subsequently, the Gene Ontology (GO) database was used to analyze the genes' functions and pathway enrichment of genes in the most important modules was performed. Kyoto Encyclopedia of Genes and Genomes (KEGG) Disease and BioCyc were used to analyze the hub genes in the module to determine important sub-pathways. In addition, the expression of hub genes was confirmed by reverse transcription-quantitative PCR in AMI and control specimens. In the present study, 52 DEGs, including 26 upregulated and 26 downregulated genes, were identified. As key hub genes, three upregulated genes (AKR1C3, RPS24 and P2RY12) and three downregulated genes (ACSL1, B3GNT5 and MGAM) were identified from the co-expression network. Furthermore, GO enrichment analysis of all AMI co-expression network genes revealed functional enrichment mainly in ‘RAGE receptor binding’ and ‘negative regulation of T cell cytokine production’. In addition, KEGG Disease and BioCyc analysis indicated functional enrichment of the genes RPS24 and P2RY12 in ‘cardiovascular diseases’, of AKR1C3 in ‘cardenolide biosynthesis’, of MGAM in ‘glycogenolysis’, of B3GNT5 in ‘glycosphingolipid biosynthesis’ and of ACSL1 in ‘icosapentaenoate biosynthesis II’. In conclusion, the hub genes AKR1C3, RPS24, P2RY12, ACSL1, B3GNT5 and MGAM are potential markers of AMI, and have potential application value in the diagnosis of AMI.

QTL mapping and candidate gene mining of flag leaf size traits in Japonica rice based on linkage mapping and genome-wide association study

BACKGROUND

As one of the most important factors of the japonica rice plant, leaf shape affects the photosynthesis and carbohydrate accumulation directly. Mining and using new leaf shape related genes/QTLs can further enrich the theory of molecular breeding and accelerate the breeding process of japonica rice.

METHODS

In the present study, 2 RILs and a natural population with 295 japonica rice varieties were used to map QTLs for flag leaf length (FL), flag leaf width (FW) and flag leaf area (FLA) by linkage analysis and genome-wide association study (GWAS) throughout 2 years.

RESULTS

A total of 64 QTLs were detected by 2 ways, and pleiotropic QTLs qFL2 (Chr2_33,332,579) and qFL10 (Chr10_10,107,835; Chr10_10,230,100) consisted of overlapping QTLs mapped by linkage analysis and GWAS throughout the 2 years were identified.

CONCLUSIONS

The candidate genes LOC_Os02g54254, LOC_Os02g54550, LOC_Os10g20160, LOC_Os10g20240, LOC_Os10g20260 were obtained, filtered by linkage disequilibrium (LD), and haplotype analysis. LOC_Os10g20160 (SD-RLK-45) showed outstanding characteristics in quantitative real-time PCR (qRT-PCR) analysis in leaf development period, belongs to S-domain receptor-like protein kinases gene and probably to be a main gene regulating flag leaf width of japonica rice. The results of this study provide valuable resources for mining the main genes/QTLs of japonica rice leaf development and molecular breeding of japonica rice ideal leaf shape.

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.

Rice transcription factor GAMYB modulates bHLH142 and is homeostatically regulated by TDR during anther tapetal and pollen development

GIBBERELLIN MYB GENE (GAMYB), UNDEVELOPED TAPETUM1 (UDT1), TDR INTERACTING PROTEIN2 (TIP2/bHLH142), TAPETUM DEGENERATION RETARDATION (TDR), and ETERNAL TAPETUM 1/DELAYED TAPETUM DEGENERATION (EAT1/DTD) are important transcription factors that play a crucial role during pollen development in rice. This study demonstrates that bHLH142 acts downstream of UDT1 and GAMYB and works as a 'hub' in these two pollen pathways. We show that GAMYB modulates bHLH142 expression through specific binding to the MYB motif of the bHLH142 promoter during the early stage of pollen development, while TDR acts as a transcriptional repressor of the GAMYB modulation of bHLH142 by binding to the E-box close to the MYB motif on the promoter. Altered expression of these transcription factors highlights that a tight, precise, and coordinated regulation among them is essential for normal pollen development. Most notably, we show that the regulatory pathways of GAMYB and UDT1 rely on bHLH142 in a direct and indirect manner, respectively, and function in different tissues with distinct biological roles during pollen development. This study advances our understanding of the molecular mechanisms of rice pollen development.

Transcriptome analyses and weighted gene coexpression network analysis reveal key pathways and genes involved in the rapid cold resistance of the Chinese white wax scale insect

The Chinese white wax scale insect, Ericerus pela, is an important resource insect in China. The rapid response of E. pela to decreasing temperatures plays key roles in the population distribution. In this study, we analyzed the gene expression of E. pela treated with low temperature using transcriptome analyses and weighted gene coexpression network analysis (WGCNA). The results showed that the cold resistance of E. pela involved changes in the expression of many genes. The genes were mainly involved in alcohol formation activity, lipid metabolism, membrane and structure, and oxidoreductase activity. According to the WGCNA results, some pathways related to cold resistance were found in the genes in the modules, such as cytoskeleton proteins, cytoskeleton protein pathway, biosynthesis of unsaturated fatty acids, glycerophospholipid metabolism, ether lipid metabolism, and thermogenesis. Some of the hub genes were nonspecific lipid-transfer proteins, DnaJ homolog subfamily C member 13, paramyosin, tropomodulin, and tubulin beta chain. In particular, the hub genes of the tan module included the heat shock protein (hsp) 10, hsp 60, hsp 70, and hsp 90 genes. Thirty-five antifreeze protein (afp) genes were identified according to the annotation results. Three afp genes were further identified among the hub genes. Six of these genes were selected for heterogeneous protein expression. One of them was expressed successfully. The thermal hysteresis activity (THA) analyses showed that the THA was 1.73°C. These results showed that the cytoskeleton, lipid metabolism, thermogenesis, HSPs and AFPs may play important roles in the cold resistance of E. pela.

Beyond gene ontology (GO): using biocuration approach to improve the gene nomenclature and functional annotation of rice S-domain kinase subfamily

The S-domain subfamily of receptor-like kinases (SDRLKs) in plants is poorly characterized. Most members of this subfamily are currently assigned gene function based on the S-locus Receptor Kinase from Brassica that acts as the female determinant of self-incompatibility (SI). However, Brassica like SI mechanisms does not exist in most plants. Thus, automated Gene Ontology (GO) pipelines are not sufficient for functional annotation of SDRLK subfamily members and lead to erroneous association with the GO biological process of SI. Here, we show that manual bio-curation can help to correct and improve the gene annotations and association with relevant biological processes. Using publicly available genomic and transcriptome datasets, we conducted a detailed analysis of the expansion of the rice (Oryza sativa) SDRLK subfamily, the structure of individual genes and proteins, and their expression.The 144-member SDRLK family in rice consists of 82 receptor-like kinases (RLKs) (67 full-length, 15 truncated),12 receptor-like proteins, 14 SD kinases, 26 kinase-like and 10 GnK2 domain-containing kinases and RLKs. Except for nine genes, all other SDRLK family members are transcribed in rice, but they vary in their tissue-specific and stress-response expression profiles. Furthermore, 98 genes show differential expression under biotic stress and 98 genes show differential expression under abiotic stress conditions, but share 81 genes in common.Our analysis led to the identification of candidate genes likely to play important roles in plant development, pathogen resistance, and abiotic stress tolerance. We propose a nomenclature for 144 SDRLK gene family members based on gene/protein conserved structural features, gene expression profiles, and literature review. Our biocuration approach, rooted in the principles of findability, accessibility, interoperability and reusability, sets forth an example of how manual annotation of large-gene families can fill in the knowledge gap that exists due to the implementation of automated GO projections, thereby helping to improve the quality and contents of public databases.

OsMYB80 Regulates Anther Development and Pollen Fertility by Targeting Multiple Biological Pathways

Pollen development is critical to the reproductive success of flowering plants, but how it is regulated is not well understood. Here, we isolated two allelic male-sterile mutants of OsMYB80 and investigated how OsMYB80 regulates male fertility in rice. OsMYB80 was barely expressed in tissues other than anthers, where it initiated the expression during meiosis, reached the peak at the tetrad-releasing stage and then quickly declined afterward. The osmyb80 mutants exhibited premature tapetum cell death, lack of Ubisch bodies, no exine and microspore degeneration. To understand how OsMYB80 regulates anther development, RNA-seq analysis was conducted to identify genes differentially regulated by OsMYB80 in rice anthers. In addition, DNA affinity purification sequencing (DAP-seq) analysis was performed to identify DNA fragments interacting with OsMYB80 in vitro. Overlap of the genes identified by RNA-seq and DAP-seq revealed 188 genes that were differentially regulated by OsMYB80 and also carried an OsMYB80-interacting DNA element in the promoter. Ten of these promoter elements were randomly selected for gel shift assay and yeast one-hybrid assay, and all showed OsMYB80 binding. The 10 promoters also showed OsMYB80-dependent induction when co-expressed in rice protoplast. Functional annotation of the 188 genes suggested that OsMYB80 regulates male fertility by directly targeting multiple biological processes. The identification of these genes significantly enriched the gene networks governing anther development and provided much new information for the understanding of pollen development and male fertility.

Fertile Tetraploids: New Resources for Future Rice Breeding?

Ploidy manipulation is an efficient technique for the development of novel phenotypes in plant breeding. However, in rice (Oryza sativa L.), severe seed sterility has been considered a barrier preventing cultivation of autotetraploids since the 1930s. Recently, a series of studies identified two fertile autotetraploids, identified herein as the PMeS (Polyploid Meiosis Stability) and Neo-Tetraploid lines. Here, we summarize their characteristics, focusing on the recovery of seed fertility, and discuss potential future directions of study in this area, providing a comprehensive understanding of current progress in the study of fertile tetraploid rice, a classical, but promising, concept for rice breeding.

A Co-Expression Network in Hexaploid Wheat Reveals Mostly Balanced Expression and Lack of Significant Gene Loss of Homeologous Meiotic Genes Upon Polyploidization

Polyploidization has played an important role in plant evolution. However, upon polyploidization, the process of meiosis must adapt to ensure the proper segregation of increased numbers of chromosomes to produce balanced gametes. It has been suggested that meiotic gene (MG) duplicates return to a single copy following whole genome duplication to stabilize the polyploid genome. Therefore, upon the polyploidization of wheat, a hexaploid species with three related (homeologous) genomes, the stabilization process may have involved rapid changes in content and expression of MGs on homeologous chromosomes (homeologs). To examine this hypothesis, sets of candidate MGs were identified in wheat using co-expression network analysis and orthology informed approaches. In total, 130 RNA-Seq samples from a range of tissues including wheat meiotic anthers were used to define co-expressed modules of genes. Three modules were significantly correlated with meiotic tissue samples but not with other tissue types. These modules were enriched for GO terms related to cell cycle, DNA replication, and chromatin modification and contained orthologs of known MGs. Overall, 74.4% of genes within these meiosis-related modules had three homeologous copies which was similar to other tissue-related modules. Amongst wheat MGs identified by orthology, rather than co-expression, the majority (93.7%) were either retained in hexaploid wheat at the same number of copies (78.4%) or increased in copy number (15.3%) compared to ancestral wheat species. Furthermore, genes within meiosis-related modules showed more balanced expression levels between homeologs than genes in non-meiosis-related modules. Taken together, our results do not support extensive gene loss nor changes in homeolog expression of MGs upon wheat polyploidization. The construction of the MG co-expression network allowed identification of hub genes and provided key targets for future studies.

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.

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.

UDP-arabinopyranose mutase gene expressions are required for the biosynthesis of the arabinose side chain of both pectin and arabinoxyloglucan, and normal leaf expansion in Nicotiana tabacum

Plant cell walls are composed of polysaccharides such as cellulose, hemicelluloses, and pectins, whose location and function differ depending on plant type. Arabinose is a constituent of many different cell wall components, including pectic rhamnogalacturonan I (RG-I) and II (RG-II), glucuronoarabinoxylans (GAX), and arabinoxyloglucan (AXG). Arabinose is found predominantly in the furanose rather than in the thermodynamically more stable pyranose form. The UDP-arabinopyranose mutases (UAMs) have been demonstrated to convert UDP-arabinopyranose (UDP-Arap) to UDP-arabinofuranose (UDP-Araf) in rice (Oryza sativa L.). The UAMs have been implicated in polysaccharide biosynthesis and developmental processes. Arabinose residues could be a component of many polysaccharides, including branched (1→5)-α-arabinans, arabinogalactans in pectic polysaccharides, and arabinoxyloglucans, which are abundant in the cell walls of solanaceous plants. Therefore, to elucidate the role of UAMs and arabinan side chains, we analyzed the UAM RNA interference transformants in tobacco (Nicotiana tabacum L.). The tobacco UAM gene family consists of four members. We generated RNAi transformants (NtUAM-KD) to down-regulate all four of the UAM members. The NtUAM-KD showed abnormal leaf development in the form of a callus-like structure and many holes in the leaf epidermis. A clear reduction in the pectic arabinan content was observed in the tissue of the NtUAM-KD leaf. The arabinose/xylose ratio in the xyloglucan-rich cell wall fraction was drastically reduced in NtUAM-KD. These results suggest that UAMs are required for Ara side chain biosynthesis in both RG-I and AXG in Solanaceae plants, and that arabinan-mediated cell wall networks might be important for normal leaf expansion.

Phylogenomic detection and functional prediction of genes potentially important for plant meiosis

Meiosis is a specialized type of cell division necessary for sexual reproduction in eukaryotes. A better understanding of the cytological procedures of meiosis has been achieved by comprehensive cytogenetic studies in plants, while the genetic mechanisms regulating meiotic progression remain incompletely understood. The increasing accumulation of complete genome sequences and large-scale gene expression datasets has provided a powerful resource for phylogenomic inference and unsupervised identification of genes involved in plant meiosis. By integrating sequence homology and expression data, 164, 131, 124 and 162 genes potentially important for meiosis were identified in the genomes of Arabidopsis thaliana, Oryza sativa, Selaginella moellendorffii and Pogonatum aloides, respectively. The predicted genes were assigned to 45 meiotic GO terms, and their functions were related to different processes occurring during meiosis in various organisms. Most of the predicted meiotic genes underwent lineage-specific duplication events during plant evolution, with about 30% of the predicted genes retaining only a single copy in higher plant genomes. The results of this study provided clues to design experiments for better functional characterization of meiotic genes in plants, promoting the phylogenomic approach to the evolutionary dynamics of the plant meiotic machineries.

Affinity proteomics reveals extensive phosphorylation of the Brassica chromosome axis protein ASY1 and a network of associated proteins at prophase I of meiosis

During meiosis, the formation of crossovers (COs) generates genetic variation and provides physical links that are essential for accurate chromosome segregation. COs occur in the context of a proteinaceous chromosome axis. The transcriptomes and proteomes of anthers and meiocytes comprise several thousand genes and proteins, but because of the level of complexity relatively few have been functionally characterized. Our understanding of the physical and functional interactions between meiotic proteins is also limited. Here we use affinity proteomics to analyse the proteins that are associated with the meiotic chromosome axis protein, ASY1, in Brassica oleracea anthers and meiocytes. We show that during prophase I ASY1 and its interacting partner, ASY3, are extensively phosphorylated, and we precisely assign phosphorylation sites. We identify 589 proteins that co-immunoprecipitate with ASY1. These correspond to 492 Arabidopsis orthologues, over 90% of which form a coherent protein-protein interaction (PPI) network containing known and candidate meiotic proteins, including proteins more usually associated with other cellular processes such as DNA replication and proteolysis. Mutant analysis confirms that affinity proteomics is a viable strategy for revealing previously unknown meiotic proteins, and we show how the PPI network can be used to prioritise candidates for analysis. Finally, we identify another axis-associated protein with a role in meiotic recombination. Data are available via ProteomeXchange with identifier PXD006042.

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.

TGMS in Rapeseed (Brassica napus) Resulted in Aberrant Transcriptional Regulation, Asynchronous Microsporocyte Meiosis, Defective Tapetum, and Fused Sexine

The thermo-sensitive genic male sterility (TGMS) line SP2S is a spontaneous rapeseed mutation with several traits that are favorable for the production of two-line hybrids. To uncover the key cellular events and genetic regulation associated with TGMS expression, a combined study using cytological observation, transcriptome profiling, and gene expression analysis was conducted for SP2S and its near-isogenic line SP2F grown under warm conditions. Asynchronous microsporocyte meiosis and abnormal tapetal plastids and elaioplasts were demonstrated in the anther of SP2S. The tetrad microspore did not undergo mitosis before the cytoplasm degenerated. Delayed degradation of the tetrad wall, which led to tetrad microspore aggregation, resulted in postponement of sexine (outer layer of pollen exine) formation and sexine fusion in the tetrad. The nexine (foot layer of exine) was also absent. The delay of tetrad wall degradation and abnormality of the exine structure suggested that the defective tapetum lost important functions. Based on transcriptomic comparisons between young flower buds of SP2S and SP2F plants, a total of 465 differentially expressed transcripts (DETs) were identified, including 303 up-regulated DETs and 162 down-regulated DETs in SP2S. Several genes encoding small RNA degrading nuclease 2, small RNA 2'-O-methyltransferase, thioredoxin reductase 2, regulatory subunit A alpha isoform of serine/threonine-protein phosphatase 2A, glycine rich protein 1A, transcription factor bHLH25, leucine-rich repeat receptor kinase At3g14840 like, and fasciclin-like arabinogalactan proteins FLA19 and FLA20 were greatly depressed in SP2S. Interestingly, a POLLENLESS3-LIKE 2 gene encoding the Arabidopsis MS5 homologous protein, which is necessary for microsporocyte meiosis, was down-regulated in SP2S. Other genes that were up-regulated in SP2S encoded glucanase A6, ethylene-responsive transcription factor 1A-like, pollen-specific SF3, stress-associated endoplasmic reticulum protein 2, WRKY transcription factors and pentatricopeptide repeat (PPR) protein At1g07590. The tapetum-development-related genes, including BnEMS1, BnDYT1, and BnAMS, were slightly up-regulated in 3-mm-long flower buds or their anthers, and their downstream genes, BnMS1 and BnMYB80, which affect callose dissolution and exine formation, were greatly up-regulated in SP2S. This aberrant genetic regulation corresponded well with the cytological abnormalities. The results suggested that expression of TGMS associates with complex transcriptional regulation.

Tightly Controlled Expression of bHLH142 Is Essential for Timely Tapetal Programmed Cell Death and Pollen Development in Rice

Male sterility is important for hybrid seed production. Pollen development is regulated by a complex network. We previously showed that knockout of bHLH142 in rice (Oryza sativa) causes pollen sterility by interrupting tapetal programmed cell death (PCD) and bHLH142 coordinates with TDR to modulate the expression of EAT1. In this study, we demonstrated that overexpression of bHLH142 (OE142) under the control of the ubiquitin promoter also leads to male sterility in rice by triggering the premature onset of PCD. Protein of bHLH142 was found to accumulate specifically in the OE142 anthers. Overexpression of bHLH142 induced early expression of several key regulatory transcription factors in pollen development. In particular, the upregulation of EAT1 at the early stage of pollen development promoted premature PCD in the OE142 anthers, while its downregulation at the late stage impaired pollen development by suppressing genes involved in pollen wall biosynthesis, ROS scavenging and PCD. Collectively, these events led to male sterility in OE142. Analyses of related mutants further revealed the hierarchy of the pollen development regulatory gene network. Thus, the findings of this study advance our understanding of the central role played by bHLH142 in the regulatory network leading to pollen development in rice and how overexpression of its expression affects pollen development. Exploitation of this novel functionality of bHLH142 may confer a big advantage to hybrid seed production.

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.

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.

Gibberellic Acid: A Key Phytohormone for Spikelet Fertility in Rice Grain Production

The phytohormone gibberellic acid (GA) has essential signaling functions in multiple processes during plant development. In the "Green Revolution", breeders developed high-yield rice cultivars that exhibited both semi-dwarfism and altered GA responses, thus improving grain production. Most studies of GA have concentrated on germination and cell elongation, but GA also has a pivotal role in floral organ development, particularly in stamen/anther formation. In rice, GA signaling plays an important role in spikelet fertility; however, the molecular genetic and biochemical mechanisms of GA in male fertility remain largely unknown. Here, we review recent progress in understanding the network of GA signaling and its connection with spikelet fertility, which is tightly associated with grain productivity in cereal crops.

Cytoplasmic male sterility (CMS) in hybrid breeding in field crops

A comprehensive understanding of CMS/Rf system enabled by modern omics tools and technologies considerably improves our ability to harness hybrid technology for enhancing the productivity of field crops. Harnessing hybrid vigor or heterosis is a promising approach to tackle the current challenge of sustaining enhanced yield gains of field crops. In the context, cytoplasmic male sterility (CMS) owing to its heritable nature to manifest non-functional male gametophyte remains a cost-effective system to promote efficient hybrid seed production. The phenomenon of CMS stems from a complex interplay between maternally-inherited (mitochondrion) and bi-parental (nucleus) genomic elements. In recent years, attempts aimed to comprehend the sterility-inducing factors (orfs) and corresponding fertility determinants (Rf) in plants have greatly increased our access to candidate genomic segments and the cloned genes. To this end, novel insights obtained by applying state-of-the-art omics platforms have substantially enriched our understanding of cytoplasmic-nuclear communication. Concomitantly, molecular tools including DNA markers have been implicated in crop hybrid breeding in order to greatly expedite the progress. Here, we review the status of diverse sterility-inducing cytoplasms and associated Rf factors reported across different field crops along with exploring opportunities for integrating modern omics tools with CMS-based hybrid breeding.

Pre-Meiotic Anther Development: Cell Fate Specification and Differentiation

Research into anther ontogeny has been an active and developing field, transitioning from a strictly lineage-based view of cellular differentiation events to a more complex understanding of cell fate specification. Here we describe the modern interpretation of pre-meiotic anther development, from the earliest cell specifications within the anther lobes through SPL/NZZ-, MSP1-, and MEL1-dependent pathways as well as the initial setup of the abaxial and adaxial axes and outgrowth of the anther lobes. We then continue with a look at the known information regarding further differentiation of the somatic layers of the anther (the epidermis, endothecium, middle layer, and tapetum), with an emphasis on male-sterile mutants identified as defective in somatic cell specification. We also describe the differences in developmental stages among species and use this information to discuss molecular studies that have analyzed transcriptome, proteome, and small-RNA information in the anther.

Regulatory Networks in Pollen Development under Cold Stress

Cold stress modifies anthers' metabolic pathways to induce pollen sterility. Cold-tolerant plants, unlike the susceptible ones, produce high proportion of viable pollen. Anthers in susceptible plants, when exposed to cold stress, increase abscisic acid (ABA) metabolism and reduce ABA catabolism. Increased ABA negatively regulates expression of tapetum cell wall bound invertase and monosaccharide transport genes resulting in distorted carbohydrate pool in anther. Cold-stress also reduces endogenous levels of the bioactive gibberellins (GAs), GA4 and GA7, in susceptible anthers by repression of the GA biosynthesis genes. Here, we discuss recent findings on mechanisms of cold susceptibility in anthers which determine pollen sterility. We also discuss differences in regulatory pathways between cold-stressed anthers of susceptible and tolerant plants that decide pollen sterility or viability.

The role of OsMSH4 in male and female gamete development in rice meiosis

Meiosis is essential for gametogenesis in sexual reproduction in rice (Oryza sativa L.). We identified a MutS-homolog (MSH) family gene OsMSH4 in a trisomic plant. Cytological analysis showed that developments of both pollen and embryo sacs in an Osmsh4 mutant were blocked due to defective chromosome pairing. Compared with the wild type, the Osmsh4 mutant displayed a significant ~21.9% reduction in chiasma frequency, which followed a Poisson distribution, suggesting that class I crossover formation in the mutant was impaired. Temporal and spatial expression pattern analyses showed that OsMSH4 was preferentially expressed in meiocytes during their meiosis, indicating a critical role in gametogenesis. Subcellular localization showed that OsMSH4-green fluorescent protein was predominantly located in the nucleus. OsMSH4 could interact with another MSH member (OsMSH5) through the N-terminus and C-terminus, respectively. Direct physical interaction between OsMSH5, OsRPA1a, OsRPA2b, OsRPA1c, and OsRPA2c was identified by yeast two-hybrid assays and further validated by pull-down assays. Our results supported the conclusion that the OsMSH4/5 heterodimer plays a key role in regulation of crossover formation during rice meiosis by interaction with the RPA complex.

Trait Specific Expression Profiling of Salt Stress Responsive Genes in Diverse Rice Genotypes as Determined by Modified Significance Analysis of Microarrays

Stress responsive gene expression is commonly profiled in a comparative manner involving different stress conditions or genotypes with contrasting reputation of tolerance/resistance. In contrast, this research exploited a wide natural variation in terms of taxonomy, origin and salt sensitivity in eight genotypes of rice to identify the trait specific patterns of gene expression under salt stress. Genome wide transcptomic responses were interrogated by the weighted continuous morpho-physiological trait responses using modified Significance Analysis of Microarrays. More number of genes was found to be differentially expressed under salt stressed compared to that of under unstressed conditions. Higher numbers of genes were observed to be differentially expressed for the traits shoot Na(+)/K(+), shoot Na(+), root K(+), biomass and shoot Cl(-), respectively. The results identified around 60 genes to be involved in Na(+), K(+), and anion homeostasis, transport, and transmembrane activity under stressed conditions. Gene Ontology (GO) enrichment analysis identified 1.36% (578 genes) of the entire transcriptome to be involved in the major molecular functions such as signal transduction (>150 genes), transcription factor (81 genes), and translation factor activity (62 genes) etc., under salt stress. Chromosomal mapping of the genes suggests that majority of the genes are located on chromosomes 1, 2, 3, 6, and 7. The gene network analysis showed that the transcription factors and translation initiation factors formed the major gene networks and are mostly active in nucleus, cytoplasm and mitochondria whereas the membrane and vesicle bound proteins formed a secondary network active in plasma membrane and vacuoles. The novel genes and the genes with unknown functions thus identified provide picture of a synergistic salinity response representing the potentially fundamental mechanisms that are active in the wide natural genetic background of rice and will be of greater use once their roles are functionally verified.

Learning from Co-expression Networks: Possibilities and Challenges

Plants are fascinating and complex organisms. A comprehensive understanding of the organization, function and evolution of plant genes is essential to disentangle important biological processes and to advance crop engineering and breeding strategies. The ultimate aim in deciphering complex biological processes is the discovery of causal genes and regulatory mechanisms controlling these processes. The recent surge of omics data has opened the door to a system-wide understanding of the flow of biological information underlying complex traits. However, dealing with the corresponding large data sets represents a challenging endeavor that calls for the development of powerful bioinformatics methods. A popular approach is the construction and analysis of gene networks. Such networks are often used for genome-wide representation of the complex functional organization of biological systems. Network based on similarity in gene expression are called (gene) co-expression networks. One of the major application of gene co-expression networks is the functional annotation of unknown genes. Constructing co-expression networks is generally straightforward. In contrast, the resulting network of connected genes can become very complex, which limits its biological interpretation. Several strategies can be employed to enhance the interpretation of the networks. A strategy in coherence with the biological question addressed needs to be established to infer reliable networks. Additional benefits can be gained from network-based strategies using prior knowledge and data integration to further enhance the elucidation of gene regulatory relationships. As a result, biological networks provide many more applications beyond the simple visualization of co-expressed genes. In this study we review the different approaches for co-expression network inference in plants. We analyse integrative genomics strategies used in recent studies that successfully identified candidate genes taking advantage of gene co-expression networks. Additionally, we discuss promising bioinformatics approaches that predict networks for specific purposes.

Polyploidy Enhances F1 Pollen Sterility Loci Interactions That Increase Meiosis Abnormalities and Pollen Sterility in Autotetraploid Rice

Intersubspecific autotetraploid rice (Oryza sativa ssp. indica × japonica) hybrids have greater biological and yield potentials than diploid rice. However, the low fertility of intersubspecific autotetraploid hybrids, which is largely caused by high pollen abortion rates, limits their commercial utility. To decipher the cytological and molecular mechanisms underlying allelic interactions in autotetraploid rice, we developed an autotetraploid rice hybrid that was heterozygous (S(i)S(j)) at F1 pollen sterility loci (Sa, Sb, and Sc) using near-isogenic lines. Cytological studies showed that the autotetraploid had higher percentages (>30%) of abnormal chromosome behavior and aberrant meiocytes (>50%) during meiosis than did the diploid rice hybrid control. Analysis of gene expression profiles revealed 1,888 genes that were differentially expressed between the autotetraploid and diploid hybrid lines at the meiotic stage, among which 889 and 999 were up- and down-regulated, respectively. Of the 999 down-regulated genes, 940 were associated with the combined effect of polyploidy and pollen sterility loci interactions (IPE). Gene Ontology enrichment analysis identified a prominent functional gene class consisting of seven genes related to photosystem I (Gene Ontology 0009522). Moreover, 55 meiosis-related or meiosis stage-specific genes were associated with IPE in autotetraploid rice, including Os02g0497500, which encodes a DNA repair-recombination protein, and Os02g0490000, which encodes a component of the ubiquitin-proteasome pathway. These results suggest that polyploidy enhances epistatic interactions between alleles of pollen sterility loci, thereby altering the expression profiles of important meiosis-related or meiosis stage-specific genes and resulting in high pollen sterility.

The Rice Floral Repressor Early flowering1 Affects Spikelet Fertility By Modulating Gibberellin Signaling

BACKGROUND

Gibberellic acid (GA; or gibberellin) affects the development of floral organs, especially anthers and pollen, and perturbation of development of male floral organs can cause sterility. Many studies of GA signaling have concentrated on anther development, but the effect of GA on grain production remains to be examined.

RESULTS

Using a cross of 'Milyang23 (M23)', which has a functional allele of Early flowering1 (EL1), and 'H143', which has a nonfunctional el1 allele, we generated heterogeneous inbred family-near isogenic lines (HNILs) that are homozygous for EL1 [HNIL(M23)] or el1 [HNIL(H143)]. Here, we found that HNIL(H143) exhibited anther deformities and low pollen viability. The expression of GAMYB, a major activator of GA signaling, and its downstream genes CYP703A3 and KAR, mainly involved in pollen formation, increased abnormally during spikelet development; this activation of GA signaling may cause the sterility. To confirm the negative effect of the el1 mutation on spikelet fertility, we examined a line carrying a T-DNA insertion el1 mutant [hereafter ZH11(el1)] and its parental cultivar 'Zhonghua11 (ZH11)'. ZH11(el1) showed nearly identical defects in anther development and pollen viability as HNIL(H143), leading to decreased seed setting rate. However, the elite japonica cultivar Koshihikari, which has a nonfunctional el1 allele for early flowering in long days, produces fertile spikelets and normal grain yields, like other elite japonica cultivars. This indicates that as-yet-unknown regulator(s) that can overcome the male sterile phenotype of the el1 mutation must have been introduced into Koshihikari.

CONCLUSIONS

The el1 mutation contributes to early flowering in japonica rice under long days but fails to limit GA signaling, thus negatively affecting spikelet fertility, which results in a loss of grain yield. Thus, EL1 is essential for photoperiod sensitivity in flowering as well as spikelet fertility in grain production.

Flowers and mycorrhizal roots--closer than we think?

Roots and flowers are formed at the extreme ends of plants and they differ in almost every aspect of their development and function; even so, they exhibit surprising molecular commonalities. For example, the calcium and calmodulin-dependent protein kinase (CCaMK) plays a central role in root symbioses with fungi and bacteria, but is also highly expressed in developing anthers. Moreover, independent evidence from transcriptomics, phylogenomics, and genetics reveals common developmental elements in root symbioses and reproductive development. We discuss the significance of these overlaps, and we argue that an integrated comparative view of the two phenomena will stimulate research and provide new insight, not only into shared components, but also into the specific aspects of anther development and root symbioses.

A decade of pollen transcriptomics

Overview of pollen transcriptome studies. Pollen development is driven by gene expression, and knowledge of the molecular events underlying this process has undergone a quantum leap in the last decade through studies of the transcriptome. Here, we outline historical evidence for male haploid gene expression and review the wealth of pollen transcriptome data now available. Knowledge of the transcriptional capacity of pollen has progressed from genetic studies to the direct analysis of RNA and from gene-by-gene studies to analyses on a genomic scale. Microarray and/or RNA-seq data can now be accessed for all phases and cell types of developing pollen encompassing 10 different angiosperms. These growing resources have accelerated research and will undoubtedly inspire new directions and the application of system-based research into the mechanisms that govern the development, function and evolution of angiosperm pollen.

Comprehensive network analysis of genes expressed in human oropharyngeal cancer

PURPOSE

Oropharyngeal cancer (OPC) is the eighth most common cancer worldwide, however the genes involved in the development of OPC have been reported few. We constructed a co-expression network to extend knowledge of the molecular biomarkers in OPC development.

MATERIALS AND METHODS

Microarray data of HPV-active, -inactive, -negative OPC and normal benign tissue (uvula, tonsil) (Series GSE55550) were retrieved from NCBI GEO DataSets. We performed co-expression analysis of OPC transcriptome data by the Pearson correlation coefficient (PCC) method with the mutual rank (MR)-based cut-off using 13 guide genes.

RESULTS

The OPC subnetwork contained three clusters: cell cycle (62 node genes and 125 edge genes), immune system (44 node genes and 70 edge genes) and organ morphogenesis (128 node gene and 215 edge genes) process separately.

CONCLUSION

Our co-expression analysis includes separated transcriptomes of OPC, which is a useful resource for OPC researchers to elucidate important and complex biological events, to prevent and to predict cancer.

Comparative cytological and transcriptomic analysis of pollen development in autotetraploid and diploid rice

Autotetraploid rice has greater genetic variation and higher vigor than diploid rice, but low pollen fertility is one of the major reasons for low yield of autotetraploid rice. Very little is known about the molecular mechanisms of low pollen fertility of autotetraploid rice. In this study, cytological observations and microarray analysis were used to assess the genetic variation during pollen development in autotetraploid and diploid rice. Many abnormal chromosome behaviors, such as mutivalents, lagged chromosomes, asynchronous cell division, and so on, were found during meiosis in autotetraploid. Microsporogenesis and microgametogenesis in autotetraploid rice was similar to diploid rice, but many different kinds of abnormalities, including microspores degeneration, multi-aperture, and abnormal cell walls, were found in autotetraploid rice. Compared with diploid rice, a total of 1,251 genes were differentially expressed in autotetraploid rice in pollen transcriptome, among them 1,011 and 240 genes were up-regulated and down-regulated, respectively. 124 and 6 genes were co-up-regulated and co-down-regulated during three pollen development stages, respectively. These results suggest that polyploidy induced up-regulation for most of the genes during pollen development. Quantitative RT-PCR was done to validate 12 differentially expressed genes selected from functional categories based on the gene ontology analysis. These stably expressed genes not only related to the pollen development genes, but also involved in cell metabolism, cell physiology, binding, catalytic activity, molecular transducer activity, and transcription regulator activity. The present study suggests that differential expression of some key genes may lead to complex gene regulation and abnormal pollen development in autotetraploid rice.

Hypericum perforatum hydroxyalkylpyrone synthase involved in sporopollenin biosynthesis--phylogeny, site-directed mutagenesis, and expression in nonanther tissues

Anther-specific chalcone synthase-like enzyme (ASCL), an ancient plant type III polyketide synthase, is involved in the biosynthesis of sporopollenin, the stable biopolymer found in the exine layer of the wall of a spore or pollen grain. The gene encoding polyketide synthase 1 from Hypericum perforatum (HpPKS1) was previously shown to be expressed mainly in young flower buds, but also in leaves and other tissues at lower levels. Angiosperm ASCLs, identified by sequence and phylogenetic analyses, are divided into two sister clades, the Ala-clade and the Val-clade, and HpPKS1 belongs to the Ala-clade. Recombinant HpPKS1 produced triketide and, to a lesser extent, tetraketide alkylpyrones from medium-chain (C6) to very long-chain (C24) fatty acyl-CoA substrates. Like other ASCLs, HpPKS1 also preferred hydroxyl fatty acyl-CoA esters over the analogous unsubstituted fatty acyl-CoA esters. To study the structural basis of the substrate preference, mutants of Ala200 and Ala215 at the putative active site and Arg202 and Asp211 at the modeled acyl-binding tunnel were constructed. The A200T/A215Q mutant accepted decanoyl-CoA, a poor substrate for the wild-type enzyme, possibly because of active site constriction by bulkier substitutions. The substrate preference of the A215V and A200T/A215Q mutants shifted toward nonhydroxylated, medium-chain to long-chain fatty acyl-CoA substrates. The R202L/D211V double mutant was selective for acyl-CoA with chain lengths of C16-C18, and showed a diminished preference for the hydroxylated acyl-CoA substrates. Transient upregulation by abscisic acid and downregulation by jasmonic acid and wounding suggested that HpPKS1, and possibly other Ala-clade ASCLs, may be involved in the biosynthesis of minor cell wall components in nonanther tissues.

Transcriptomes and proteomes define gene expression progression in pre-meiotic maize anthers

Plants lack a germ line; consequently, during reproduction adult somatic cells within flowers must switch from mitotic proliferation to meiosis. In maize (Zea mays L.) anthers, hypoxic conditions in the developing tassel trigger pre-meiotic competence in the column of pluripotent progenitor cells in the center of anther lobes, and within 24 hr these newly specified germinal cells have patterned their surrounding neighbors to differentiate as the first somatic niche cells. Transcriptomes were analyzed by microarray hybridization in carefully staged whole anthers during initial specification events, after the separation of germinal and somatic lineages, during the subsequent rapid mitotic proliferation phase, and during final pre-meiotic germinal and somatic cell differentiation. Maize anthers exhibit a highly complex transcriptome constituting nearly three-quarters of annotated maize genes, and expression patterns are dynamic. Laser microdissection was applied to begin assigning transcripts to tissue and cell types and for comparison to transcriptomes of mutants defective in cell fate specification. Whole anther proteomes were analyzed at three developmental stages by mass spectrometric peptide sequencing using size-fractionated proteins to evaluate the timing of protein accumulation relative to transcript abundance. New insights include early and sustained expression of meiosis-associated genes (77.5% of well-annotated meiosis genes are constitutively active in 0.15 mm anthers), an extremely large change in transcript abundances and types a few days before meiosis (including a class of 1340 transcripts absent specifically at 0.4 mm), and the relative disparity between transcript abundance and protein abundance at any one developmental stage (based on 1303 protein-to-transcript comparisons).

The meiotic transcriptome architecture of plants

Although a number of genes that play key roles during the meiotic process have been characterized in great detail, the whole process of meiosis is still not completely unraveled. To gain insight into the bigger picture, large-scale approaches like RNA-seq and microarray can help to elucidate the transcriptome landscape during plant meiosis, discover co-regulated genes, enriched processes, and highly expressed known and unknown genes which might be important for meiosis. These high-throughput studies are gaining more and more popularity, but their beginnings in plant systems reach back as far as the 1960's. Frequently, whole anthers or post-meiotic pollen were investigated, while less data is available on isolated cells during meiosis, and only few studies addressed the transcriptome of female meiosis. For this review, we compiled meiotic transcriptome studies covering different plant species, and summarized and compared their key findings. Besides pointing to consistent as well as unique discoveries, we finally draw conclusions what can be learned from these studies so far and what should be addressed next.

Polyamine oxidase 7 is a terminal catabolism-type enzyme in Oryza sativa and is specifically expressed in anthers

Polyamine oxidase (PAO), which requires FAD as a cofactor, functions in polyamine catabolism. Plant PAOs are classified into two groups based on their reaction modes. The terminal catabolism (TC) reaction always produces 1,3-diaminopropane (DAP), H2O2, and the respective aldehydes, while the back-conversion (BC) reaction produces spermidine (Spd) from tetraamines, spermine (Spm) and thermospermine (T-Spm) and/or putrescine from Spd, along with 3-aminopropanal and H2O2. The Oryza sativa genome contains seven PAO-encoded genes termed OsPAO1-OsPAO7. To date, we have characterized four OsPAO genes. The products of these genes, i.e. OsPAO1, OsPAO3, OsPAO4 and OsPAO5, catalyze BC-type reactions. Whereas OsPAO1 remains in the cytoplasm, the other three PAOs localize to peroxisomes. Here, we examined OsPAO7 and its gene product. OsPAO7 shows high identity to maize ZmPAO1, the best characterized plant PAO having TC-type activity. OsPAO7 seems to remain in a peripheral layer of the plant cell with the aid of its predicted signal peptide and transmembrane domain. Recombinant OsPAO7 prefers Spm and Spd as substrates, and it produces DAP from both substrates in a time-dependent manner, indicating that OsPAO7 is the first TC-type enzyme identified in O. sativa. The results clearly show that two types of PAOs co-exist in O. sativa. Furthermore, OsPAO7 is specifically expressed in anthers, with an expressional peak at the bicellular pollen stage. The physiological function of OsPAO7 in anthers is discussed.

The bHLH142 Transcription Factor Coordinates with TDR1 to Modulate the Expression of EAT1 and Regulate Pollen Development in Rice

Male sterility plays an important role in F1 hybrid seed production. We identified a male-sterile rice (Oryza sativa) mutant with impaired pollen development and a single T-DNA insertion in the transcription factor gene bHLH142. Knockout mutants of bHLH142 exhibited retarded meiosis and defects in tapetal programmed cell death. RT-PCR and in situ hybridization analyses showed that bHLH142 is specifically expressed in the anther, in the tapetum, and in meiocytes during early meiosis. Three basic helix-loop-helix transcription factors, UDT1 (bHLH164), TDR1 (bHLH5), and EAT1/DTD1 (bHLH141) are known to function in rice pollen development. bHLH142 acts downstream of UDT1 and GAMYB but upstream of TDR1 and EAT1 in pollen development. In vivo and in vitro assays demonstrated that bHLH142 and TDR1 proteins interact. Transient promoter assays demonstrated that regulation of the EAT1 promoter requires bHLH142 and TDR1. Consistent with these results, 3D protein structure modeling predicted that bHLH142 and TDR1 form a heterodimer to bind to the EAT1 promoter. EAT1 positively regulates the expression of AP37 and AP25, which induce tapetal programmed cell death. Thus, in this study, we identified bHLH142 as having a pivotal role in tapetal programmed cell death and pollen development.

Evaluation of rice promoters conferring pollen-specific expression in a heterologous system, Arabidopsis

Promoters can direct gene expression specifically to targeted tissues or cells. Effective with both crop species and model plant systems, these tools can help researchers overcome the practical obstacles associated with transgenic protocols. Here, we identified promoters that allow one to target the manipulation of gene expression during pollen development. Utilizing published transcriptomic databases for rice, we investigated the promoter activity of selected genes in Arabidopsis. From various microarray datasets, including those for anthers and pollen grains at different developmental stages, we selected nine candidate genes that showed high levels of expression in the late stages of rice pollen development. We named these Oryza sativa late pollen-specific genes. Their promoter regions contained various cis-acting elements that could be responsible for anther-/pollen-specific expression. Promoter::GUS-GFP reporters were constructed and introduced into Arabidopsis plants. Histochemical GUS staining revealed that six of the nine rice promoters conferred strong GUS expression that was restricted to the anthers in Arabidopsis. Further analysis showed that although the GUS signals were not detected at the unicellular stage, they strengthened in the bicellular or tricellular stages, peaking at the mature pollen stage. This paralleled their transcriptomic profiles in rice. Based on our results, we proposed that these six rice promoters, which are active in the late stages of pollen formation in the dicot Arabidopsis, can aid molecular breeders in generating new varieties of a monocot plant, rice.

Maize germinal cell initials accommodate hypoxia and precociously express meiotic genes

In flowering plants, anthers are the site of de novo germinal cell specification, male meiosis, and pollen development. Atypically, anthers lack a meristem. Instead, both germinal and somatic cell types differentiate from floral stem cells packed into anther lobes. To better understand anther cell fate specification and to provide a resource for the reproductive biology community, we isolated cohorts of germinal and somatic initials from maize anthers within 36 h of fate acquisition, identifying 815 specific and 1714 significantly enriched germinal transcripts, plus 2439 specific and 2112 significantly enriched somatic transcripts. To clarify transcripts involved in cell differentiation, we contrasted these profiles to anther primordia prior to fate specification and to msca1 anthers arrested in the first step of fate specification and hence lacking normal cell types. The refined cell-specific profiles demonstrated that both germinal and somatic cell populations differentiate quickly and express unique transcription factor sets; a subset of transcript localizations was validated by in situ hybridization. Surprisingly, germinal initials starting 5 days of mitotic divisions were enriched significantly in >100 transcripts classified in meiotic processes that included recombination and synapsis, along with gene sets involved in RNA metabolism, redox homeostasis, and cytoplasmic ATP generation. Enrichment of meiotic-specific genes in germinal initials challenges current dogma that the mitotic to meiotic transition occurs later in development during pre-meiotic S phase. Expression of cytoplasmic energy generation genes suggests that male germinal cells accommodate hypoxia by diverting carbon away from mitochondrial respiration into alternative pathways that avoid producing reactive oxygen species (ROS).

Genetic mechanisms of postzygotic reproductive isolation: An epistatic network in rice

Products of interspecific crosses often show abnormal phenotypes such as sterility, weakness and inviability. These phenomena play an important role in speciation as mechanisms of postzygotic reproductive isolation (RI). During the past two decades, genetics studies in rice have characterized a number of gene loci responsible for postzygotic RI. I have identified 10 loci including three sets of epistatic networks in a single inter-subspecific cross (Oryza sativa ssp. indica × japonica). These results suggest that RI genes cause developmental dysfunction of vegetative and/or reproductive organs through a variety of molecular pathways. The latest molecular studies demonstrated that hybrid incompatibility is mainly due to deleterious interactions caused by species-specific mutations of two or more genes, mediated by proteins acting within the same molecular pathway. Because genetic interactions provide a perspective on gene function, epistatic networks are a key to the understanding of the molecular basis of postzygotic RI. In this review, I focus on recent progress in postzygotic RI studies in rice and discuss the evolutionary significance as well as implications for improving rice productivity.

The ATP-binding cassette transporter OsABCG15 is required for anther development and pollen fertility in rice

Plant male reproductive development is a complex biological process, but the underlying mechanism is not well understood. Here, we characterized a rice (Oryza sativa L.) male sterile mutant. Based on map-based cloning and sequence analysis, we identified a 1,459-bp deletion in an adenosine triphosphate (ATP)-binding cassette (ABC) transporter gene, OsABCG15, causing abnormal anthers and male sterility. Therefore, we named this mutant osabcg15. Expression analysis showed that OsABCG15 is expressed specifically in developmental anthers from stage 8 (meiosis II stage) to stage 10 (late microspore stage). Two genes CYP704B2 and WDA1, involved in the biosynthesis of very-long-chain fatty acids for the establishment of the anther cuticle and pollen exine, were downregulated in osabcg15 mutant, suggesting that OsABCG15 may play a key function in the processes related to sporopollenin biosynthesis or sporopollenin transfer from tapetal cells to anther locules. Consistently, histological analysis showed that osabcg15 mutants developed obvious abnormality in postmeiotic tapetum degeneration, leading to rapid degredation of young microspores. The results suggest that OsABCG15 plays a critical role in exine formation and pollen development, similar to the homologous gene of AtABCG26 in Arabidopsis. This work is helpful to understand the regulatory network in rice anther development.

Developmental transcriptional profiling reveals key insights into Triticeae reproductive development

Despite their importance, there remains a paucity of large-scale gene expression-based studies of reproductive development in species belonging to the Triticeae. As a first step to address this deficiency, a gene expression atlas of triticale reproductive development was generated using the 55K Affymetrix GeneChip(®) wheat genome array. The global transcriptional profiles of the anther/pollen, ovary and stigma were analyzed at concurrent developmental stages, and co-expressed as well as preferentially expressed genes were identified. Data analysis revealed both novel and conserved regulatory factors underlying Triticeae floral development and function. This comprehensive resource rests upon detailed gene annotations, and the expression profiles are readily accessible via a web browser.

Transcriptome comparison and gene coexpression network analysis provide a systems view of citrus response to 'Candidatus Liberibacter asiaticus' infection

BACKGROUND

Huanglongbing (HLB) is arguably the most destructive disease for the citrus industry. HLB is caused by infection of the bacterium, Candidatus Liberibacter spp. Several citrus GeneChip studies have revealed thousands of genes that are up- or down-regulated by infection with Ca. Liberibacter asiaticus. However, whether and how these host genes act to protect against HLB remains poorly understood.

RESULTS

As a first step towards a mechanistic view of citrus in response to the HLB bacterial infection, we performed a comparative transcriptome analysis and found that a total of 21 Probesets are commonly up-regulated by the HLB bacterial infection. In addition, a number of genes are likely regulated specifically at early, late or very late stages of the infection. Furthermore, using Pearson correlation coefficient-based gene coexpression analysis, we constructed a citrus HLB response network consisting of 3,507 Probesets and 56,287 interactions. Genes involved in carbohydrate and nitrogen metabolic processes, transport, defense, signaling and hormone response were overrepresented in the HLB response network and the subnetworks for these processes were constructed. Analysis of the defense and hormone response subnetworks indicates that hormone response is interconnected with defense response. In addition, mapping the commonly up-regulated HLB responsive genes into the HLB response network resulted in a core subnetwork where transport plays a key role in the citrus response to the HLB bacterial infection. Moreover, analysis of a phloem protein subnetwork indicates a role for this protein and zinc transporters or zinc-binding proteins in the citrus HLB defense response.

CONCLUSION

Through integrating transcriptome comparison and gene coexpression network analysis, we have provided for the first time a systems view of citrus in response to the Ca. Liberibacter spp. infection causing HLB.

Molecular control of male reproductive development and pollen fertility in rice

Anther development and male fertility are essential biological processes for flowering plants and are important for crop seed production. Genetic manipulation of male fertility/sterility is critical for crop hybrid breeding. Rice (Oryza sativa L.) male sterility phenotypes, including genic male sterility, hybrid male sterility, and cytoplasmic male sterility, are generally caused by mutations of fertility-related genes, by incompatible interactions between divergent allelic or non-allelic genes, or by genetic incompatibilities between cytoplasmic and nuclear genomes. Here, we review the recent advances in the molecular basis of anther development and male fertility-sterility conversion in specific genetic backgrounds, and the interactions with certain environmental factors. The highlighted findings in this review have significant implications in both basic studies and rice genetic improvement.

Molecular genetics, physiology and biology of self-incompatibility in Brassicaceae

Self-incompatibility (SI) is defined as the inability to produce zygotes after self-pollination in a fertile hermaphrodite plant, which has stamens and pistils in the same flower. This structural organization of the hermaphrodite flower increases the risk of self-pollination, leading to low genetic diversity. To avoid this problem plants have established several pollination systems, among which the most elegant system is surely SI. The SI trait can be observed in Brassica crops, including cabbage, broccoli, turnip and radish. To produce hybrid seed of these crops efficiently, the SI trait has been employed in an agricultural context. From another point of view, the recognition reaction of SI during pollen-stigma interaction is an excellent model system for cell-cell communication and signal transduction in higher plants. In this review, we describe the molecular mechanisms of SI in Brassicaceae, which have been dissected by genetic, physiological, and biological approaches, and we discuss the future prospects in relation to associated scientific fields and new technologies.

Inferring the Brassica rapa Interactome Using Protein-Protein Interaction Data from Arabidopsis thaliana

Following successful completion of the Brassica rapa sequencing project, the next step is to investigate functions of individual genes/proteins. For Arabidopsis thaliana, large amounts of protein-protein interaction (PPI) data are available from the major PPI databases (DBs). It is known that Brassica crop species are closely related to A. thaliana. This provides an opportunity to infer the B. rapa interactome using PPI data available from A. thaliana. In this paper, we present an inferred B. rapa interactome that is based on the A. thaliana PPI data from two resources: (i) A. thaliana PPI data from three major DBs, BioGRID, IntAct, and TAIR. (ii) ortholog-based A. thaliana PPI predictions. Linking between B. rapa and A. thaliana was accomplished in three complementary ways: (i) ortholog predictions, (ii) identification of gene duplication based on synteny and collinearity, and (iii) BLAST sequence similarity search. A complementary approach was also applied, which used known/predicted domain-domain interaction data. Specifically, since the two species are closely related, we used PPI data from A. thaliana to predict interacting domains that might be conserved between the two species. The predicted interactome was investigated for the component that contains known A. thaliana meiotic proteins to demonstrate its usability.