Genome-scale identification of soybean BURP domain-containing genes and their expression under stress treatments

Plant peptides - redefining an area of ribosomally synthesized and post-translationally modified peptides

Covering 1965 to February 2024Plants are prolific peptide chemists and are known to make thousands of different peptidic molecules. These peptides vary dramatically in their size, chemistry, and bioactivity. Despite their differences, all plant peptides to date are biosynthesized as ribosomally synthesized and post-translationally modified peptides (RiPPs). Decades of research in plant RiPP biosynthesis have extended the definition and scope of RiPPs from microbial sources, establishing paradigms and discovering new families of biosynthetic enzymes. The discovery and elucidation of plant peptide pathways is challenging due to repurposing and evolution of housekeeping genes as both precursor peptides and biosynthetic enzymes and due to the low rates of gene clustering in plants. In this review, we highlight the chemistry, biosynthesis, and function of the known RiPP classes from plants and recommend a nomenclature for the recent addition of BURP-domain-derived RiPPs termed burpitides. Burpitides are an emerging family of cyclic plant RiPPs characterized by macrocyclic crosslinks between tyrosine or tryptophan side chains and other amino acid side chains or their peptide backbone that are formed by copper-dependent BURP-domain-containing proteins termed burpitide cyclases. Finally, we review the discovery of plant RiPPs through bioactivity-guided, structure-guided, and gene-guided approaches.

An intramolecular macrocyclase in plant ribosomal peptide biosynthesis

The biosynthetic dogma of ribosomally synthesized and posttranslationally modified peptides (RiPP) involves enzymatic intermolecular modification of core peptide motifs in precursor peptides. The plant-specific BURP-domain protein family, named after their four founding members, includes autocatalytic peptide cyclases involved in the biosynthesis of side-chain-macrocyclic plant RiPPs. Here we show that AhyBURP, a representative of the founding Unknown Seed Protein-type BURP-domain subfamily, catalyzes intramolecular macrocyclizations of its core peptide during the sequential biosynthesis of monocyclic lyciumin I via glycine-tryptophan crosslinking and bicyclic legumenin via glutamine-tyrosine crosslinking. X-ray crystallography of AhyBURP reveals the BURP-domain fold with two type II copper centers derived from a conserved stapled-disulfide and His motif. We show the macrocyclization of lyciumin-C(sp3)-N-bond formation followed by legumenin-C(sp3)-O-bond formation requires dioxygen and radical involvement based on enzyme assays in anoxic conditions and isotopic labeling. Our study expands enzymatic intramolecular modifications beyond catalytic moiety and chromophore biogenesis to RiPP biosynthesis.

Transcriptome Profiling of a Soybean Mutant with Salt Tolerance Induced by Gamma-ray Irradiation

Salt stress is a significant abiotic stress that reduces crop yield and quality globally. In this study, we utilized RNA sequencing (RNA-Seq) to identify differentially expressed genes (DEGs) in response to salt stress induced by gamma-ray irradiation in a salt-tolerant soybean mutant. The total RNA library samples were obtained from the salt-sensitive soybean cultivar Kwangan and the salt-tolerant mutant KA-1285. Samples were taken at three time points (0, 24, and 72 h) from two tissues (leaves and roots) under 200 mM NaCl. A total of 967,719,358 clean reads were generated using the Illumina NovaSeq 6000 platform, and 94.48% of these reads were mapped to 56,044 gene models of the soybean reference genome (Glycine_max_Wm82.a2.v1). The DEGs with expression values were compared at each time point within each tissue between the two soybeans. As a result, 296 DEGs were identified in the leaves, while 170 DEGs were identified in the roots. In the case of the leaves, eight DEGs were related to the phenylpropanoid biosynthesis pathway; however, in the roots, Glyma.03G171700 within GmSalt3, a major QTL associated with salt tolerance in soybean plants, was differentially expressed. Overall, these differences may explain the mechanisms through which mutants exhibit enhanced tolerance to salt stress, and they may provide a basic understanding of salt tolerance in soybean plants.

Genome-wide identification and expression analysis of the BURP domain-containing genes in Malus domestica

The conserved BURP-containing proteins are specific to plants and play a crucial role in plant growth, development, and response to abiotic stresses. However, less is known about the systematic characterization of BURP-containing proteins in apple. This study aimed to identify and analyze all BURP-containing genes in the apple genome, as well as to examine their expression patterns through various bioinformatics methods. Eighteen members of BURP-containing genes were identified in apple, six members lacked signal peptides, and the secondary structure was mainly a Random coil of BURP-containing genes. Gene structure and Motif analysis showed that proteins have similar structures and are conserved at the C-terminal. Cis-acting element analysis revealed that the proteins contain phytohormone and stress response elements, and chromosomal localization revealed that the family is unevenly distributed across eight chromosomes, with duplication of fragments leading to the expansion of family proteins. Tissue expression showed that MdPG3 and MdPG4 were expressed in different tissues and different varieties, MdRD2 and MdRD7 were highly expressed in 'M74' fruits and MdRD7 in 'M49' leaves, while MdUSP1 was highly expressed in 'GD' roots. The quantitative real-time PCR analysis showed that the expressions of six and seven genes were significantly up-regulated under NaCl and PEG treatments, respectively, whereas MdRD7 was significantly up-regulated under NaCl and PEG treatment over time. This study offers a comprehensive identification and expression analysis of BURP-containing proteins in apple. The findings provide a theoretical foundation for further exploration of the functions of this protein family.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01393-7.

Cre-mediated autoexcision of selectable marker genes in soybean, cotton, canola and maize transgenic plants

Efficient selectable marker gene autoexcision in transgenic plants of soybean, cotton, canola, and maize is achieved by effective Cre recombinase expression. Selectable marker genes are often required for efficient generation of transgenic plants in plant transformation but are not desired once the transgenic events are obtained. We have developed Cre/loxP autoexcision systems to remove selectable marker genes in soybean, cotton, canola and maize. We tested a set of vectors with diverse promoters and identified promising promoters to drive cre expression for each of the four crops. We evaluated both the efficiency of generating primary transgenic events with low transgene copy numbers, and the frequency of marker-free progeny in the next generation. The best performing vectors gave no obvious decrease in the transformation frequency in each crop and generated homozygous marker-free progeny in the next generation. We found that effective expression of Cre recombinase for marker gene autoexcision can be species dependent. Among the vectors tested, the best autoexcision frequency (41%) in soybean transformation came from using the soybean RSP1 promoter for cre expression. The cre gene expressed by soybean RSP1 promoter with an Arabidopsis AtpE intron delivered the best autoexcision frequency (69%) in cotton transformation. The cre gene expressed by the embryo-specific eUSP88 promoter from Vicia faba conferred the best marker excision frequency (32%) in canola transformation. Finally, the cre gene expressed by the rice CDC45-1 promoter resulted in 44% autoexcision in maize transformation. The Cre/loxP recombinase system enables the generation of selectable marker-free transgenic plants for commercial product development in four agriculturally important crops and provides further improvement opportunities for more specific and better marker excision efficiency.

Identification and expression analysis of BURP domain-containing genes in jujube and their involvement in low temperature and drought response

Background

Plant-specific BURP domain-containing genes are involved in plant development and stress responses. However, the role of BURP family in jujube (Ziziphus jujuba Mill.) has not been investigated.

Results

In this study, 17 BURP genes belonging to four subfamilies were identified in jujube based on homology analysis, gene structures, and conserved motif confirmation. Gene duplication analysis indicated both tandem duplication and segmental duplication had contributed to ZjBURP expansion. The ZjBURPs were extensively expressed in flowers, young fruits, and jujube leaves. Transcriptomic data and qRT-PCR analysis further revealed that ZjBURPs also significantly influence fruit development, and most genes could be induced by low temperature, salinity, and drought stresses. Notably, several BURP genes significantly altered expression in response to low temperature (ZjPG1) and drought stresses (ZjBNM7, ZjBNM8, and ZjBNM9).

Conclusions

These results provided insights into the possible roles of ZjBURPs in jujube development and stress response. These findings would help selecting candidate ZjBURP genes for cold- and drought-tolerant jujube breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08907-9.

Whole-Genome Identification of APX and CAT Gene Families in Cultivated and Wild Soybeans and Their Regulatory Function in Plant Development and Stress Response

