Structure-based prediction of protein-protein interaction network in rice.

Comprehensive protein-protein interaction (PPI) maps are critical for understanding the functional organization of the proteome, but challenging to produce experimentally. Here, we developed a computational method for predicting PPIs based on protein docking. Evaluation of performance on benchmark sets demonstrated the ability of the docking-based method to accurately identify PPIs using predicted protein structures. By employing the docking-based method, we constructed a structurally resolved PPI network consisting of 24,653 interactions between 2,131 proteins, which greatly extends the current knowledge on the rice protein-protein interactome. Moreover, we mapped the trait-associated single nucleotide polymorphisms (SNPs) to the structural interactome, and computationally identified 14 SNPs that had significant consequences on PPI network. The protein structural interactome map provided a resource to facilitate functional investigation of PPI-perturbing alleles associated with agronomically important traits in rice.

Natural variation of Alfin-like family affects seed size and drought tolerance in rice.

The Alfin-like (AL) family is a group of small plant-specific transcriptional factors involved in abiotic stresses in dicotyledon. In an early study, we found an AL gene in rice that was associated with grain yield under drought stress. However, little information is known about the AL family in rice. In this study, AL genes in the rice genome were identified, and the OsAL proteins were found to locate in the nucleus and have no transcriptional self-activation activity. The expression of the OsALs was regulated by different environmental stimulations and plant hormones. Association and domestication analysis revealed that natural variation of most OsALs was significantly associated with yield traits, drought resistance and divergence in grain size in indica and japonica rice varieties. Hap1 of OsAL7.1 and Hap7 of OsAL11 were favorable haplotypes of seed weight and germination under osmotic stress. Furthermore, osal7.1 and osal11 mutants have larger seeds and are more sensitive to abscisic acid and mannitol during germination stage. Overexpressing of OsAL7.1 and OsAL11 in rice weakened the tolerance to drought in the adult stage. Thus, our work provides informative knowledge for exploring and harnessing haplotype diversity of OsALs to improve yield stability under drought stress.

Coordination of growth and drought responses by GA-ABA signaling in rice.

The drought caused by global warming seriously affects the crop growth and agricultural production. Plants have evolved distinct strategies to cope with the drought environment. Under drought stress, energy and resources should be diverted from growth toward stress management. However, the molecular mechanism underlying coordination of growth and drought response remains largely elusive. Here, we discovered that most of the gibberellin (GA) metabolic genes were regulated by water scarcity in rice, leading to the lower GA contents and hence inhibited plant growth. Low GA contents resulted in the accumulation of more GA signaling negative regulator SLENDER RICE 1, which inhibited the degradation of abscisic acid (ABA) receptor PYL10 by competitively binding to the co-activator of anaphase-promoting complex TAD1, resulting in the enhanced ABA response and drought tolerance. These results elucidate the synergistic regulation of crop growth inhibition and promotion of drought tolerance and survival, and provide useful genetic resource in breeding improvement of crop drought resistance.

A novel rice grain size gene OsSNB was identified by genome-wide association study in natural population.

Increasing agricultural productivity is one of the most important goals of plant science research and imperative to meet the needs of a rapidly growing population. Rice (Oryza sativa L.) is one of the most important staple crops worldwide. Grain size is both a major determinant of grain yield in rice and a target trait for domestication and artificial breeding. Here, a genome-wide association study of grain length and grain width was performed using 996,722 SNP markers in 270 rice accessions. Five and four quantitative trait loci were identified for grain length and grain width, respectively. In particular, the novel grain size gene OsSNB was identified from qGW7, and further results showed that OsSNB negatively regulated grain size. Most notably, knockout mutant plants by CRISPR/Cas9 technology showed increased grain length, width, and weight, while overexpression of OsSNB yielded the opposite. Sequencing of this gene from the promoter to the 3'-untranslated region in 168 rice accessions from a wide geographic range identified eight haplotypes. Furthermore, Hap 3 has the highest grain width discovered in japonica subspecies. Compared to other haplotypes, Hap 3 has a 225 bp insertion in the promoter. Based on the difference between Hap 3 and other haplotypes, OsSNB_Indel2 was designed as a functional marker for the improvement of rice grain width. This could be directly used to assist selection toward an improvement of grain width. These findings suggest OsSNB as useful for further improvements in yield characteristics in most cultivars.

