Dissecting Tumor Antigens and Immune Subtypes of Glioma to Develop mRNA Vaccine.

Background: Nowadays, researchers are leveraging the mRNA-based vaccine technology used to develop personalized immunotherapy for cancer. However, its application against glioma is still in its infancy. In this study, the applicable candidates were excavated for mRNA vaccine treatment in the perspective of immune regulation, and suitable glioma recipients with corresponding immune subtypes were further investigated. Methods: The RNA-seq data and clinical information of 702 and 325 patients were recruited from TCGA and CGGA, separately. The genetic alteration profile was visualized and compared by cBioPortal. Then, we explored prognostic outcomes and immune correlations of the selected antigens to validate their clinical relevance. The prognostic index was measured via GEPIA2, and infiltration of antigen-presenting cells (APCs) was calculated and visualized by TIMER. Based on immune-related gene expression, immune subtypes of glioma were identified using consensus clustering analysis. Moreover, the immune landscape was visualized by graph learning-based dimensionality reduction analysis. Results: Four glioma antigens, namely ANXA5, FKBP10, MSN, and PYGL, associated with superior prognoses and infiltration of APCs were selected. Three immune subtypes IS1-IS3 were identified, which fundamentally differed in molecular, cellular, and clinical signatures. Patients in subtypes IS2 and IS3 carried immunologically cold phenotypes, whereas those in IS1 carried immunologically hot phenotype. Particularly, patients in subtypes IS3 and IS2 demonstrated better outcomes than that in IS1. Expression profiles of immune checkpoints and immunogenic cell death (ICD) modulators showed a difference among IS1-IS3 tumors. Ultimately, the immune landscape of glioma elucidated considerable heterogeneity not only between individual patients but also within the same immune subtype. Conclusions: ANXA5, FKBP10, MSN, and PYGL are identified as potential antigens for anti-glioma mRNA vaccine production, specifically for patients in immune subtypes 2 and 3. In summary, this study may shed new light on the promising approaches of immunotherapy, such as devising mRNA vaccination tailored to applicable glioma recipients.

Genetic Diversity Relationship Between Grain Quality and Appearance in Rice.

Grain quality is an important breeding objective in rice, and the appearance of the grain also affects its commercial value in the market. The aim of this study was to decode the rice grain qualities and appearances, such as gelatinization temperature (GT), amylose content (AC), grain protein content (GPC), pericarp color (PC), length/width ratio (LWR), and grain volume (GV) at phenotypic and genetic levels, as well as the relationships among them. A genome-wide association study (GWAS) was used to identify the quantitative trait locus (QTLs) associated with the target traits using mixed linear model (MLM) and Bayesian-information and linkage-disequilibrium iteratively nested keyway (BLINK) methods. In general, AC was negatively correlated with GPC and GV, while it was positively correlated with LWR and PC. GPC was positively correlated with LWR. Using the rice diversity panel 1 (RDP1) population, we identified 11, 6, 2, 7, 11, and 6 QTLs associated with GT, AC, GPC, PC, LWR, and GV, respectively. Five germplasm lines, superior in grain qualities and appearances for basic breeding materials or improvement, were identified. Notably, an F-box gene OsFbox394 was located in the linkage disequilibrium (LD) block of qLWR7-2, which specifically expresses in endosperm and seed tissues, suggesting that this gene may regulate the seed development in rice grain. Besides, different haplotypes of OsHyPRP45 showed significant differences in AC, indicating that this gene may be related to AC in rice grain.

Multi-locus genome-wide association studies for five yield-related traits in rice.

BACKGROUND: Improving the overall production of rice with high quality is a major target of breeders. Mining potential yield-related loci have been geared towards developing efficient rice breeding strategies. In this study, one single-locus genome-wide association studies (SL-GWAS) method (MLM) in conjunction with five multi-locus genome-wide association studies (ML-GWAS) approaches (mrMLM, FASTmrMLM, pLARmEB, pKWmEB, and ISIS EM-BLASSO) were conducted in a panel consisting of 529 rice core varieties with 607,201 SNPs. RESULTS: A total of 152, 106, 12, 111, and 64 SNPs were detected by the MLM model associated with the five yield-related traits, namely grain length (GL), grain width (GW), grain thickness (GT), thousand-grain weight (TGW), and yield per plant (YPP), respectively. Furthermore, 74 significant quantitative trait nucleotides (QTNs) were presented across at least two ML-GWAS methods to be associated with the above five traits successively. Finally, 20 common QTNs were simultaneously discovered by both SL-GWAS and ML-GWAS methods. Based on genome annotation, gene expression analysis, and previous studies, two candidate key genes (LOC_Os09g02830 and LOC_Os07g31450) were characterized to affect GW and TGW, separately. CONCLUSIONS: These outcomes will provide an indication for breeding high-yielding rice varieties in the immediate future.

