Functional analysis of the ABA-responsive protein family in ABA and stress signal transduction in Arabidopsis

The melon Fom-1-Prv resistance gene pair: Correlated spatial expression and interaction with a viral protein

The head-to-head oriented pair of melon resistance genes, Fom-1 and Prv, control resistance to Fusarium oxysporum races 0 and 2 and papaya ringspot virus (PRSV), respectively. They encode, via several RNA splice variants, TIR-NBS-LRR proteins, and Prv has a C-terminal extra domain with a second NBS homologous sequence. In other systems, paired R-proteins were shown to operate by "labor division," with one protein having an extra integrated domain that directly binds the pathogen's Avr factor, and the second protein executing the defense response. We report that the expression of the two genes in two pairs of near-isogenic lines was higher in the resistant isoline and inducible by F. oxysporum race 2 but not by PRSV. The intergenic DNA region separating the coding sequences of the two genes acted as a bi-directional promoter and drove GUS expression in transgenic melon roots and transgenic tobacco plants. Expression of both genes was strong in melon root tips, around the root vascular cylinder, and the phloem and xylem parenchyma of tobacco stems and petioles. The pattern of GUS expression suggests coordinated expression of the two genes. In agreement with the above model, Prv's extra domain was shown to interact with the cylindrical inclusion protein of PRSV both in yeast cells and in planta.

A novel insight into transcriptional and epigenetic regulation underlying sex expression and flower development in melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important cucurbit and has been considered as a model plant for studying sex determination. The four most common sexual morphotypes in melon are monoecious (A-G-M), gynoecious (--ggM-), andromonoecious (A-G-mm), and hermaphrodite (--ggmm). Sex expression in melons is complex, as the genes and associated networks that govern the sex expression are not fully explored. Recently, RNA-seq transcriptomic profiling, ChIP-qPCR analysis integrated with gene ontology annotation and Kyoto Encyclopedia of Genes and Genomes pathways predicted the differentially expressed genes including sex-specific ACS and ACO genes, in regulating the sex-expression, phytohormonal cross-talk, signal transduction, and secondary metabolism in melons. Integration of transcriptional control through genetic interaction in between the ACS7, ACS11, and WIP1 in epistatic or hypostatic manner, along with the recruitment of H3K9ac and H3K27me3, epigenetically, overall determine sex expression. Alignment of protein sequences for establishing phylogenetic evolution, motif comparison, and protein-protein interaction supported the structural conservation while presence of the conserved hydrophilic and charged residues across the diverged evolutionary group predicted the functional conservation of the ACS protein. Presence of the putative cis-binding elements or DNA motifs, and its further comparison with DAP-seq-based cistrome and epicistrome of Arabidopsis, unraveled strong ancestry of melons with Arabidopsis. Motif comparison analysis also characterized putative genes and transcription factors involved in ethylene biosynthesis, signal transduction, and hormonal cross-talk related to sex expression. Overall, we have comprehensively reviewed research findings for a deeper insight into transcriptional and epigenetic regulation of sex expression and flower development in melons.

Pentatricopeptide repeat protein PHOTOSYSTEM I BIOGENESIS FACTOR2 is required for splicing of ycf3

To gain a better understanding of the molecular mechanisms of photosystem I (PSI) biogenesis, we characterized the Arabidopsis thaliana photosystem I biogenesis factor 2 (pbf2) mutant, which lacks PSI complex. PBF2 encodes a P-class pentatricopeptide repeat (PPR) protein. In the pbf2 mutants, we observed a striking decrease in the transcript level of only one gene, the chloroplast gene ycf3, which is essential for PSI assembly. Further analysis of ycf3 transcripts showed that PBF2 is specifically required for the splicing of ycf3 intron 1. Computational prediction of binding sequences and electrophoretic mobility shift assays reveal that PBF2 specifically binds to a sequence in ycf3 intron 1. Moreover, we found that PBF2 interacted with two general factors for group II intron splicing CHLOROPLAST RNA SPLICING2-ASSOCIATED FACTOR1 (CAF1) and CAF2, and facilitated the association of these two factors with ycf3 intron 1. Our results suggest that PBF2 is specifically required for the splicing of ycf3 intron 1 through cooperating with CAF1 and CAF2. Our results also suggest that additional proteins are required to contribute to the specificity of CAF-dependent group II intron splicing.