Plants coevolved with their antioxidant defense systems, which detoxify and adjust levels of reactive oxygen species (ROS) under multiple plant stresses. We performed whole-genome identification of ascorbate peroxidase (APX) and catalase (CAT) families in cultivated and wild soybeans. In cultivated and wild soybean genomes, we identified 11 and 10 APX genes, respectively, whereas the numbers of identified CAT genes were four in each species. Comparative phylogenetic analysis revealed more homology among cultivated and wild soybeans relative to other legumes. Exon/intron structure, motif and synteny blocks are conserved in cultivated and wild species. According to the Ka/Ks value, purifying selection is a major force for evolution of these gene families in wild soybean; however, the APX gene family was evolved by both positive and purifying selection in cultivated soybean. Segmental duplication was a major factor involved in the expansion of APX and CAT genes. Expression patterns revealed that APX and CAT genes are differentially expressed across fourteen different soybean tissues under water deficit (WD), heat stress (HS) and combined drought plus heat stress (WD + HS). Altogether, the current study provides broad insights into these gene families in soybeans. Our results indicate that APX and CAT gene families modulate multiple stress response in soybeans.

Arabidopsis ETHYLENE INSENSITIVE 3 directly regulates the expression of PG1β-like family genes in response to aluminum stress

The genes in the subfamily PG1β (beta subunit of poly-galacturonase isoenzyme 1) have a clear effect on the biosynthesis pathway of pectin, a main component of the cell wall. However, the detailed functions of the PG1β-like gene members in Arabidopsis (AtPG1-3) have not yet been determined. In this study, we investigated their functional roles in response to aluminum (Al) stress. Our results indicate that the PG1β-like gene members are indeed involved in the Al-stress response and they can modulate its accumulation in roots to achieve optimum root elongation and hence better seedling growth. We found that transcription factor EIN3 (ETHYLENE INSENSITIVE 3) alters pectin metabolism and the EIN3 gene responds to Al stress to affect the pectin content in the root cell walls, leading to exacerbation of the inhibition of root growth, as reflected by the phenotypes of overexpressing lines. We determined that EIN3 can directly bind to the promoter regions of PG1-3, which act downstream of EIN3. Thus, our results show that EIN3 responds to Al stress in Arabidopsis directly through regulating the expression of PG1-3. Hence, EIN3 mediates their functions by acting as a biomarker in their molecular biosynthesis pathways, and consequently orchestrates their biological network in response to Al stress.

In silico analysis and expression profiling of Expansin A4, BURP domain protein RD22-like and E6-like genes associated with fiber quality in cotton

BACKGROUND

To supply high-quality cotton fibre for the textile industry, the development of long, strong and fine fibre cotton varieties is imperative. An interlinked approach was used to comprehend the role of fibre genes by analyzing interspecific progenies of cotton species. Wild Gossypium species and races are rich source of genetic polymorphism due to environmental dispersal and continuous natural selection. These genetic resources hold mass of outclass genes that can be used in cotton improvement breeding programs to exploit possible traits such as fibre quality, abiotic stress tolerance, and disease and insect resistance. Therefore, use of new molecular techniques such as genomics, transcriptomics and bioinformatics is very important to utilize the genetic potential of wild species in cotton improvement programs.

METHODS

Interspecific lines and Gossypium species used in the study were grown at Central Cotton Research Institute (CCRI), Multan. After retrieving DNA sequence of the genes from NCBI, the primers for gene expression and full-length gene sequence were designed. Expression profiling of Expansin A4, BURP Domain protein RD22-like and E6-like fibre genes was performed through Real Time PCR. BLAST and DNA sequence alignment was conducted for sequence comparison of interspecific lines and Gossypium species. Different in silico analysis were used for characterization of fibre genes and identification of cis acting promoter elements in promoter region.

RESULTS

Variable expression of genes related to fibre development was observed at different stages. BLAST and DNA sequence alignment demonstrated resemblance of interspecific lines with G. hirsutum. In silico analysis on the sequence data also confirmed the role of Expansin A4, BURP Domain protein RD22-like and E6-like fibre genes in fibre development. Genetic engineering is also recommended by transferring E6-like, Expansin A4 and BURP Domain RD22-like genes in local cotton cultivars. Similarly, several stress tolerant and light responsive cis acting elements were identified through promotor analysis, which may contribute for fibre development in the breeding programs.

CONCLUSION

Expansin A4, BURP Domain RD22-like and E6-like have positive role in fibre development with variable expression at fiber length and strength associated stages.

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.

Genome-wide analysis of BURP genes and identification of a BURP-V gene RcBURP4 in Rosa chinensis

Nine RcBURPs have been identified in Rosa chinensis, and overexpression of RcBURP4 increased ABA, NaCl sensitivity, and drought tolerance in transgenic Arabidopsis. BURP proteins are unique to plants and may contribute greatly to growth, development, and stress responses of plants. Despite the vital role of BURP proteins, little is known about these proteins in rose (Rosa spp.). In the present study, nine genes belonging to the BURP family in R. chinensis were identified using multiple bioinformatic approaches against the rose genome database. The nine RcBURPs, with diverse structures, were located on all chromosomes of the rose genome, except for Chr2 and Chr3. Phylogenic analysis revealed that these RcBURPs can be classified into eight subfamilies, including BNM2-like, PG1β-like, USP-like, RD22-like, BURP-V, BURP-VI, BURP-VII, and BURP-VIII. Conserved motif and exon-intron analyses indicated a conserved pattern within the same subfamily. The presumed cis-regulatory elements (CREs) within the promoter region of each RcBURP were analyzed and the results showed that all RcBURPs contained different types of CREs, including abiotic stress-, light response-, phytohormones response-, and plant growth and development-related CREs. Transcriptomic analysis revealed that a BURP-V member, RcBURP4, was induced in rose leaves and roots under mild and severe drought treatments. We then overexpressed RcBURP4 in Arabidopsis and examined its role under abscisic acid (ABA), NaCl, polyethylene glycol (PEG), and drought treatments. Nine stress-responsive genes expression were changed in RcBURP4-overexpressing leaves and roots. Furthermore, RcBURP4-silenced rose plants exhibited decreased tolerance to dehydration. The results obtained from this study provide the first comprehensive overview of RcBURPs and highlight the importance of RcBURP4 in rose plant.

Discovery and biosynthesis of cyclic plant peptides via autocatalytic cyclases

Many bioactive plant cyclic peptides form side-chain-derived macrocycles. Lyciumins, cyclic plant peptides with tryptophan macrocyclizations, are ribosomal peptides (RiPPs) originating from repetitive core peptide motifs in precursor peptides with plant-specific BURP (BNM2, USP, RD22 and PG1beta) domains, but the biosynthetic mechanism for their formation has remained unknown. Here, we characterize precursor-peptide BURP domains as copper-dependent autocatalytic peptide cyclases and use a combination of tandem mass spectrometry-based metabolomics and plant genomics to systematically discover five BURP-domain-derived plant RiPP classes, with mono- and bicyclic structures formed via tryptophans and tyrosines, from botanical collections. As BURP-domain cyclases are scaffold-generating enzymes in plant specialized metabolism that are physically connected to their substrates in the same polypeptide, we introduce a bioinformatic method to mine plant genomes for precursor-peptide-encoding genes by detection of repetitive substrate domains and known core peptide features. Our study sets the stage for chemical, biosynthetic and biological exploration of plant RiPP natural products from BURP-domain cyclases.

Genome-wide identification of the BURP domain-containing genes in Phaseolus vulgaris

Plant-specific BURP domain-containing proteins have an essential role in the plant's development and stress responses. Although BURP domain-containing proteins have been identified in several plant species, genome-wide analysis of the BURP gene family has not been investigated in the common bean. In the present study, we identified 11 BURP family members in the common bean (Phaseolus vulgaris) genome with a comprehensive in silico analysis. Pairwise alignment and phylogenetic analyses grouped PvBURP members into four subfamilies [RD-22 like (3), PG1β-like (4), BNM2-like (3), and USP-like (1)] according to their amino acid motifs, protein domains and intron-exon structure. The physical and biochemical characteristics of amino acids, motif and intron-exon structure, and cis-regulatory elements of BURPs members were determined. Promoter regions of BURP members included stress, light, and hormone response-related cis-elements. Therefore, expression profiles of PvBURP genes were identified with in silico tools and qRT-PCR analyses under stress (salt and drought) and hormone treatment (ABA, IAA) in the current study. While significant activity changes were not observed in BURP genes in RNA-seq data sets related to salt stress, it was determined that some BURP genes were expressed differently in those with drought stress. We identified 12 different miRNA, including miRNA395, miRNA156, miRNA169, miRNA171, miRNA319, and miRNA390, targeting the nine PvBURP genes using two different in silico tools based on perfect or near-perfect complementarity to their targets. Here we present the first study to identify and characterize the BURP genes in common bean using whole-genome analysis, and the findings may serve as a reference for future functional research in common bean.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01052-9.