Global Analysis of UDP Glucose Pyrophosphorylase (UDPGP) Gene Family in Plants: Conserved Evolution Involved in Cell Death.

UDP glucose pyrophosphorylase (UDPGP) family genes have been reported to play essential roles in cell death or individual survival. However, a systematic analysis on UDPGP gene family has not been performed yet. In this study, a total of 454 UDPGP proteins from 76 different species were analyzed. The analyses of the phylogenetic tree and orthogroups divided UDPGPs into three clades, including UDP-N-acetylglucosamine pyrophosphorylase (UAP), UDP-glucose pyrophosphorylase (UGP, containing UGP-A and UGP-B), and UDP-sugar pyrophosphorylase (USP). The evolutionary history of the UDPGPs indicated that the members of UAP, USP, and UGP-B were relatively conserved while varied in UGP-A. Homologous sequences of UGP-B and USP were found only in plants. The expression profile of UDPGP genes in Oryza sativa was mainly motivated under jasmonic acid (JA), abscisic acid (ABA), cadmium, and cold treatments, indicating that UDPGPs may play an important role in plant development and environment endurance. The key amino acids regulating the activity of UDPGPs were analyzed, and almost all of them were located in the NB-loop, SB-loop, or conserved motifs. Analysis of the natural variants of UDPGPs in rice revealed that only a few missense mutants existed in coding sequences (CDSs), and most of the resulting variations were located in the non-motif sites, indicating the conserved structure and function of UDPGPs in the evolution. Furthermore, alternative splicing may play a key role in regulating the activity of UDPGPs. The spatial structure prediction, enzymatic analysis, and transgenic verification of UAP isoforms illustrated that the loss of N- and C-terminal sequences did not affect the overall 3D structures, but the N- and C-terminal sequences are important for UAP genes to maintain their enzymatic activity. These results revealed a conserved UDPGP gene family and provided valuable information for further deep functional investigation of the UDPGP gene family in plants.

Uncovering the genetic mechanisms regulating panicle architecture in rice with GPWAS and GWAS.

BACKGROUND: The number of panicles per plant, number of grains per panicle, and 1000-grain weight are important factors contributing to the grain yield per plant in rice. The Rice Diversity Panel 1 (RDP1) contains a total of 421 purified, homozygous rice accessions representing diverse genetic variations within O. sativa. The release of High-Density Rice Array (HDRA, 700 k SNPs) dataset provides a new opportunity to discover the genetic variants of panicle architectures in rice. RESULTS: In this report, a new method genome-phenome wide association study (GPWAS) was performed with 391 individuals and 27 traits derived from RDP1 to scan the relationship between the genes and multi-traits. A total of 1985 gene models were linked to phenomic variation with a p-value cutoff of 4.49E-18. Besides, 406 accessions derived from RDP1 with 411,066 SNPs were used to identify QTLs associated with the total spikelets number per panicle (TSNP), grain number per panicle (GNP), empty grain number per panicle (EGNP), primary branch number (PBN), panicle length (PL), and panicle number per plant (PN) by GLM, MLM, FarmCPU, and BLINK models for genome-wide association study (GWAS) analyses. A total of 18, 21, 18, 17, 15, and 17 QTLs were identified tightly linked with TSNP, GNP, EGNP, PBN, PL, and PN, respectively. Then, a total of 23 candidate genes were mapped simultaneously using both GWAS and GPWAS methods, composed of 6, 4, 5, 4, and 4 for TSNP, GNP, EGNP, PBN, and PL. Notably, one overlapped gene (Os01g0140100) were further investigated based on the haplotype and gene expression profile, indicating this gene might regulate the TSNP or panicle architecture in rice. CONCLUSIONS: Nearly 30 % (30/106) QTLs co-located with the previous published genes or QTLs, indicating the power of GWAS. Besides, GPWAS is a new method to discover the relationship between genes and traits, especially the pleiotropy genes. Through comparing the results from GWAS and GPWAS, we identified 23 candidate genes related to panicle architectures in rice. This comprehensive study provides new insights into the genetic basis controlling panicle architectures in rice, which lays a foundation in rice improvement.