Temporal responses of conserved miRNAs to drought and their associations with drought tolerance and productivity in rice.

BACKGROUND: Plant miRNAs play crucial roles in responses to drought and developmental processes. It is essential to understand the association of miRNAs with drought-tolerance (DT), as well as their impacts on growth, development, and reproduction (GDP). This will facilitate our utilization of rice miRNAs in breeding. RESULTS: In this study, we investigated the time course of miRNA responses to a long-term drought among six rice genotypes by high-throughput sequencing. In total, 354 conserved miRNAs were drought responsive, representing obvious genotype- and stage-dependent patterns. The drought-responsive miRNAs (DRMs) formed complex regulatory network via their coexpression and direct/indirect impacts on the rice transcriptome. Based on correlation analyses, 211 DRMs were predicted to be associated with DT and/or GDP. Noticeably, 14.2% DRMs were inversely correlated with DT and GDP. In addition, 9 pairs of mature miRNAs, each derived from the same pre-miRNAs, were predicted to have opposite roles in regulating DT and GDP. This suggests a potential yield penalty if an inappropriate miRNA/pre-miRNA is utilized. miRNAs have profound impacts on the rice transcriptome reflected by great number of correlated drought-responsive genes. By regulating these genes, a miRNA could activate diverse biological processes and metabolic pathways to adapt to drought and have an influence on its GDP. CONCLUSION: Based on the temporal pattern of miRNAs in response to drought, we have described the complex network between DRMs. Potential associations of DRMs with DT and/or GDP were disclosed. This knowledge provides valuable information for a better understanding in the roles of miRNAs play in rice DT and/or GDP, which can facilitate our utilization of miRNA in breeding.

An APETALA2/ethylene responsive factor, OsEBP89 knockout enhances adaptation to direct-seeding on wet land and tolerance to drought stress in rice.

Water stress is the most important adverse factor limiting rice production. Too much water leads to flood and too little leads to drought. Floods and droughts can severely damage crop at different times of the rice life cycle. So the research on submergence tolerance and drought resistance of rice is particularly urgent. In this study, we reported that OsEBP89 (Oryza sativa Ethylene-responsive element binding protein, clone 89), a member of the AP2/ERF subfamily, is involved in a novel signal transduction associated with the tolerance to drought and submergence stress. OsEBP89 was found to be strongly inhibited by drought stress and promoted by submergence. The OsEBP89 protein was located at the nucleus in the rice protoplast. Loss of OsEBP89 was found to improve the seed germination under submerged conditions and also enhanced the tolerance to drought stress throughout growth stage. Additionally, OsEBP89 knockout rice plants increased the accumulation of proline, improved the ability to scavenge ROS compared to overexpression lines and wild type after PEG treatment. Transcriptome data indicates that knockout of OsEBP89 improved the expression of specific genes in response to adverse factors, such as OsAPX1, OsHsfA3, and OsP5CS. Further results indicate that OsEBP89 can interact with and be phosphorylated by SnRK1α (sucrose non-fermenting-1-related protein kinase-1 gene). These findings provide insight into the mechanism of abiotic stress tolerance, and suggest OsEBP89 as a new genetic engineering resource to improve abiotic stress tolerance in rice.

Drought stress modulates diurnal oscillations of circadian clock and drought-responsive genes in Oryza sativa L.