Local adaptation contributes to gene expression divergence in maize

Gene expression links genotypes to phenotypes, so identifying genes whose expression is shaped by selection will be important for understanding the traits and processes underlying local adaptation. However, detecting local adaptation for gene expression will require distinguishing between divergence due to selection and divergence due to genetic drift. Here, we adapt a QST−FST framework to detect local adaptation for transcriptome-wide gene expression levels in a population of diverse maize genotypes. We compare the number and types of selected genes across a wide range of maize populations and tissues, as well as selection on cold-response genes, drought-response genes, and coexpression clusters. We identify a number of genes whose expression levels are consistent with local adaptation and show that genes involved in stress response show enrichment for selection. Due to its history of intense selective breeding and domestication, maize evolution has long been of interest to researchers, and our study provides insight into the genes and processes important for in local adaptation of maize.

Genome-wide identification and expression analysis of the BURP domain-containing genes in Gossypium hirsutum

Background

Many BURP domain-containing proteins, which are unique to plants, have been identified. They performed diverse functions in plant development and the stress response. To date, only a few BURP domain-containing genes have been studied, and no comprehensive analysis of the gene family in cotton has been reported.

Results

In this study, 18, 17 and 30 putative BURP genes were identified in G. raimondii (D₅), G. arboreum (A₂) and G. hirsutum (AD₁), respectively. These BURP genes were phylogenetically classified into eight subfamilies, which were confirmed by analyses of gene structures, motifs and protein domains. The uneven distribution of BURPs in chromosomes and gene duplication analysis indicated that segmental duplication might be the main driving force of the GhBURP family expansion. Promoter regions of all GhBURPs contained at least one putative stress-related cis-elements. Analysis of transcriptomic data and qRT-PCR showed that GhBURPs showed different expression patterns in different organs, and all of them, especially the members of the RD22-like subfamily, could be induced by different stresses, such as abscisic acid (ABA) and salicylic acid (SA), which indicated that the GhBURPs may performed important functions in cotton’s responses to various abiotic stresses.

Conclusions

Our study comprehensively analyzed BURP genes in G. hirsutum, providing insight into the functions of GhBURPs in cotton development and adaptation to stresses.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5948-y) contains supplementary material, which is available to authorized users.

RNA sequencing analysis of salt tolerance in soybean (Glycine max)

Salt stress causes foliar chlorosis and scorch, plant stunting, and eventually yield reduction in soybean. There are differential responses, namely tolerance (excluder) and intolerance (includer), among soybean germplasm. However, the genetic and physiological mechanisms for salt tolerance is complex and not clear yet. Based on the results from the screening of the RA-452 x Osage mapping population, two F4:6 lines with extreme responses, most tolerant and most sensitive, were selected for a time-course gene expression study in which the 250 mM NaCl treatment was initially imposed at the V1 stage and continued for 24 h (hrs). Total RNA was isolated from the leaves harvested at 0, 6, 12, 24 h after the initiation of salt treatment, respectively. The RNA-Seq analysis was conducted to compare the salt tolerant genotype with salt sensitive genotype at each time point using RNA-Seq pipeline method. A total of 2374, 998, 1746, and 630 differentially expressed genes (DEGs) between salt-tolerant line and salt-sensitive line, were found at 0, 6, 12, and 24 h, respectively. The expression patterns of 154 common DEGs among all the time points were investigated, of which, six common DEGs were upregulated and seven common DEGs were downregulated in salt-tolerant line. Moreover, 13 common DEGs were dramatically expressed at all the time points. Based on Log2 (fold change) of expression level of salt-tolerant line to salt-sensitive line and gene annotation, Glyma.02G228100, Glyma.03G226000, Glyma.03G031000, Glyma.03G031400, Glyma.04G180300, Glyma.04G180400, Glyma.05 g204600, Glyma.08G189600, Glyma.13G042200, and Glyma.17G173200, were considered to be the key potential genes involving in the salt-tolerance mechanism in the soybean salt-tolerant line.

A quantitative proteomic analysis of the molecular mechanism underlying fertility conversion in thermo-sensitive genetic male sterility line AnnongS-1

BACKGROUND

Thermo-sensitive genetic male sterile (TGMS) lines have been widely used in two-line hybrid rice breeding. The two-line hybrids have increased rice yields substantially. However, the effect of environmental temperatures on the fertility conversion is still not fully clear. In this study, we performed a tandem mass tag (TMT)-based proteomic analysis on the anthers of the TGMS line AnnongS-1 grown under permissive (low) temperature (21 °C) and restrictive (high) temperature (> 26 °C) conditions in an attempt to explore the effect of temperature on the fertility of the male sterile line.

RESULTS

After the AnnongS-1 plants were induced under either permissive or restrictive conditions, morphological observations and I₂-KI staining confirmed that the pollen grains formed under high temperature conditions were abortive while those formed under low temperature developed normally. In comparison to the plants grown under permissive conditions, the restrictive high-temperature conditions led to the differential accumulation of 89 proteins in the anthers, of which 46 were increased in abundance and 43 were decreased in abundance. Most of the subcellular compartments of the anther cells had one or more proteins that had been differentially accumulated, with the cytoplasm and chloroplast having the greatest accumulations. More than 40% of the differentially abundant proteins (DAPs) were enzymes involved in photosynthesis, energy metabolism, biosynthesis and catabolism of cellular components, metabolic regulation, defense and stress, etc. The DAPs related to protein metabolism accounted for the largest proportion (21.35%), followed by those related to defense and stress (12.36%), metabolic regulation (10.11%) and carbohydrate metabolism (8.99%), indicating that such biological processes in anther cells were more susceptible to high temperature stress.

CONCLUSIONS

The restrictive temperature induction caused fertility-sterility conversion in the TGMS line AnnongS-1 mainly by adversely affecting the metabolism of protein, carbohydrate and energy, and decreasing the abundances of important proteins closely related to defense and stress, thereby impeding the growth and development of the pollen and weakening the overall defense and ability to endure stress of AnnongS-1. These data are helpful for deepening our understanding of the molecular mechanism underlying fertility conversion in TGMS lines.

Over-expression of GmKR3, a TIR-NBS-LRR type R gene, confers resistance to multiple viruses in soybean

KEY MESSAGE

That overexpression of GmKR3 enhances innate virus resistance by stimulating. Soybean mosaic virus (SMV) is found in many soybean production areas, and SMV infection is one of the prevalent viral diseases that can cause significant yield losses in soybean. In plants, resistance (R) genes are involved in pathogen reorganization and innate immune response activation. Most R proteins have nucleotide-binding site and leucine-rich repeat (NBS-LRR) domains, and some of the NBS-LRR type R proteins in dicots have Toll/Interleukin-1 Receptor (TIR) motifs. We report here the analysis of the over-expression of GmKR3, a soybean TIR-NBS-LRR type R gene on virus resistance in soybean. When over-expressed in soybean, GmKR3 enhanced the plant's resistance to several strains of SMV, the closely related potyviruses bean common mosaic virus (BCMV) and watermelon mosaic virus (WMV), and the secovirus bean pod mottle virus (BPMV). Importantly, over-expression of GmKR3 did not affect plant growth and development, including yield and qualities of the seeds. HPLC analysis showed that abscisic acid (ABA) content increased in the 35S:GmKR3 transgenic plants, and in both wild-type and 35S:GmKR3 transgenic plants in response to virus inoculation. Consistent with this observation, we found that the expression of two ABA catabolism genes was down-regulated in 35S:GmKR3 transgenic plants. We also found that the expression of Gm04.3, an ABA responsive gene encoding BURP domain-containing protein, was up-regulated in 35S:GmKR3 transgenic plants. Taken together, our results suggest that overexpression of GmKR3 enhanced virus resistance in soybean, which was achieved at least in part via ABA signaling.

Assessment of shared alleles in drought-associated candidate genes among southern California white oak species (Quercus sect. Quercus)

BACKGROUND