Feeding Arsenic-Containing Rice Bran to Growing Pigs: Growth Performance, Arsenic Tissue Distribution, and Arsenic Excretion.

This research was conducted to study the growth performance, arsenic (As) tissue distribution, and As excretion of pigs fed As-containing rice bran. Twenty gilts (26.3 kg) were randomly assigned to 3 dietary treatments (n = 6 or 7) with Diets I, II, and III containing 0, 36.7, and 73.5% rice bran and 0, 306, and 612 ppb As, respectively. Pigs were fed for 6 weeks, and their growth performance and daily activities were examined. Fecal, blood, and hair samples were collected immediately before and after the 6-weeks. At the end of the 6-weeks, pigs were slaughtered; the liver, kidney, muscle, and urine samples were collected. No pig showed any unhealthy signs throughout the trial. The average daily feed intake, average daily gain, and final body weight of Diet III pigs were lower (p ≤ 0.001) than Diet I pigs. The gain to feed ratios were not different among the treatments. The fecal, hair, kidney, and urinary As concentrations of both Diets II and III pigs were higher than Diet I pigs. The hair As concentration of Diet III pigs was higher than Diet II pigs, but no difference was found in the fecal, urinary, kidney, or muscle As concentrations between Diets II and III pigs. The blood and muscle As concentrations were below 10 ppb. These results suggest that 73.5% dietary rice bran inclusion compromised growth performance, whereas the 36.7% inclusion did not. The fecal As data imply that dietary As was poorly absorbed by the gastrointestinal tract. The tissue As data indicate that the absorbed As was rapidly cleared from the blood with some retained in various organs and others eliminated via urine. The hair As concentration was much higher than that of liver and kidney. The muscle As data suggest that the pork produced from the pigs fed a typical As-containing rice bran as used in this study is safe for human consumption.

Genome-wide association studies of ionomic and agronomic traits in USDA mini core collection of rice and comparative analyses of different mapping methods.

BACKGROUND: Rice is an important human staple food vulnerable to heavy metal contamination leading to serious concerns. High yield with low heavy metal contamination is a common but highly challenging goal for rice breeders worldwide due to lack of genetic knowledge and markers. RESULTS: To identify candidate QTLs and develop molecular markers for rice yield and heavy metal content, a total of 191 accessions from the USDA Rice mini-core collection with over 3.2 million SNPs were employed to investigate the QTLs. Sixteen ionomic and thirteen agronomic traits were analyzed utilizing two univariate (GLM and MLM) and two multivariate (MLMM and FarmCPU) GWAS methods. 106, 47, and 97 QTLs were identified for ionomics flooded, ionomics unflooded, and agronomic traits, respectively, with the criterium of p-value < 1.53 × 10- 8, which was determined by the Bonferroni correction for p-value of 0.05. While 49 (~ 20%) of the 250 QTLs were coinciding with previously reported QTLs/genes, about 201 (~ 80%) were new. In addition, several new candidate genes involved in ionomic and agronomic traits control were identified by analyzing the DNA sequence, gene expression, and the homologs of the QTL regions. Our results further showed that each of the four GWAS methods can identify unique as well as common QTLs, suggesting that using multiple GWAS methods can complement each other in QTL identification, especially by combining univariate and multivariate methods. CONCLUSIONS: While 49 previously reported QTLs/genes were rediscovered, over 200 new QTLs for ionomic and agronomic traits were found in the rice genome. Moreover, multiple new candidate genes for agronomic and ionomic traits were identified. This research provides novel insights into the genetic basis of both ionomic and agronomic variations in rice, establishing the foundation for marker development in breeding and further investigation on reducing heavy-metal contamination and improving crop yields. Finally, the comparative analysis of the GWAS methods showed that each method has unique features and different methods can complement each other.

Comprehensive Analysis of the Lysine Succinylome and Protein Co-modifications in Developing Rice Seeds.