Global analysis of expression profile of members of DnaJ gene families involved in capsaicinoids synthesis in pepper (Capsicum annuum L)

BACKGROUND

The DnaJ proteins play critical roles in plant development and stress responses. Recently, seventy-six DnaJ genes were identified through a comprehensive bioinformatics analysis in the pepper genome. However, there were no reports on understanding of phylogenetic relationships and diverse expression profile of pepper DnaJ genes to date. Herein, we performed the systemic analysis of the phylogenetic relationships and expression profile of pepper DnaJ genes in different tissues and in response to both abiotic stress and plant hormones.

RESULTS

Phylogenetic analysis showed that all the pepper DnaJ genes were grouped into 7 sub-families (sub-family I, II, III, IV, V, VI and VII) according to sequence homology. The expression of pepper DnaJs in different tissues revealed that about 38% (29/76) of pepper DnaJs were expressed in at least one tissue. The results demonstrate the potentially critical role of DnaJs in pepper growth and development. In addition, to gain insight into the expression difference of pepper DnaJ genes in placenta between pungent and non-pungent, their expression patterns were also analyzed using RNA-seq data and qRT-PCR. Comparison analysis revealed that eight genes presented distinct expression profiles in pungent and non-pungent pepper. The CaDnaJs co-expressed with genes involved in capsaicinoids synthesis during placenta development. What is more, our study exposed the fact that these eight DnaJ genes were probably regulated by stress (heat, drought and salt), and were also regulated by plant hormones (ABA, GA3, MeJA and SA).

CONCLUSIONS

In summary, these results showed that some DnaJ genes expressed in placenta may be involved in plant response to abiotic stress during biosynthesis of compounds related with pungency. The study provides wide insights to the expression profiles of pepper DanJ genes and contributes to our knowledge about the function of DnaJ genes in pepper.

OsABAR1, a novel GRAM domain-containing protein, confers drought and salt tolerance via an ABA-dependent pathway in rice

Glucosyltransferases-like GTPase activators and Myotubularin (GRAM) domain-containing proteins are important for plant development and responses to biotic stresses. However, the effects of GRAM proteins on abiotic stress responses remain unclear. In this study, we identified a novel GRAM protein-encoding gene, OsABAR1, and characterized its regulatory functions related to rice drought and salt tolerance. The OsABAR1 protein was localized in the cytoplasm and nucleus. Among all examined organs, the OsABAR1 transcript level was highest in the roots. Moreover, OsABAR1 expression was up-regulated by drought and salinity stresses. The OsABAR1-overexpressing (OsABAR1-OX) lines exhibited enhanced tolerance to drought and salinity, whereas the knock-out lines (Osabar1) had the opposite phenotypes. We further analyzed the involvement of OsABAR1 in the abscisic acid (ABA) signaling pathway. The OsABAR1 expression level was up-regulated by ABA. In turn, OsABAR1 regulated the expression of ABA metabolic genes and responsive genes. Furthermore, OsABAR1-OX seedlings were hypersensitive to exogenous ABA, whereas Osabar1 seedlings were hyposensitive. These results imply that OsABAR1 is a positive regulator of the ABA pathway and confirm that OsABAR1 improves rice drought and salt tolerance via an ABA-dependent pathway. This study is the first to clarify the regulatory roles of GRAM proteins in rice responses to abiotic stresses.

Genomic Organization and Comparative Phylogenic Analysis of NBS-LRR Resistance Gene Family in Solanum pimpinellifolium and Arabidopsis thaliana

NBS-LRR (nucleotide-binding site and leucine-rich repeat) is one of the largest resistance gene families in plants. The completion of the genome sequencing of wild tomato Solanum pimpinellifolium provided an opportunity to conduct a comprehensive analysis of the NBS-LRR gene superfamily at the genome-wide level. In this study, gene identification, chromosome mapping, and phylogenetic analysis of the NBS-LRR gene family were analyzed using the bioinformatics methods. The results revealed 245 NBS-LRRs in total, similar to that in the cultivated tomato. These genes are unevenly distributed on 12 chromosomes, and ~59.6% of them form gene clusters, most of which are tandem duplications. Phylogenetic analysis divided the NBS-LRRs into 2 subfamilies (CNL-coiled-coil NBS-LRR and TNL-TIR NBS-LRR), and the expansion of the CNL subfamily was more extensive than the TNL subfamily. Novel conserved structures were identified through conserved motif analysis between the CNL and TNL subfamilies. Compared with the NBS-LRR sequences from the model plant Arabidopsis thaliana, wide genetic variation occurred after the divergence of S. pimpinellifolium and A thaliana. Species-specific expansion was also found in the CNL subfamily in S. pimpinellifolium. The results of this study provide the basis for the deeper analysis of NBS-LRR resistance genes and contribute to mapping and isolation of candidate resistance genes in S. pimpinellifolium.