Hybridization and introgression are common phenomena among oak species. These processes can be beneficial by introducing favorable genetic variants across species (adaptive introgression). Given that drought is an important stress, impacting physiological and morphological variation and limiting distributions, our goal was to identify drought-related genes that might exhibit patterns of introgression influenced by natural selection. Using RNAseq, we sequenced whole transcriptomes of 24 individuals from three oaks in southern California: (Quercus engelmannii, Quercus berberidifolia, Quercus cornelius-mulleri) and identified genetic variants to estimate admixture rates of all variants and those in drought genes.

RESULTS

We found 398,042 variants across all loci and 4352 variants in 139 drought candidate genes. STRUCTURE analysis of all variants revealed the majority of our samples were assignable to a single species, but with several highly admixed individuals. When using drought-associated variants, the same individuals exhibited less admixture and their allele frequencies were more polarized between Engelmann and scrub oaks than when using the total gene set. These findings are consistent with the hypothesis that selection may act differently on functional genes, such as drought-associated genes, and point to candidate genes that are suggestive of divergent selection among species maintaining adaptive differences. For example, the drought genes that showed the strongest bias against engelmannii-fixed oak variants in scrub oaks were related to sugar transporter, coumarate-coA ligases, glutathione S-conjugation, and stress response.

CONCLUSION

This pilot study illustrates that whole transcriptomes of individuals will provide useful data for identifying functional genes that contribute to adaptive divergence among hybridizing species.

Transcriptomic analysis of boron hyperaccumulation mechanisms in Puccinellia distans

Puccinellia distans, common alkali grass, is found throughout the world and can survive in soils with boron concentrations that are lethal for other plant species. Indeed, P. distans accumulates very high levels of this element. Despite these interesting features, very little research has been performed to elucidate the boron tolerance mechanism in this species. In this study, P. distans samples were treated for three weeks with normal (0.5 mg L^-1) and elevated (500 mg L^-1) boron levels in hydroponic solution. Expressed sequence tags (ESTs) derived from shoot tissue were analyzed by RNA sequencing to identify genes up and down-regulated under boron stress. In this way, 3312 differentially expressed transcripts were detected, 67.7% of which were up-regulated and 32.3% of which were down-regulated in boron-treated plants. To partially confirm the RNA sequencing results, 32 randomly selected transcripts were analyzed for their expression levels in boron-treated plants. The results agreed with the expected direction of change (up or down-regulation). A total of 1652 transcripts had homologs in A. thaliana and/or O. sativa and mapped to 1107 different proteins. Functional annotation of these proteins indicated that the boron tolerance and hyperaccumulation mechanisms of P. distans involve many transcriptomic changes including: alterations in the malate pathway, changes in cell wall components that may allow sequestration of excess boron without toxic effects, and increased expression of at least one putative boron transporter and two putative aquaporins. Elucidation of the boron accumulation mechanism is important in developing approaches for bioremediation of boron contaminated soils.

An endoplasmic reticulum-localized Coffea arabica BURP domain-containing protein affects the response of transgenic Arabidopsis plants to diverse abiotic stresses

The Coffea arabica BURP domain-containing gene plays an important role in the response of transgenic Arabidopsis plants to abiotic stresses via regulating the level of diverse proteins. Although the functions of plant-specific BURP domain-containing proteins (BDP) have been determined for a few plants, their roles in the growth, development, and stress responses of most plant species, including coffee plant (Coffea arabica), are largely unknown. In this study, the function of a C. arabica BDP, designated CaBDP1, was investigated in transgenic Arabidopsis plants. The expression of CaBDP1 was highly modulated in coffee plants subjected to drought, cold, salt, or ABA. Confocal analysis of CaBDP1-GFP fusion proteins revealed that CaBDP1 is localized in the endoplasmic reticulum. The ectopic expression of CaBDP1 in Arabidopsis resulted in delayed germination of the transgenic plants under abiotic stress and in the presence of ABA. Cotyledon greening and seedling growth of the transgenic plants were inhibited in the presence of ABA due to the upregulation of ABA signaling-related genes like ABI3, ABI4, and ABI5. Proteome analysis revealed that the levels of several proteins are modulated in CaBDP1-expressing transgenic plants. The results of this study underscore the importance of BURP domain proteins in plant responses to diverse abiotic stresses.

Drought and exogenous abscisic acid alter hydrogen peroxide accumulation and differentially regulate the expression of two maize RD22-like genes

Increased biosynthesis of abscisic acid (ABA) occurs in plants in response to water deficit, which is mediated by changes in the levels of reactive oxygen species such as H₂O₂. Water deficit and ABA induce expression of some RD22-like proteins. This study aimed to evaluate the effect of water deficit and exogenous ABA (50 µM ABA applied every 24 hours for a total of 72 hours) on H₂O₂ content in Zea mays (maize) and to characterise genes encoding two putative maize RD22-like proteins (designated ZmRD22A and ZmRD22B). The expression profiles of the two putative maize RD22-like genes in response to water deficit and treatment with ABA were examined in leaves. In silico analyses showed that the maize RD22-like proteins share domain organisation with previously characterized RD22-like proteins. Both water deficit and exogenous ABA resulted in increased H₂O₂ content in leaves but the increase was more pronounced in response to water deficit than to exogenous ABA. Lignin content was not affected by exogenous ABA, whereas it was decreased by water deficit. Expression of both RD22-like genes was up-regulated by drought but the ZmRD22A gene was not influenced by exogenous ABA, whereas ZmRD22B was highly responsive to exogenous ABA.

Didehydrophenylalanine, an abundant modification in the beta subunit of plant polygalacturonases

The structure and the activity of proteins are often regulated by transient or stable post- translational modifications (PTM). Different from well-known, abundant modifications such as phosphorylation and glycosylation some modifications are limited to one or a few proteins across a broad range of related species. Although few examples of the latter type are known, the evolutionary conservation of these modifications and the enzymes responsible for their synthesis suggest an important physiological role. Here, the first observation of a new, fold-directing PTM is described. During the analysis of alfalfa cell wall proteins a -2Da mass shift was observed on phenylalanine residues in the repeated tetrapeptide FxxY of the beta-subunit of polygalacturonase. This modular protein is known to be involved in developmental and stress-responsive processes. The presence of this modification was confirmed using in-house and external datasets acquired by different commonly used techniques in proteome studies. Based on these analyses it was found that all identified phenylalanine residues in the sequence FxxY of this protein were modified to α,β-didehydro-Phe (ΔPhe). Besides showing the reproducible identification of ΔPhe in different species arguments that substantiate the fold-determining role of ΔPhe are given.

Genome-wide analysis and expression profiling of glyoxalase gene families in soybean (Glycine max) indicate their development and abiotic stress specific response

BACKGROUND

Glyoxalase pathway consists of two enzymes, glyoxalase I (GLYI) and glyoxalase II (GLYII) which detoxifies a highly cytotoxic metabolite methylglyoxal (MG) to its non-toxic form. MG may form advanced glycation end products with various cellular macro-molecules such as proteins, DNA and RNA; that ultimately lead to their inactivation. Role of glyoxalase enzymes has been extensively investigated in various plant species which showed their crucial role in salinity, drought and heavy metal stress tolerance. Previously genome-wide analysis of glyoxalase genes has been conducted in model plants Arabidopsis and rice, but no such study was performed in any legume species.

RESULTS

In the present study, a comprehensive genome database analysis of soybean was performed and identified a total of putative 41 GLYI and 23 GLYII proteins encoded by 24 and 12 genes, respectively. Detailed analysis of these identified members was conducted including their nomenclature and classification, chromosomal distribution and duplication, exon-intron organization, and protein domain(s) and motifs identification. Expression profiling of these genes has been performed in different tissues and developmental stages as well as under salinity and drought stresses using publicly available RNAseq and microarray data. The study revealed that GmGLYI-7 and GmGLYII-8 have been expressed intensively in all the developmental stages and tissues; while GmGLYI-6, GmGLYI-9, GmGLYI-20, GmGLYII-5 and GmGLYII-10 were highly abiotic stress responsive members.

CONCLUSIONS

The present study identifies the largest family of glyoxalase proteins to date with 41 GmGLYI and 23 GmGLYII members in soybean. Detailed analysis of GmGLYI and GmGLYII genes strongly indicates the genome-wide segmental and tandem duplication of the glyoxalase members. Moreover, this study provides a strong basis about the biological role and function of GmGLYI and GmGLYII members in soybean growth, development and stress physiology.