Endogenous circadian rhythms play a key role in regulating plant growth and development, and in allowing plants to respond and adapt to changing environments. To understand how drought regulates upland rice(Oryza sativa L.) IRAT109, we examined the expression levels of circadian clock and drought-responsive genes through real-time PCR. The results revealed that, first, drought reduced the relative expression level and amplitude of peak expression of several morning circadian clock components (such as OsPRRs, OsLHY and OsZTL1), increased the relative expression level and amplitude of some evening circadian clock components (such as OsTOC1, OsGI and OsELF3), but did not influence OsFKF1. Secondly, the relative expression level of most drought-responsive genes was generally increased, except for OsDST, a negative regulator. Lastly, expression rhythms of most drought-responsive genes were disturbed, but not that of OsCIPK12, OsCDPK7 and OsDREB1A. The results indicate that drought stress modulates the expression of circadian clock components and the interplay regulates diurnal oscillations of relative genes.

Autotetraploidy of rice does not potentiate the tolerance to drought stress in the seedling stage.

Polyploid is considered an advantage that has evolved to be more environmentally adaptable than its diploid. To understand if doubled chromosome of diploid rice can improve drought tolerance, we evaluated the diploid (2X) and autotetraploid (4X) plants of three indica and three japonica varieties. Drought stress in the plastic bucket of four-leaf stage revealed that the drought tolerance of 4X plants was lower than that of its diploid donor plants. The assay of photosynthetic rate of all varieties showed that all 4X varieties had lower rates than their diploid donors. The capacity for reactive oxygen species production and scavenging varied among different 2X and 4X varieties. Further, transcriptomic analysis of 2X and 4X plants of four varieties under normal and drought condition showed that the wide variation of gene expression was caused by difference of varieties, not by chromosome ploidy. However, weighted gene co-expression network analysis (WGCNA) revealed that the severe interference of photosynthesis-related genes in tetraploid plants under drought stress is the primary reason for the decrease of drought tolerance in autotetraploid lines. Consistently, new transcripts analysis in autotetraploid revealed that the gene transcription related with mitochondrion and plastid of cell component was influenced most significantly. The results indicated that chromosome doubling of diploid rice weakened their drought tolerance, primarily due to disorder of photosynthesis-related genes in tetraploid plants under drought stress. Maintain tetraploid drought tolerance through chromosome doubling breeding in rice needs to start with the selection of parental varieties and more efforts.

Overexpressing PpBURP2 in Rice Increases Plant Defense to Abiotic Stress and Bacterial Leaf Blight.

Mosses are one of the earliest diverging land plants that adapted to living on land. The BURP domain-containing proteins (BURP proteins) are plant-specific proteins that appeared when plants shifted from aquatic environments to land. Phylogenetic analysis revealed that the BURP domain of higher plants is originated from lower land plants and divergent because of motif conversion. To discover the function of BURP protein in moss, rice transgenics with ectopic expression of PpBURP2 were subjected to different abiotic stresses treatments. The results revealed that the ectopic expression of PpBURP2 enhanced the tolerance to osmotic and saline stresses at the seedling stage and drought stress at the adult stage. Further ectopic expression of PpBURP2 improved the cadmium (2+) (Cd2+) tolerance and reduced Cd2+ accumulation in rice leaves. Transcriptomic analysis of the transgenic PpBURP2 plants showed that the differentially expressed genes were involved in the metabolism of secondary metabolites, energy, oxidation-reduction process, and defense-related genes. Further experiments showed that the photosynthetic efficiency and resistance against bacterial leaf blight were obviously improved in transgenic plants. Yeast two-hybrid and bimolecular fluorescence complementation (BiFC) assays revealed the physical interaction of BURP domain protein from rice and moss with mitogen-activated protein kinase kinase (MKK) from rice. Therefore, our findings demonstrate that overexpressing PpBURP2 in rice confers resistance to abiotic stresses and bacterial leaf blight. They also suggested that the regulatory role of BURP-like proteins across lower and higher plants was evolutionary conservation of responses of different classes of plants to different environmental challenges.