New Insights into Evolution of Plant Heat Shock Factors (Hsfs) and Expression Analysis of Tea Genes in Response to Abiotic Stresses

Heat shock transcription factor (Hsf) is one of key regulators in plant abotic stress response. Although the Hsf gene family has been identified from several plant species, original and evolution relationship have been fragmented. In addition, tea, an important crop, genome sequences have been completed and function of the Hsf family genes in response to abiotic stresses was not illuminated. In this study, a total of 4208 Hsf proteins were identified within 163 plant species from green algae (Gonium pectorale) to angiosperm (monocots and dicots), which were distributed unevenly into each of plant species tested. The result indicated that Hsf originated during the early evolutionary history of chlorophytae algae and genome-wide genetic varies had occurred during the course of evolution in plant species. Phylogenetic classification of Hsf genes from the representative nine plant species into ten subfamilies, each of which contained members from different plant species, imply that gene duplication had occurred during the course of evolution. In addition, based on RNA-seq data, the member of the Hsfs showed different expression levels in the different organs and at the different developmental stages in tea. Expression patterns also showed clear differences among Camellia species, indicating that regulation of Hsf genes expression varied between organs in a species-specific manner. Furthermore, expression of most Hsfs in response to drought, cold and salt stresses, imply a possible positive regulatory role under abiotic stresses. Expression profiles of nineteen Hsf genes in response to heat stress were also analyzed by quantitative real-time RT-PCR. Several stress-responsive Hsf genes were highly regulated by heat stress treatment. In conclusion, these results lay a solid foundation for us to elucidate the evolutionary origin of plant Hsfs and Hsf functions in tea response to abiotic stresses in the future.

Metabolome and Transcriptome Analysis of Hexaploid Solidago canadensis Roots Reveals its Invasive Capacity Related to Polyploidy

Polyploid plants are more often invasive species than their diploid counterparts. As the invasiveness of a species is often linked to its production of allelopathic compounds, we hypothesize that differences in invasive ability between cytotypes may be due to their different ability to synthesize allelopathic metabolites. We test this using two cytotypes of Solidago canadensis as the model and use integrated metabolome and transcriptome data to resolve the question. Metabolome analysis identified 122 metabolites about flavonoids, phenylpropanoids and terpenoids, of which 57 were differentially accumulated between the two cytotypes. Transcriptome analysis showed that many differentially expressed genes (DEGs) were enriched in ‘biosynthesis of secondary metabolites’, ‘plant hormone signal transduction’, and ‘MAPK signaling’, covering most steps of plant allelopathic metabolite synthesis. Importantly, the differentially accumulated flavonoids, phenylpropanoids and terpenoids were closely correlated with related DEGs. Furthermore, 30 miRNAs were found to be negatively associated with putative targets, and they were thought to be involved in target gene expression regulation. These miRNAs probably play a vital role in the regulation of metabolite synthesis in hexaploid S. canadensis. The two cytotypes of S. canadensis differ in the allelopathic metabolite synthesis and this difference is associated with regulation of expression of a range of genes. These results suggest that changes in gene expression may underlying the increased invasive potential of the polyploidy.

Evolutionary Conservation and Expression Patterns of Neutral/Alkaline Invertases in Solanum