Identification and Expression Analysis of BURP Domain-Containing Genes in Medicago truncatula

BURP domain-containing proteins belong to a newly identified protein class that is unique to plants and plays an important role in plant development and metabolism. Although systematic characterization of BURP domain-containing proteins have been carried out in many species, such as rice, poplar and maize, little is known about BURP domain-containing proteins in Medicago. In this study, multiple bioinformatics approaches were employed to identify all the members of BURP family genes in Medicago. A complete set of 39 BURP family genes were identified. These genes have diverse structures and were distributed on chromosome 1-8 except 7. According to phylogenetic analysis, these BURP family genes could be classified into eight classes. Motif and exon-intron organization, stress-related cis-elements in promoter regions and microarray analysis of MtBURPs were also performed. Furthermore, transcript level analysis of MtBURP genes in response to drought stress revealed that all of the 39 BURP genes were regulated by drought stress. The results of this study reveal a comprehensive overview of the Medicago BURP gene family and provide the first step toward the selection of MtBURP genes for cloning and functional analysis of the BURP gene family in Medicago truncatula.

The divergence and positive selection of the plant-specific BURP-containing protein family

BURP domain-containing proteins belong to a plant-specific protein family and have diverse roles in plant development and stress responses. However, our understanding about the genetic divergence patterns and evolutionary rates of these proteins remain inadequate. In this study, 15 plant genomes were explored to elucidate the genetic origins, divergence, and functions of these proteins. One hundred and twenty-five BURP protein-encoding genes were identified from four main plant lineages, including 13 higher plant species. The absence of BURP family genes in unicellular and multicellular algae suggests that this family (1) appeared when plants shifted from relatively stable aquatic environments to land, where conditions are more variable and stressful, and (2) is critical in the adaptation of plants to adverse environments. Promoter analysis revealed that several responsive elements to plant hormones and external environment stresses are concentrated in the promoter region of BURP protein-encoding genes. This finding confirms that these genes influence plant stress responses. Several segmentally and tandem-duplicated gene pairs were identified from eight plant species. Thus, in general, BURP domain-containing genes have been subject to strong positive selection, even though these genes have conformed to different expansion models in different species. Our study also detected certain critical amino acid sites that may have contributed to functional divergence among groups or subgroups. Unexpectedly, all of the critical amino acid residues of functional divergence and positive selection were exclusively located in the C-terminal region of the BURP domain. In conclusion, our results contribute novel insights into the genetic divergence patterns and evolutionary rates of BURP proteins.

Daytime soybean transcriptome fluctuations during water deficit stress

BACKGROUND

Since drought can seriously affect plant growth and development and little is known about how the oscillations of gene expression during the drought stress-acclimation response in soybean is affected, we applied Illumina technology to sequence 36 cDNA libraries synthesized from control and drought-stressed soybean plants to verify the dynamic changes in gene expression during a 24-h time course. Cycling variables were measured from the expression data to determine the putative circadian rhythm regulation of gene expression.

RESULTS

We identified 4866 genes differentially expressed in soybean plants in response to water deficit. Of these genes, 3715 were differentially expressed during the light period, from which approximately 9.55% were observed in both light and darkness. We found 887 genes that were either up- or down-regulated in different periods of the day. Of 54,175 predicted soybean genes, 35.52% exhibited expression oscillations in a 24 h period. This number increased to 39.23% when plants were submitted to water deficit. Major differences in gene expression were observed in the control plants from late day (ZT16) until predawn (ZT20) periods, indicating that gene expression oscillates during the course of 24 h in normal development. Under water deficit, dissimilarity increased in all time-periods, indicating that the applied stress influenced gene expression. Such differences in plants under stress were primarily observed in ZT0 (early morning) to ZT8 (late day) and also from ZT4 to ZT12. Stress-related pathways were triggered in response to water deficit primarily during midday, when more genes were up-regulated compared to early morning. Additionally, genes known to be involved in secondary metabolism and hormone signaling were also expressed in the dark period.

CONCLUSIONS

Gene expression networks can be dynamically shaped to acclimate plant metabolism under environmental stressful conditions. We have identified putative cycling genes that are expressed in soybean leaves under normal developmental conditions and genes whose expression oscillates under conditions of water deficit. These results suggest that time of day, as well as light and temperature oscillations that occur considerably affect the regulation of water deficit stress response in soybean plants.

Characterization of the cork oak transcriptome dynamics during acorn development

BACKGROUND

Cork oak (Quercus suber L.) has a natural distribution across western Mediterranean regions and is a keystone forest tree species in these ecosystems. The fruiting phase is especially critical for its regeneration but the molecular mechanisms underlying the biochemical and physiological changes during cork oak acorn development are poorly understood. In this study, the transcriptome of the cork oak acorn, including the seed, was characterized in five stages of development, from early development to acorn maturation, to identify the dominant processes in each stage and reveal transcripts with important functions in gene expression regulation and response to water.

RESULTS

A total of 80,357 expressed sequence tags (ESTs) were de novo assembled from RNA-Seq libraries representative of the several acorn developmental stages. Approximately 7.6 % of the total number of transcripts present in Q. suber transcriptome was identified as acorn specific. The analysis of expression profiles during development returned 2,285 differentially expressed (DE) transcripts, which were clustered into six groups. The stage of development corresponding to the mature acorn exhibited an expression profile markedly different from other stages. Approximately 22 % of the DE transcripts putatively code for transcription factors (TF) or transcriptional regulators, and were found almost equally distributed among the several expression profile clusters, highlighting their major roles in controlling the whole developmental process. On the other hand, carbohydrate metabolism, the biological pathway most represented during acorn development, was especially prevalent in mid to late stages as evidenced by enrichment analysis. We further show that genes related to response to water, water deprivation and transport were mostly represented during the early (S2) and the last stage (S8) of acorn development, when tolerance to water desiccation is possibly critical for acorn viability.

CONCLUSIONS

To our knowledge this work represents the first report of acorn development transcriptomics in oaks. The obtained results provide novel insights into the developmental biology of cork oak acorns, highlighting transcripts putatively involved in the regulation of the gene expression program and in specific processes likely essential for adaptation. It is expected that this knowledge can be transferred to other oak species of great ecological value.

AtPGL3 is an Arabidopsis BURP domain protein that is localized to the cell wall and promotes cell enlargement

The BURP domain is a plant-specific domain that has been identified in secretory proteins, and some of these are involved in cell wall modification. The tomato polygalacturonase I complex involved in pectin degradation in ripening fruits has a non-catalytic subunit that has a BURP domain. This protein is called polygalacturonase 1 beta (PG1β) and the Arabidopsis genome encodes three proteins that exhibit strong amino acid similarities with PG1β? We generated Arabidopsis lines in which expression levels of AtPGLs are altered in order to investigate the biological roles of the Arabidopsis PG1β-like proteins (AtPGLs). Among the three AtPGLs (AtPGL1-3), AtPGL3 exhibited the highest transcriptional activity throughout all developmental stages. AtPGL triple mutant plants have smaller rosette leaves than those of wild type plants because the leaf cells are smaller in the mutant plants. Interestingly, when we overexpressed AtPGL3 using a 35S promoter, leaf cells in transgenic plants grew larger than those of the wild type. A C-terminal GFP fusion protein of AtPGL3 complemented phenotypes of the triple mutant plants and it localized to the cell wall. A truncated AtPGL3-GFP fusion protein lacking the BURP domain failed to rescue the mutant phenotypes even though the GFP protein was targeted to the cell wall, indicating that the BURP domain is required for the protein's effect on cell expansion. Quantitative RT-PCR and immunoblot analyses indicated that the α-expansin 6 gene is up-regulated in the overexpressor plants. Taken together, these results indicate that AtPGL3 is an apoplastic BURP domain protein playing a role in cell expansion.

Inspection of the grapevine BURP superfamily highlights an expansion of RD22 genes with distinctive expression features in berry development and ABA-mediated stress responses