Overexpression of OsHMGB707, a High Mobility Group Protein, Enhances Rice Drought Tolerance by Promoting Stress-Related Gene Expression.

Drought stress adversely affects crop growth and productivity worldwide. In response, plants have evolved several strategies in which numerous genes are induced to counter stress. High mobility group (HMG) proteins are the second most abundant family of chromosomal proteins. They play a crucial role in gene transcriptional regulation by modulating the chromatin/DNA structure. In this study, we isolated a novel HMG gene, OsHMGB707, one of the candidate genes localized in the quantitative trait loci (QTL) interval of rice drought tolerance, and examined its function on rice stress tolerance. The expression of OsHMGB707 was up-regulated by dehydration and high salt treatment. Its overexpression significantly enhanced drought tolerance in transgenic rice plants, whereas its knockdown through RNA interference (RNAi) did not affect the drought tolerance of the transgenic rice plants. Notably, OsHMGB707-GFP is localized in the cell nucleus, and OsHMGB707 is protein-bound to the synthetic four-way junction DNA. Several genes were up-regulated in OsHMGB707-overexpression (OE) rice lines compared to the wild-type rice varieties. Some of the genes encode stress-related proteins (e.g., DREB transcription factors, heat shock protein 20, and heat shock protein DnaJ). In summary, OsHMGB707 encodes a stress-responsive high mobility group protein and regulates rice drought tolerance by promoting the expression of stress-related genes.

Isolation of a novel lipase gene from Serratia liquefaciens S33 DB-1, functional expression in Pichia pastoris and its properties.

A new lipase gene designated as SlLipA was isolated from Serratia liquefaciens S33 DB-1 by the genomic-walking method. The cloned gene contained an open reading frame (ORF) of 1,845 bp encoding 615 amino acids with a conserved GXSXG motif. Genome sequence analysis showed that an aldo/keto reductase gene closed to the SlLipA gene. The lipase gene was cloned into the expression vector pPICZalphaA and successfully integrated into the heterologous host, methylotrophic yeast Pichia pastoris GS115. Five transformants could be expressed as secreted recombinant proteins with the high activity on Triglyceride-Agarose plate and as candidates to produce the recombinant enzyme. A C-terminal His tag was used for its purification. The lipase activity of different transformants against substrate para-nitrophenyl laurate (p-NPL) varied from 14 to 16 U ml(-1). For the substrates para-nitrophenyl caprate (p-NPC), p-NPL, para-nitrophenyl myristate (p-NPM), para-nitrophenyl palmitate (p-NPP), and para-nitrophenyl stearate (p-NPS), the specific activity was shown to be preferred to long acyl chain length of p-NPS.

Expression analysis of a novel pyridoxal kinase messenger RNA splice variant, PKL, in oil rape suffering abiotic stress and phytohormones.

Pyridoxal kinase is key enzyme for the biosynthesis of pyridoxal 5'-phosphate, the biologically active form of vitamin B6, in the salvage pathway. A pyridoxal kinase gene, BnPKL (GenBank accession No. DQ463962), was isolated from oilseed rape (Brassica napus L.) following water stress through rapid amplification of complementary DNA (cDNA) ends. The results showed that the gene had two splice variants: PKL and PKL2. PKL, the long cDNA, encodes a 334 amino acid protein with a complete ATP-binding site, pyridoxal kinase-binding site and dimer interface site of a pyridoxal kinase, while PKL2, the short cDNA, lacked a partial domain. Southern blot showed that there were two copies in Brassica napus. The expression of BnPKL cDNA could rescue the mutant phenotype of Escherichia coli defective in pyridoxal kinase. Real-time reverse transcription-polymerase chain reaction revealed that the relative abundance of two transcripts are modulated by development and environmental stresses. Abscisic acid and NaCl were inclined to decrease PKL expression, but H2O2 and cold temperatures induced the PKL expression. In addition, the PKL expression could be transiently induced by jasmonate acid at an early stage, abscisic acid, salicylic acid and jasmonate acid enhanced the PKL expression in roots. Our results demonstrated that BnPKL was a pyridoxal kinase involved in responses to biotic and abiotic stresses.