Lysine succinylation has been recognized as a post-translational modification (PTM) in recent years. It is plausible that succinylation may have a vaster functional impact than acetylation because of bulkier structural changes and more significant charge differences on the modified lysine residue. Currently, however, the quantity and identity of succinylated proteins and their corresponding functions in cereal plants remain largely unknown. In this study, we estimated the native succinylation occupancy on lysine was between 2% to 10% in developing rice seeds. Eight hundred fifty-four lysine succinylation sites on 347 proteins have been identified by a thorough investigation in developing rice seeds. Six motifs were revealed as preferred amino acid sequence arrangements for succinylation sites, and a noteworthy motif preference was identified in proteins associated with different biological processes, molecular functions, pathways, and domains. Remarkably, heavy succinylation was detected on major seed storage proteins, in conjunction with critical enzymes involved in central carbon metabolism and starch biosynthetic pathways for rice seed development. Meanwhile, our results showed that the modification pattern of in vitro nonenzymatically succinylated proteins was different from those of the proteins isolated from cells in Western blots, suggesting that succinylation is not generated via nonenzymatic reaction in the cells, at least not completely. Using the acylation data obtained from the same rice tissue, we mapped many sites harboring lysine succinylation, acetylation, malonylation, crotonylation, and 2-hydroxisobutyrylation in rice seed proteins. A striking number of proteins with multiple modifications were shown to be involved in critical metabolic events. Given that these modification moieties are intermediate products of multiple cellular metabolic pathways, these targeted lysine residues may mediate the crosstalk between different metabolic pathways via modifications by different moieties. Our study exhibits a platform for extensive investigation of molecular networks administrating cereal seed development and metabolism via PTMs.

Proteome-wide lysine acetylation identification in developing rice (Oryza sativa) seeds and protein co-modification by acetylation, succinylation, ubiquitination, and phosphorylation.

Protein lysine acetylation is a highly conserved post-translational modification with various biological functions. However, only a limited number of acetylation sites have been reported in plants, especially in cereals, and the function of non-histone protein acetylation is still largely unknown. In this report, we identified 1003 lysine acetylation sites in 692 proteins of developing rice seeds, which greatly extended the number of known acetylated sites in plants. Seven distinguished motifs were detected flanking acetylated lysines. Functional annotation analyses indicated diverse biological processes and pathways engaged in lysine acetylation. Remarkably, we found that several key enzymes in storage starch synthesis pathway and the main storage proteins were heavily acetylated. A comprehensive comparison of the rice acetylome, succinylome, ubiquitome and phosphorylome with available published data was conducted. A large number of proteins carrying multiple kinds of modifications were identified and many of these proteins are known to be key enzymes of vital metabolic pathways. Our study provides extending knowledge of protein acetylation. It will have critical reference value for understanding the mechanisms underlying PTM mediated multiple signal integration in the regulation of metabolism and development in plants.

Malonylome analysis in developing rice (Oryza sativa) seeds suggesting that protein lysine malonylation is well-conserved and overlaps with acetylation and succinylation substantially.

In recent years, lysine malonylation has garnered wide spread interest due to its potential regulatory roles. While studies have been performed in bacteria, mouse, and human, the involvement and the biological function of this modification in plant are still largely unknown. We examined the global proteome profile of lysine malonylation in developing rice seeds using affinity enrichment followed by LC-MS/MS analysis. We identified 421 malonylated lysine sites across 247 proteins. Functional analyses showed predominant presence of malonylated proteins in metabolic processes, including carbon metabolism, glycolysis/gluconeogenesis, TCA cycle, as well as photosynthesis. Malonylation was also detected on enzymes in starch biosynthesis pathway in developing rice seeds. In addition, we found a remarkable overlap among the malonylated, succinylated and acetylated sites identified in rice. Furthermore, malonylation at conserved sites of homologous proteins was observed across organisms of different kingdoms, including mouse, human, and bacteria. Finally, distinct motifs were identified when the rice malonylation sites were analyzed and conserved motifs were observed from bacterium to human and rice. Our results provide an initial understanding of the lysine malonylome in plants. The study has critical reference value for future understanding of the biological function of protein lysine malonylation in plants. BIOLOGICAL SIGNIFICANCE: Lysine malonylation is a newly discovered acylation with functional potential in regulating cellular metabolisms and activities. However, the malonylation status has not been reported in plants. Grain yield and quality, mainly determined during cereal seed development, are closely related to food security, human health and economic value. To evaluate malonylation level in plants and the possible regulatory functions of malonylation in seed development, we conducted comprehensive analyses of malonylome in developing rice seeds. A total of 421 malonylated lysine sites from 247 proteins were identified, which involved in multiple critical metabolic processes, including central carbon metabolism, lipid metabolism, photosynthesis, and starch biosynthesis. We found that charged amino acids, lysine and arginine, were the preferred residues in positions flanking the modified lysines. Highly conserved modification sites on both histone and non-histone proteins were observed among different organisms through sequence alignment analysis. More interestingly, a large number of modification sites shared by malonylation, acetylation and succinylation were identified in rice. The study presents a comprehensive understanding of malonylome in plants, which will serve as an initial platform for further investigation of the functions of lysine malonylation, especially in cereal seeds development.