The invertase gene family in plants is composed of two subfamilies of enzymes, namely, acid- and neutral/alkaline invertases (cytosolic invertase, CIN). Both can irreversibly cleave sucrose into fructose and glucose, which are thought to play key roles in carbon metabolism and plant growth. CINs are widely found in plants, but little is reported about this family. In this paper, a comparative genomic approach was used to analyze the CIN gene family in Solanum, including Solanum tuberosum, Solanum lycopersicum, Solanum pennellii, Solanum pimpinellifolium, and Solanum melongena. A total of 40 CINs were identified in five Solanum plants, and sequence features, phylogenetic relationships, motif compositions, gene structure, collinear relationship, and expression profile were further analyzed. Sequence analysis revealed a remarkable conservation of CINs in sequence length, gene number, and molecular weight. The previously verified four amino acid residues (D188, E414, Arg430, and Ser547) were also observed in 39 out of 40 CINs in our study, showing to be deeply conserved. The CIN gene family could be distinguished into groups α and β, and α is further subdivided into subgroups α1 and α2 in our phylogenetic tree. More remarkably, each species has an average of four CINs in the α and β groups. Marked interspecies conservation and collinearity of CINs were also further revealed by chromosome mapping. Exon-intron configuration and conserved motifs were consistent in each of these α and β groups on the basis of in silico analysis. Expression analysis indicated that CINs were constitutively expressed and share similar expression profiles in all tested samples from S. tuberosum and S. lycopersicum. In addition, in CIN genes of the tomato and potato in response to abiotic and biotic stresses, phytohormones also performed. Overall, CINs in Solanum were encoded by a small and highly conserved gene family, possibly reflecting structural and functional conservation in Solanum. These results lay the foundation for further expounding the functional characterization of CIN genes and are also significant for understanding the evolutionary profiling of the CIN gene family in Solanum.

Genome-Wide Association Study Reveals Novel Genomic Regions for Grain Yield and Yield-Related Traits in Drought-Stressed Synthetic Hexaploid Wheat

Synthetic hexaploid wheat (SHW; 2n = 6x = 42, AABBDD, Triticum aestivum L.) is produced from an interspecific cross between durum wheat (2n = 4x = 28, AABB, T. turgidum L.) and goat grass (2n = 2x = 14, DD, Aegilops tauschii Coss.) and is reported to have significant novel alleles-controlling biotic and abiotic stresses resistance. A genome-wide association study (GWAS) was conducted to unravel these loci [marker⁻trait associations (MTAs)] using 35,648 genotyping-by-sequencing-derived single nucleotide polymorphisms in 123 SHWs. We identified 90 novel MTAs (45, 11, and 34 on the A, B, and D genomes, respectively) and haplotype blocks associated with grain yield and yield-related traits including root traits under drought stress. The phenotypic variance explained by the MTAs ranged from 1.1% to 32.3%. Most of the MTAs (120 out of 194) identified were found in genes, and of these 45 MTAs were in genes annotated as having a potential role in drought stress. This result provides further evidence for the reliability of MTAs identified. The large number of MTAs (53) identified especially on the D-genome demonstrate the potential of SHWs for elucidating the genetic architecture of complex traits and provide an opportunity for further improvement of wheat under rapidly changing climatic conditions.

Physiological and iTRAQ-Based Proteomic Analyses Reveal the Function of Spermidine on Improving Drought Tolerance in White Clover

Endogenous spermidine interacting with phytohormones may be involved in the regulation of differentially expressed proteins (DEPs) associated with drought tolerance in white clover. Plants treated with or without spermidine (50 μM) were subjected to 20% PEG 6000 nutrient solution to induce drought stress (50% leaf-relative water content). The results showed that increased endogenous spermidine induced by exogenous spermidine altered endogenous phytohormones in association with improved drought tolerance, as demonstrated by the delay in water-deficit development, improved photosynthesis and water use efficiency, and lower oxidative damage. As compared to untreated plants, Spd-treated plants maintained a higher abundance of DEPs under drought stress involved in (1) protein biosynthesis (ribosomal and chaperone proteins); (2) amino acids synthesis; (3) the carbon and energy metabolism; (4) antioxidant and stress defense (ascorbate peroxidase, glutathione peroxidase, and dehydrins); and (5) GA and ABA signaling pathways (gibberellin receptor GID1, ABA-responsive protein 17, and ABA stress ripening protein). Thus, the findings of proteome could explain the Spd-induced physiological effects associated with drought tolerance. The analysis of functional protein-protein networks further proved that the alteration of endogenous spermidine and phytohormones induced the interaction among ribosome, photosynthesis, carbon metabolism, and amino acid biosynthesis. These differences could contribute to improved drought tolerance.