Comparative Analysis of Expression Profiles of Panicle Development among Tolerant and Sensitive Rice in Response to Drought Stress.

Water deficit caused a serious threat to crops, especially panicle development at reproductive growth phase. We investigated grain yield components and gene expression profiles of panicle among tolerant and sensitive rice in response to drought stress. Panicle morphologies exhibited that secondary branches per panicle were more severely affected as compared to primary branches per panicle. Moreover, grain weight per panicle showed significant decrease for both tolerant and sensitive varieties except for MILT1444. Expression profile analysis revealed that 783 differentially expressed genes (DEGs) were identified to be drought-induced from young panicles in 2 cm length. Hierarchical clustering indicated that 76.8% of DEGs were up-regulated for all six rice varieties, and the percentage of down-regulated genes was higher in sensitive group than tolerant group. Biological process category revealed that the shared Gene Ontology (GO) terms were involved in response to abiotic stimulus and stress, whereas the specific GO terms in tolerant group were identified as regulation of biological quality, homeostatic process, cell growth, anatomical structure morphogenesis and development, and the unique terms in sensitive varieties were identified as lipid metabolic process and secondary metabolic process. Furthermore, the gene-based association analysis narrowed down list of DEGs, and four genes common to all six varieties were selected as candidate for breeders. Together, we found several shared and distinct biological processes between tolerant and sensitive varieties, and candidate stress-responsive genes. These findings provided insight into functional mechanisms regulating drought stress response in panicle development and may also help to crop tolerant improvement.

A novel gene OsAHL1 improves both drought avoidance and drought tolerance in rice.

A novel gene, OsAHL1, containing an AT-hook motif and a PPC domain was identified through genome-wide profiling and analysis of mRNAs by comparing the microarray of drought-challenged versus normally watered rice. The results indicated OsAHL1 has both drought avoidance and drought tolerance that could greatly improve drought resistance of the rice plant. Overexpression of OsAHL1 enhanced multiple stress tolerances in rice plants during both seedling and panicle development stages. Functional studies revealed that OsAHL1 regulates root development under drought condition to enhance drought avoidance, participates in oxidative stress response and also regulates the content of chlorophyll in rice leaves. OsAHL1 specifically binds to the A/T rich sequence region of promoters or introns, and hence directly regulates the expression of many stress related downstream genes.

Cloning and expression analysis of a water stress-induced gene from Brassica oleracea.

Plant growth and productivity are greatly affected by water stress, such as drought and salinity. Here we report on the cloning and expression analysis of a water stress-induced gene from Brassica oleracea (designated as BoWS, GenBank accession number AY571333) by rapid amplification of cDNA ends (RACE). The full-length cDNA of BoWS consisted of 594 bp and contained a 285 bp open reading frame (ORF) encoding a 95-amino-acid protein. The deduced protein had a calculated molecular mass of 10.53 kDa and an isoelectric point of 6.93. The sequence similarity and comparative analysis showed that BoWS was 84% identical to Arabidopsis thaliana putative water stress-induced protein (GenBank accession number AAM67282). Southern blot analysis indicated that BoWS was a low-copy gene. Northern blot analysis revealed that the expression of BoWS was upregulated by abscisic acid (ABA), mannitol, NaCl, drought, salicylic acid (SA) and hydrogen peroxide (H(2)O(2)). Our results indicate that BoWS is extremely related to the water-deficit stress in B. oleracea.

Identification of gene modules associated with drought response in rice by network-based analysis.