Proteome-wide Analysis of Lysine 2-hydroxyisobutyrylation in Developing Rice (Oryza sativa) Seeds.

Lysine 2-hydroxyisobutyrylation is a recently identified protein post-translational modification that is known to affect the association between histone and DNA. However, non-histone protein lysine 2-hydroxyisobutyrylation remains largely unexplored. Utilizing antibody-based affinity enrichment and nano-HPLC/MS/MS analyses of 2-hydroxyisobutyrylation peptides, we efficaciously identified 9,916 2-hydroxyisobutyryl lysine sites on 2,512 proteins in developing rice seeds, representing the first lysine 2-hydroxyisobutyrylome dataset in plants. Functional annotation analyses indicated that a wide variety of vital biological processes were preferably targeted by lysine 2-hydroxyisobutyrylation, including glycolysis/gluconeogenesis, TCA cycle, starch biosynthesis, lipid metabolism, protein biosynthesis and processing. Our finding showed that 2-hydroxyisobutyrylated histone sites were conserved across plants, human, and mouse. A number of 2-hydroxyisobutyryl sites were shared with other lysine acylations in both histone and non-histone proteins. Comprehensive analysis of the lysine 2-hydroxyisobutyrylation sites illustrated that the modification sites were highly sequence specific with distinct motifs, and they had less surface accessibility than other lysine residues in the protein. Overall, our study provides the first systematic analysis of lysine 2-hydroxyisobutyrylation proteome in plants, and it serves as an important resource for future investigations of the regulatory mechanisms and functions of lysine 2-hydroxyisobutyrylation.

Proteome Profile of Starch Granules Purified from Rice (Oryza sativa) Endosperm.

Starch is the most important food energy source in cereals. Many of the known enzymes involved in starch biosynthesis are partially or entirely granule-associated in the endosperm. Studying the proteome of rice starch granules is critical for us to further understand the mechanisms underlying starch biosynthesis and packaging of starch granules in rice amyloplasts, consequently for the improvement of rice grain quality. In this article, we developed a protocol to purify starch granules from mature rice endosperm and verified the quality of purified starch granules by microscopy observations, I2 staining, and Western blot analyses. In addition, we found the phenol extraction method was superior to Tris-HCl buffer extraction method with respect to the efficiency in recovery of starch granule associated proteins. LC-MS/MS analysis showed identification of already known starch granule associated proteins with high confidence. Several proteins reported to be involved in starch synthesis in prior genetic studies in plants were also shown to be enriched with starch granules, either directly or indirectly, in our studies. In addition, our results suggested that a few additional candidate proteins may also be involved in starch synthesis. Furthermore, our results indicated that some starch synthesis pathway proteins are subject to protein acetylation modification. GO analysis and KEGG pathway enrichment analysis showed that the identified proteins were mainly located in plastids and involved in carbohydrate metabolism. This study substantially advances the understanding of the starch granule associated proteome in rice and post translational regulation of some starch granule associated proteins.

Genome-wide identification and functional prediction of nitrogen-responsive intergenic and intronic long non-coding RNAs in maize (Zea mays L.).