The RESPONSIVE TO DEHYDRATION 22 (RD22) gene is a molecular link between abscisic acid (ABA) signalling and abiotic stress responses. Its expression has been used as a reliable ABA early response marker. In Arabidopsis, the single copy RD22 gene possesses a BURP domain also located at the C-terminus of USP embryonic proteins and the beta subunit of polygalacturonases. In grapevine, a RD22 gene has been identified but putative paralogs are also found in the grape genome, possibly forming a large RD22 family in this species. In this work, we searched for annotations containing BURP domains in the Vitis vinifera genome. Nineteen proteins were defined by a comparative analysis between the two genome predictions and RNA-Seq data. These sequences were compared to other plant BURPs identified in previous genome surveys allowing us to reconceive group classifications based on phylogenetic relationships and protein motif occurrence. We observed a lineage-specific evolution of the RD22 family, with the biggest expansion in grapevine and poplar. In contrast, rice, sorghum and maize presented highly expanded monocot-specific groups. The Vitis RD22 group may have expanded from segmental duplications as most of its members are confined to a region in chromosome 4. The inspection of transcriptomic data revealed variable expression of BURP genes in vegetative and reproductive organs. Many genes were induced in specific tissues or by abiotic and biotic stresses. Three RD22 genes were further studied showing that they responded oppositely to ABA and to stress conditions. Our results show that the inclusion of RNA-Seq data is essential while describing gene families and improving gene annotations. Robust phylogenetic analyses including all BURP members from other sequenced species helped us redefine previous relationships that were erroneously established. This work provides additional evidence for RD22 genes serving as marker genes for different organs or stresses in grapevine.

AtRD22 and AtUSPL1, members of the plant-specific BURP domain family involved in Arabidopsis thaliana drought tolerance

Crop plants are regularly challenged by a range of environmental stresses which typically retard their growth and ultimately compromise economic yield. The stress response involves the reprogramming of approximately 4% of the transcriptome. Here, the behavior of AtRD22 and AtUSPL1, both members of the Arabidopsis thaliana BURP (BNM2, USP, RD22 and polygalacturonase isozyme) domain-containing gene family, has been characterized. Both genes are up-regulated as part of the abscisic acid (ABA) mediated moisture stress response. While AtRD22 transcript was largely restricted to the leaf, that of AtUSPL1 was more prevalent in the root. As the loss of function of either gene increased the plant's moisture stress tolerance, the implication was that their products act to suppress the drought stress response. In addition to the known involvement of AtUSPL1 in seed development, a further role in stress tolerance was demonstrated. Based on transcriptomic data and phenotype we concluded that the enhanced moisture stress tolerance of the two loss-of-function mutants is a consequence of an enhanced basal defense response.

Screening and identification of soybean seed-specific genes by using integrated bioinformatics of digital differential display, microarray, and RNA-seq data

Soybean is one of the most economically important crops in the world. Soybean seeds have abundant protein and lipid content and very high economic value. In this study, a total of 184 seed-specific genes were obtained using online microarray databases, DDD, and RNA-seq data. The reported seed-specific genes in soybean and the 184 seed-specific genes analyzed in this paper were compared. Of the screened genes, 26 were common to both previous reports and the current screening. Meanwhile, 90 of the 184 genes have homologous counterparts in Arabidopsis, among which 24 have seed-specific expression, as indicated by microarray data for Arabidopsis. Furthermore, promoter analysis showed that almost all seed-specific genes contain at least one seed specific-related element. Seed-specific element Skn-1 motif exists in most, if not all, of the seed-specific genes screened. Five genes were randomly selected from 184 soybean seed specific gene pool and their expressions were quantified using quantitative real time polymerase chain reaction (qRT-PCR) to further confirm the specificity of the screened genes. The results indicated that all five genes showed seed-specific expression. Moreover, the identification of genes with seed-specific expression screened in this study provides information valuable to the in-depth study of soybean.

Study of seed hair growth in Populus tomentosa, an important character of female floral bud development

Background

Poplar seed hair is an environmental annoyance in northern China due to its abundance and widespread airborne distribution after maturation. The morphogenesis and molecular mechanisms of its development are not well understood, and little attention has been focused on the dynamics of its development. To better understand the mechanism of poplar seed hair development, paraffin sections were used to examine the initiation and elongation of poplar seed hairs. RNA-seq technology was also employed to provide a comprehensive overview of transcriptional changes that occur during seed hair development.

Results

The placenta at the base of ovary, was identified as the origin of seed hair development, which is in sharp contrast to cotton fibers that originate from epidermal cells of the seed coat. An enlarged cell nucleus in seed hair cells was also observed, which was supported by our gene ontology enrichment analysis. The significant enriched GO term of “endoreduplication” indicated that cycles of endoreduplication, bypassing normal mitosis, is the underlying mechanisms for the maintenance of the uni-cellular structure of seed hairs. By analyzing global changes in the transcriptome, many genes regulating cell cycle, cell elongation, cell well modification were identified. Additionally, in an analysis of differential expression, cellulose synthesis and cell wall biosynthesis-related biological processes were enriched, indicating that this component of fiber structure in poplar seed hairs is consistent with what is found in cotton fibers. Differentially expressed transcription factors exhibited a stage-specific up-regulation. A dramatic down-regulation was also revealed during the mid-to-late stage of poplar seed hair development, which may point to novel mechanisms regulating cell fate determination and cell elongation.

Conclusions

This study revealed the initiation site of poplar seed hairs and also provided a comprehensive overview of transcriptome dynamics during the process of seed hair development. The high level of resolution on dynamic changes in the transcriptome provided in this study may serve as a valuable resource for developing a more complete understanding of this important biological process.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-475) contains supplementary material, which is available to authorized users.

Expression of a vacuole-localized BURP-domain protein from soybean (SALI3-2) enhances tolerance to cadmium and copper stresses

The plant-specific BURP family proteins play diverse roles in plant development and stress responses, but the function mechanism of these proteins is still poorly understood. Proteins in this family are characterized by a highly conserved BURP domain with four conserved Cys-His repeats and two other Cys, indicating that these proteins potentially interacts with metal ions. In this paper, an immobilized metal affinity chromatography (IMAC) assay showed that the soybean BURP protein SALI3-2 could bind soft transition metal ions (Cd(2+), Co(2+), Ni(2+), Zn(2+) and Cu(2+)) but not hard metal ions (Ca(2+) and Mg(2+)) in vitro. A subcellular localization analysis by confocal laser scanning microscopy revealed that the SALI3-2-GFP fusion protein was localized to the vacuoles. Physiological indexes assay showed that Sali3-2-transgenic Arabidopsis thaliana seedlings were more tolerant to Cu(2+) or Cd(2+) stresses than the wild type. An inductively coupled plasma optical emission spectrometry (ICP-OES) analysis illustrated that, compared to the wild type seedlings the Sali3-2-transgenic seedlings accumulated more cadmium or copper in the roots but less in the upper ground tissues when the seedlings were exposed to excessive CuCl2 or CdCl2 stress. Therefore, our findings suggest that the SALI3-2 protein may confer cadmium (Cd(2+)) and copper (Cu(2+)) tolerance to plants by helping plants to sequester Cd(2+) or Cu(2+) in the root and reduce the amount of heavy metals transported to the shoots.

Genome-wide analysis of soybean HD-Zip gene family and expression profiling under salinity and drought treatments

BACKGROUND

Homeodomain-leucine zipper (HD-Zip) proteins, a group of homeobox transcription factors, participate in various aspects of normal plant growth and developmental processes as well as environmental responses. To date, no overall analysis or expression profiling of the HD-Zip gene family in soybean (Glycine max) has been reported.

METHODS AND FINDINGS

An investigation of the soybean genome revealed 88 putative HD-Zip genes. These genes were classified into four subfamilies, I to IV, based on phylogenetic analysis. In each subfamily, the constituent parts of gene structure and motif were relatively conserved. A total of 87 out of 88 genes were distributed unequally on 20 chromosomes with 36 segmental duplication events, indicating that segmental duplication is important for the expansion of the HD-Zip family. Analysis of the Ka/Ks ratios showed that the duplicated genes of the HD-Zip family basically underwent purifying selection with restrictive functional divergence after the duplication events. Analysis of expression profiles showed that 80 genes differentially expressed across 14 tissues, and 59 HD-Zip genes are differentially expressed under salinity and drought stress, with 20 paralogous pairs showing nearly identical expression patterns and three paralogous pairs diversifying significantly under drought stress. Quantitative real-time RT-PCR (qRT-PCR) analysis of six paralogous pairs of 12 selected soybean HD-Zip genes under both drought and salinity stress confirmed their stress-inducible expression patterns.

CONCLUSIONS