Understanding the molecular mechanisms that underlie plant responses to drought stress is challenging due to the complex interplay of numerous different genes. Here, we used network-based gene clustering to uncover the relationships between drought-responsive genes from large microarray datasets. We identified 2,607 rice genes that showed significant changes in gene expression under drought stress; 1,392 genes were highly intercorrelated to form 15 gene modules. These drought-responsive gene modules are biologically plausible, with enrichments for genes in common functional categories, stress response changes, tissue-specific expression and transcription factor binding sites. We observed that a gene module (referred to as module 4) consisting of 134 genes was significantly associated with drought response in both drought-tolerant and drought-sensitive rice varieties. This module is enriched for genes involved in controlling the response of the plant to water and embryonic development, including a heat shock transcription factor as the key regulator in the expression of ABRE-containing genes. These results suggest that module 4 is highly conserved in the ABA-mediated drought response pathway in different rice varieties. Moreover, our study showed that many hub genes clustered in rice chromosomes had significant associations with QTLs for drought stress tolerance. The relationship between hub gene clusters and drought tolerance QTLs may provide a key to understand the genetic basis of drought tolerance in rice.

Identification of rice transcription factors associated with drought tolerance using the Ecotilling method.

The drought tolerance (DT) of plants is a complex quantitative trait. Under natural and artificial selection, drought tolerance represents the crop survival ability and production capacity under drought conditions (Luo, 2010). To understand the regulation mechanism of varied drought tolerance among rice genotypes, 95 diverse rice landraces or varieties were evaluated within a field screen facility based on the 'line-source soil moisture gradient', and their resistance varied from extremely resistant to sensitive. The method of Ecotype Targeting Induced Local Lesions in Genomes (Ecotilling) was used to analyze the diversity in the promoters of 24 transcription factor families. The bands separated by electrophoresis using Ecotilling were converted into molecular markers. STRUCTURE analysis revealed a value of K = 2, namely, the population with two subgroups (i.e., indica and japonica), which coincided very well with the UPGMA clusters (NTSYS-pc software) using distance-based analysis and InDel markers. Then the association analysis between the promoter diversity of these transcription factors and the DT index/level of each variety was performed. The results showed that three genes were associated with the DT index and that five genes were associated with the DT level. The sequences of these associated genes are complex and variable, especially at approximately 1000 bp upstream of the transcription initiation sites. The study illuminated that association analysis aimed at Ecotilling diversity of natural groups could facilitate the isolation of rice genes related to complex quantitative traits.

Distributional gradient of amino acid repeats in plant proteins.

A computer-based analysis was conducted to assess the characteristics of amino acid repeats in Arabidopsis and rice. Our analysis showed a negative gradient in amino acid repeat distribution along the direction of translation in plants. Repeat occurrences are obviously associated with position in plant proteins but are not consistent with the corresponding amino acid contents. These repeats are encoded by the mixed synonymous codons rather than the uninterrupted reiterations of a single codon, and both Arabidopsis and rice have gradients in their distribution. Functional investigation showed that these repeat-containing proteins are preferentially involved in transcription regulation and protein ubiquitination but significantly underrepresented in the processes of DNA recombination and DNA replication. These data reveal that the direction-related mutation bias and functional selection have influenced the distribution of amino acid repeats in plants.

Preference of simple sequence repeats in coding and non-coding regions of Arabidopsis thaliana.

MOTIVATION: Simple sequence repeats or microsatellites have been found abundantly in many genomes. However, the significance of distribution preference has not been completely understood. Completion of the Arabidopsis genome sequencing allows us to better understand and characterize microsatellites. RESULTS: Microsatellite distribution was more abundant in 5'-flanking regions of genes compared with that expected in the whole genome, with an over-representation of AG and AAG repeats; there were clear differences from distributions in 3'-flanks and coding fractions, where triplet frequencies evidently corresponded to codon usage. We identified 1140 full-length genes that contained at least one locus of AG or AAG repeats in their upstream sequences, and whose functional characteristics were significantly associated with the repeats. This observation indicates that selective pressure markedly differed in the three transcribed regions, with positive selection of AG and AAG repeats in 5'-flanks close to those genes whose products are preferentially involved in transcription.