BACKGROUND: Nitrogen (N) is an essential and often limiting nutrient to plant growth and development. Previous studies have shown that the mRNA expressions of numerous genes are regulated by nitrogen supplies; however, little is known about the expressed non-coding elements, for example long non-coding RNAs (lncRNAs) that control the response of maize (Zea mays L.) to nitrogen. LncRNAs are a class of non-coding RNAs larger than 200 bp, which have emerged as key regulators in gene expression. RESULTS: In this study, we surveyed the intergenic/intronic lncRNAs in maize B73 leaves at the V7 stage under conditions of N-deficiency and N-sufficiency using ribosomal RNA depletion and ultra-deep total RNA sequencing approaches. By integration with mRNA expression profiles and physiological evaluations, 7245 lncRNAs and 637 nitrogen-responsive lncRNAs were identified that exhibited unique expression patterns. Co-expression network analysis showed that the nitrogen-responsive lncRNAs were enriched mainly in one of the three co-expressed modules. The genes in the enriched module are mainly involved in NADH dehydrogenase activity, oxidative phosphorylation and the nitrogen compounds metabolic process. CONCLUSIONS: We identified a large number of lncRNAs in maize and illustrated their potential regulatory roles in response to N stress. The results lay the foundation for further in-depth understanding of the molecular mechanisms of lncRNAs' role in response to nitrogen stresses.

Comparative Proteomic Analysis of Cotton Fiber Development and Protein Extraction Method Comparison in Late Stage Fibers.

The distinct stages of cotton fiber development and maturation serve as a single-celled model for studying the molecular mechanisms of plant cell elongation, cell wall development and cellulose biosynthesis. However, this model system of plant cell development is compromised for proteomic studies due to a lack of an efficient protein extraction method during the later stages of fiber development, because of a recalcitrant cell wall and the presence of abundant phenolic compounds. Here, we compared the quality and quantities of proteins extracted from 25 dpa (days post anthesis) fiber with multiple protein extraction methods and present a comprehensive quantitative proteomic study of fiber development from 10 dpa to 25 dpa. Comparative analysis using a label-free quantification method revealed 287 differentially-expressed proteins in the 10 dpa to 25 dpa fiber developmental period. Proteins involved in cell wall metabolism and regulation, cytoskeleton development and carbohydrate metabolism among other functional categories in four fiber developmental stages were identified. Our studies provide protocols for protein extraction from maturing fiber tissues for mass spectrometry analysis and expand knowledge of the proteomic profile of cotton fiber development.

Global analysis of lysine acetylation suggests the involvement of protein acetylation in diverse biological processes in rice (Oryza sativa).

Lysine acetylation is a reversible, dynamic protein modification regulated by lysine acetyltransferases and deacetylases. Recent advances in high-throughput proteomics have greatly contributed to the success of global analysis of lysine acetylation. A large number of proteins of diverse biological functions have been shown to be acetylated in several reports in human cells, E.coli, and dicot plants. However, the extent of lysine acetylation in non-histone proteins remains largely unknown in monocots, particularly in the cereal crops. Here we report the mass spectrometric examination of lysine acetylation in rice (Oryza sativa). We identified 60 lysine acetylated sites on 44 proteins of diverse biological functions. Immunoblot studies further validated the presence of a large number of acetylated non-histone proteins. Examination of the amino acid composition revealed substantial amino acid bias around the acetylation sites and the amino acid preference is conserved among different organisms. Gene ontology analysis demonstrates that lysine acetylation occurs in diverse cytoplasmic, chloroplast and mitochondrial proteins in addition to the histone modifications. Our results suggest that lysine acetylation might constitute a regulatory mechanism for many proteins, including both histones and non-histone proteins of diverse biological functions.

Nuclear proteome response to cell wall removal in rice (Oryza sativa).

Plant cells are routinely exposed to various pathogens and environmental stresses that cause cell wall perturbations. Little is known of the mechanisms that plant cells use to sense these disturbances and transduce corresponding signals to regulate cellular responses to maintain cell wall integrity. Previous studies in rice have shown that removal of the cell wall leads to substantial chromatin reorganization and histone modification changes concomitant with cell wall re-synthesis. But the genes and proteins that regulate these cellular responses are still largely unknown. Here we present an examination of the nuclear proteome differential expression in response to removal of the cell wall in rice suspension cells using multiple nuclear proteome extraction methods. A total of 382 nuclear proteins were identified with two or more peptides, including 26 transcription factors. Upon removal of the cell wall, 142 nuclear proteins were up regulated and 112 were down regulated. The differentially expressed proteins included transcription factors, histones, histone domain containing proteins, and histone modification enzymes. Gene ontology analysis of the differentially expressed proteins indicates that chromatin & nucleosome assembly, protein-DNA complex assembly, and DNA packaging are tightly associated with cell wall removal. Our results indicate that removal of the cell wall imposes a tremendous challenge to the cells. Consequently, plant cells respond to the removal of the cell wall in the nucleus at every level of the regulatory hierarchy.