This study presents a thorough overview of the soybean HD-Zip gene family and provides a new perspective on the evolution of this gene family. The results indicate that HD-Zip family genes may be involved in many plant responses to stress conditions. Additionally, this study provides a solid foundation for uncovering the biological roles of HD-Zip genes in soybean growth and development.

Genome-wide analysis of soybean HD-Zip gene family and expression profiling under salinity and drought treatments

BACKGROUND

Homeodomain-leucine zipper (HD-Zip) proteins, a group of homeobox transcription factors, participate in various aspects of normal plant growth and developmental processes as well as environmental responses. To date, no overall analysis or expression profiling of the HD-Zip gene family in soybean (Glycine max) has been reported.

METHODS AND FINDINGS

An investigation of the soybean genome revealed 88 putative HD-Zip genes. These genes were classified into four subfamilies, I to IV, based on phylogenetic analysis. In each subfamily, the constituent parts of gene structure and motif were relatively conserved. A total of 87 out of 88 genes were distributed unequally on 20 chromosomes with 36 segmental duplication events, indicating that segmental duplication is important for the expansion of the HD-Zip family. Analysis of the Ka/Ks ratios showed that the duplicated genes of the HD-Zip family basically underwent purifying selection with restrictive functional divergence after the duplication events. Analysis of expression profiles showed that 80 genes differentially expressed across 14 tissues, and 59 HD-Zip genes are differentially expressed under salinity and drought stress, with 20 paralogous pairs showing nearly identical expression patterns and three paralogous pairs diversifying significantly under drought stress. Quantitative real-time RT-PCR (qRT-PCR) analysis of six paralogous pairs of 12 selected soybean HD-Zip genes under both drought and salinity stress confirmed their stress-inducible expression patterns.

CONCLUSIONS

This study presents a thorough overview of the soybean HD-Zip gene family and provides a new perspective on the evolution of this gene family. The results indicate that HD-Zip family genes may be involved in many plant responses to stress conditions. Additionally, this study provides a solid foundation for uncovering the biological roles of HD-Zip genes in soybean growth and development.

Ups and downs of a transcriptional landscape shape iron deficiency associated chlorosis of the maize inbreds B73 and Mo17

BACKGROUND

Improving nutrient homeostasis is a major challenge of a sustainable maize cultivation, and cornerstone to ensure food supply for a growing world population. Although, iron constitutes an important nutrient, iron availability is limited. In this respect, iron deficiency associated chlorosis causes severe yield losses every year. Natural variation of the latter trait has yet not been addressed in maize and was therefore studied in the present analysis.

RESULTS

In this study, we i) report about the contrasting chlorosis phenotypes of the inbreds B73 and Mo17 at 10 and 300 μM iron regime, ii) identified over 400 significantly regulated transcripts (FDR < 0.05) within both inbreds at these growth conditions by deep RNA-Sequencing, iii) linked the gained knowledge with QTL information about iron deficiency related traits within the maize intermated B73 by Mo17 (IBM) population, and iv) highlighted contributing molecular pathways. In this respect, several genes within methionine salvage pathway and phytosiderophore synthesis were found to present constitutively high expression in Mo17, even under sufficient iron supply. Moreover, the same expression pattern could be observed for two putative bHLH transcription factors. In addition, a number of differentially expressed genes showed a co-localisation with QTL confidence intervals for iron deficiency related traits within the IBM population.

CONCLUSIONS

Our study highlights differential iron deficiency associated chlorosis between B73 and Mo17 and represents a valuable resource for differentially expressed genes upon iron limitation and chlorosis response. Besides identifying two putative bHLH transcription factors, we propose that methionine salvage pathway and sterol metabolism amongst others; underlie the contrasting iron deficiency related chlorosis phenotype of both inbreds. Altogether, this study emphasizes a contribution of selected genes and pathways on natural trait variation within the IBM population.

Transcription profile of soybean-root-knot nematode interaction reveals a key role of phythormones in the resistance reaction

BACKGROUND

Root-knot nematodes (RKN- Meloidogyne genus) present extensive challenges to soybean crop. The soybean line (PI 595099) is known to be resistant against specific strains and races of nematode species, thus its differential gene expression analysis can lead to a comprehensive gene expression profiling in the incompatible soybean-RKN interaction. Even though many disease resistance genes have been studied, little has been reported about phytohormone crosstalk on modulation of ROS signaling during soybean-RKN interaction.

RESULTS

Using 454 technology to explore the common aspects of resistance reaction during both parasitism and resistance phases it was verified that hormone, carbohydrate metabolism and stress related genes were consistently expressed at high levels in infected roots as compared to mock control. Most noteworthy genes include those encoding glycosyltransferases, peroxidases, auxin-responsive proteins and gibberellin-regulated genes. Our data analysis suggests the key role of glycosyltransferases, auxins and components of gibberellin signal transduction, biosynthesis and deactivation pathways in the resistance reaction and their participation in jasmonate signaling and redox homeostasis in mediating aspects of plant growth and responses to biotic stress.

CONCLUSIONS

Based on this study we suggest a reasonable model regarding to the complex mechanisms of crosstalk between plant hormones, mainly gibberellins and auxins, which can be crucial to modulate the levels of ROS in the resistance reaction to nematode invasion. The model also includes recent findings concerning to the participation of DELLA-like proteins and ROS signaling controlling plant immune or stress responses. Furthermore, this study provides a dataset of potential candidate genes involved in both nematode parasitism and resistance, which can be tested further for their role in this biological process using functional genomics approaches.

A cotton BURP domain protein interacts with α-expansin and their co-expression promotes plant growth and fruit production

Plant growth requires cell wall extension. The cotton AtRD22-Like 1 gene GhRDL1, predominately expressed in elongating fiber cells, encodes a BURP domain-containing protein. Here, we show that GhRDL1 is localized in cell wall and interacts with GhEXPA1, an α-expansin functioning in wall loosening. Transgenic cotton overexpressing GhRDL1 showed an increase in fiber length and seed mass, and an enlargement of endopleura cells of ovules. Expression of either GhRDL1 or GhEXPA1 alone in Arabidopsis led to a substantial increase in seed size; interestingly, their co-expression resulted in the increased number of siliques, the nearly doubled seed mass, and the enhanced biomass production. Cotton plants overexpressing GhRDL1 and GhEXPA1 proteins produced strikingly more fruits (bolls), leading to up to 40% higher fiber yield per plant without adverse effects on fiber quality and vegetative growth. We demonstrate that engineering cell wall protein partners has a great potential in promoting plant growth and crop yield.

Expression of an apoplast-localized BURP-domain protein from soybean (GmRD22) enhances tolerance towards abiotic stress

The BURP-domain protein family comprises a diverse group of plant-specific proteins that share a conserved BURP domain at the C terminus. However, there have been only limited studies on the functions and subcellular localization of these proteins. Members of the RD22-like subfamily are postulated to associate with stress responses due to the stress-inducible nature of some RD22-like genes. In this report, we used different transgenic systems (cells and in planta) to show that the expression of a stress-inducible RD22-like protein from soybean (GmRD22) can alleviate salinity and osmotic stress. We also performed detailed microscopic studies using both fusion proteins and immuno-electron microscopic techniques to demonstrate the apoplast localization of GmRD22, for which the BURP domain is a critical determinant of the subcellular localization. The apoplastic GmRD22 interacts with a cell wall peroxidase and the ectopic expression of GmRD22 in both transgenic Arabidopsis thaliana and transgenic rice resulted in increased lignin production when subjected to salinity stress. It is possible that GmRD22 regulates cell wall peroxidases and hence strengthens cell wall integrity under such stress conditions.

Expression pattern of drought stress marker genes in soybean roots under two water deficit systems

The study of tolerance mechanisms for drought stress in soybean is fundamental to the understanding and development of tolerant varieties. Using in silico analysis, four marker genes involved in the classical ABA-dependent and ABA-independent pathways of drought response were identified in the Glycine max genome in the present work. The expression profiles of the marker genes ERD1-like, GmaxRD20A-like, GmaxRD22-like and GmaxRD29B-like were investigated by qPCR in root samples of drought sensitive and tolerant soybean cultivars (BR 16 and Embrapa 48, respectively), submitted to water deficit conditions in hydroponic and pot-based systems. Among the four putative soybean homologs to Arabidopsis genes investigated herein, only GmaxRD29B-like was not regulated by water deficit stress. Distinct expression profiles and different induction levels were observed among the genes, as well as between the two drought-inducing systems. Our results showed contrasting gene expression responses for the GmaxRD20A-like and GmaxRD22-like genes. GmaxRD20A-like was highly induced by continuous drought acclimating conditions, whereas GmaxRD22-like responses decreased after abrupt water deprivation. GmaxERD1-like showed a different expression profile for the cultivars in each system. Conversely, GmaxRD20A-like and GmaxRD22-like genes exhibited similar expression levels in tolerant plants in both systems.

Transcript profile of the response of two soybean genotypes to potassium deficiency

The macronutrient potassium (K) is essential to plant growth and development. Crop yield potential is often affected by lack of soluble K. The molecular regulation mechanism of physiological and biochemical responses to K starvation in soybean roots and shoots is not fully understood. In the present study, two soybean varieties were subjected to low-K stress conditions: a low-K-tolerant variety (You06-71) and a low-K-sensitive variety (HengChun04-11). Eight libraries were generated for analysis: 2 genotypes ×2 tissues (roots and shoots) ×2 time periods [short term (0.5 to 12 h) and long term (3 to 12 d)]. RNA derived from the roots and shoots of these two varieties across two periods (short term and long term) were sequenced and the transcriptomes were compared using high-throughput tag-sequencing. To this end, a large number of clean tags (tags used for analysis after removal of dirty tags) corresponding to distinct tags (all types of clean tags) were identified in eight libraries (L1, You06-71-root short term; L2, HengChun04-11-root short term; L3, You06-71-shoot short term; L4, HengChun04-11-shoot short term; L5, You06-71-root long term; L6, HengChun04-11-root long term; L7, You06-71-shoot long term; L8, HengChun04-11-shoot long term). All clean tags were mapped to the available soybean (Glycine max) transcript database (http://www.soybase.org). Many genes showed substantial differences in expression across the libraries. In total, 5,440 transcripts involved in 118 KEGG pathways were either up- or down-regulated. Fifteen genes were randomly selected and their expression levels were confirmed using quantitative RT-PCR. Our results provide preliminary information on the molecular mechanism of potassium absorption and transport under low-K stress conditions in different soybean tissues.

Systematic analysis of plant-specific B3 domain-containing proteins based on the genome resources of 11 sequenced species

B3 domain-containing proteins constitute a large transcription factor superfamily. The plant-specific B3 superfamily consists of four family members, i.e., LAV (LEC2 [LEAFY COTYLEDON 2]/ABI3 [ABSCISIC ACID INSENSITIVE 3] − VAL [VP1/ABI3-LIKE]), RAV (RELATED to ABI3/VP1), ARF (AUXIN RESPONSE FACTOR) and REM (REPRODUCTIVE MERISTEM) families. The B3 superfamily plays a central role in plant life, from embryogenesis to seed maturation and dormancy. In previous research, we have characterized ARF family, member of the B3 superfamily in silico (Wang et al., Mol Biol Rep, 2011, doi:10.1007/s11033-011-0991-z). In this study, we systematically analyzed the diversity, phylogeny and evolution of B3 domain-containing proteins based on genomic resources of 11 sequenced species. A total of 865 B3 domain-containing genes were identified from 11 sequenced species through an iterative strategy. The number of B3 domain-containing genes varies not only between species but between gene families. B3 domain-containing genes are unevenly distributed in chromosomes and tend to cluster in the genome. Numerous combinations of B3 domains and their partner domains contribute to the sequences and structural diversification of the B3 superfamiy. Phylogenetic results showed that moss VAL proteins are related to LEC2/ABI3 instead of VAL proteins from higher plants. Lineage-specific expansion of ARF and REM proteins was observed. The REM family is the most diversified member among the B3 superfamily and experiences a rapid divergence during selective sweep. Based on structural and phylogenetic analysis results, two possible evolutional modes of the B3 superfamily were presented. Results presented here provide a resource for further characterization of the B3 superfamily.

Transcript profiling reveals expression differences in wild-type and glabrous soybean lines

BACKGROUND

Trichome hairs affect diverse agronomic characters such as seed weight and yield, prevent insect damage and reduce loss of water but their molecular control has not been extensively studied in soybean. Several detailed models for trichome development have been proposed for Arabidopsis thaliana, but their applicability to important crops such as cotton and soybean is not fully known.

RESULTS

Two high throughput transcript sequencing methods, Digital Gene Expression (DGE) Tag Profiling and RNA-Seq, were used to compare the transcriptional profiles in wild-type (cv. Clark standard, CS) and a mutant (cv. Clark glabrous, i.e., trichomeless or hairless, CG) soybean isoline that carries the dominant P1 allele. DGE data and RNA-Seq data were mapped to the cDNAs (Glyma models) predicted from the reference soybean genome, Williams 82. Extending the model length by 250 bp at both ends resulted in significantly more matches of authentic DGE tags indicating that many of the predicted gene models are prematurely truncated at the 5' and 3' UTRs. The genome-wide comparative study of the transcript profiles of the wild-type versus mutant line revealed a number of differentially expressed genes. One highly-expressed gene, Glyma04g35130, in wild-type soybean was of interest as it has high homology to the cotton gene GhRDL1 gene that has been identified as being involved in cotton fiber initiation and is a member of the BURP protein family. Sequence comparison of Glyma04g35130 among Williams 82 with our sequences derived from CS and CG isolines revealed various SNPs and indels including addition of one nucleotide C in the CG and insertion of ~60 bp in the third exon of CS that causes a frameshift mutation and premature truncation of peptides in both lines as compared to Williams 82.

CONCLUSION

Although not a candidate for the P1 locus, a BURP family member (Glyma04g35130) from soybean has been shown to be abundantly expressed in the CS line and very weakly expressed in the glabrous CG line. RNA-Seq and DGE data are compared and provide experimental data on the expression of predicted soybean gene models as well as an overview of the genes expressed in young shoot tips of two closely related isolines.

Screening and identification of seed-specific genes using digital differential display tools combined with microarray data from common wheat

Background

Wheat is one of the most important cereal crops for human beings, with seeds being the tissue of highly economic value. Various morphogenetic and metabolic processes are exclusively associated with seed maturation. The goal of this study was to screen and identify genes specifically expressed in the developing seed of wheat with an integrative utilization of digital differential display (DDD) and available online microarray databases.

Results

A total of 201 unigenes were identified as the results of DDD screening and microarray database searching. The expressions of 6 of these were shown to be seed-specific by qRT-PCR analysis. Further GO enrichment analysis indicated that seed-specific genes were mainly associated with defense response, response to stress, multi-organism process, pathogenesis, extracellular region, nutrient reservoir activity, enzyme inhibitor activity, antioxidant activity and oxidoreductase activity. A comparison of this set of genes with the rice (Oryza sativa) genome was also performed and approximately three-fifths of them have rice counterparts. Between the counterparts, around 63% showed similar expression patterns according to the microarray data.

Conclusions

In conclusion, the DDD screening combined with microarray data analysis is an effective strategy for the identification of seed-specific expressed genes in wheat. These seed-specific genes screened during this study will provide valuable information for further studies about the functions of these genes in wheat.

Alteration in expression of hormone-related genes in wild emmer wheat roots associated with drought adaptation mechanisms

Transcriptomic and metabolomic profiles were used to unravel drought adaptation mechanisms in wild emmer wheat (Triticum turgidum ssp. dicoccoides), the progenitor of cultivated wheat, by comparing the response to drought stress in roots of genotypes contrasting in drought tolerance. The differences between the drought resistant (R) and drought susceptible (S) genotypes were characterized mainly by shifts in expression of hormone-related genes (e.g., gibberellins, abscisic acid (ABA) and auxin), including biosynthesis, signalling and response; RNA binding; calcium (calmodulin, caleosin and annexin) and phosphatidylinositol signalling, in the R genotype. ABA content in the roots of the R genotype was higher in the well-watered treatment and increased in response to drought, while in the S genotype ABA was invariant. The metabolomic profiling revealed in the R genotype a higher accumulation of tricarboxylic acid cycle intermediates and drought-related metabolites, including glucose, trehalose, proline and glycine. The integration of transcriptomics and metabolomics results indicated that adaptation to drought included efficient regulation and signalling pathways leading to effective bio-energetic processes, carbon metabolism and cell homeostasis. In conclusion, mechanisms of drought tolerance were identified in roots of wild emmer wheat, supporting our previous studies on the potential of this genepool as a valuable source for novel candidate genes to improve drought tolerance in cultivated wheat.