Polycomb group gene OsFIE2 regulates rice (Oryza sativa) seed development and grain filling via a mechanism distinct from Arabidopsis.

Cereal endosperm represents 60% of the calories consumed by human beings worldwide. In addition, cereals also serve as the primary feedstock for livestock. However, the regulatory mechanism of cereal endosperm and seed development is largely unknown. Polycomb complex has been shown to play a key role in the regulation of endosperm development in Arabidopsis, but its role in cereal endosperm development remains obscure. Additionally, the enzyme activities of the polycomb complexes have not been demonstrated in plants. Here we purified the rice OsFIE2-polycomb complex using tandem affinity purification and demonstrated its specific H3 methyltransferase activity. We found that the OsFIE2 gene product was responsible for H3K27me3 production specifically in vivo. Genetic studies showed that a reduction of OsFIE2 expression led to smaller seeds, partially filled seeds, and partial loss of seed dormancy. Gene expression and proteomics analyses found that the starch synthesis rate limiting step enzyme and multiple storage proteins are down-regulated in OsFIE2 reduction lines. Genome wide ChIP-Seq data analysis shows that H3K27me3 is associated with many genes in the young seeds. The H3K27me3 modification and gene expression in a key helix-loop-helix transcription factor is shown to be regulated by OsFIE2. Our results suggest that OsFIE2-polycomb complex positively regulates rice endosperm development and grain filling via a mechanism highly different from that in Arabidopsis.

Comparison of four ChIP-Seq analytical algorithms using rice endosperm H3K27 trimethylation profiling data.

Chromatin immunoprecipitation coupled with high throughput DNA Sequencing (ChIP-Seq) has emerged as a powerful tool for genome wide profiling of the binding sites of proteins associated with DNA such as histones and transcription factors. However, no peak calling program has gained consensus acceptance by the scientific community as the preferred tool for ChIP-Seq data analysis. Analyzing the large data sets generated by ChIP-Seq studies remains highly challenging for most molecular biology laboratories.Here we profile H3K27me3 enrichment sites in rice young endosperm using the ChIP-Seq approach and analyze the data using four peak calling algorithms (FindPeaks, PeakSeq, USeq, and MACS). Comparison of the four algorithms reveals that these programs produce very different peaks in terms of peak size, number, and position relative to genes. We verify the peak predictions using ChIP-PCR to evaluate the accuracy of peak prediction of the four algorithms. We discuss the approach of each algorithm and compare similarities and differences in the results. Despite their differences in the peaks identified, all of the programs reach similar conclusions about the effect of H3K27me3 on gene expression. Its presence either upstream or downstream of a gene is predominately associated with repression of the gene. Additionally, GO analysis finds that a substantially higher ratio of genes associated with H3K27me3 were involved in multicellular organism development, signal transduction, response to external and endogenous stimuli, and secondary metabolic pathways than the rest of the rice genome.

Transcriptional dynamics during cell wall removal and regeneration reveals key genes involved in cell wall development in rice.

Efficient and cost-effective conversion of plant biomass to usable forms of energy requires a thorough understanding of cell wall biosynthesis, modification and degradation. To elucidate these processes, we assessed the expression dynamics during enzymatic removal and regeneration of rice cell walls in suspension cells over time. In total, 928 genes exhibited significant up-regulation during cell wall removal, whereas, 79 genes were up-regulated during cell wall regeneration. Both gene sets are enriched for kinases, transcription factors and genes predicted to be involved in cell wall-related functions. Integration of the gene expression datasets with a catalog of known and/or predicted biochemical pathways from rice, revealed metabolic and hormonal pathways involved in cell wall degradation and regeneration. Rice lines carrying Tos17 mutations in genes up-regulated during cell wall removal exhibit dwarf phenotypes. Many of the genes up-regulated during cell wall development are also up-regulated in response to infection and environmental perturbations indicating a coordinated response to diverse types of stress.