Identification of the Hevea brasiliensis AP2/ERF superfamily by RNA sequencing

Chromosome-level wild Hevea brasiliensis genome provides new tools for genomic-assisted breeding and valuable loci to elevate rubber yield

The rubber tree (Hevea brasiliensis) is grown in tropical regions and is the major source of natural rubber. Using traditional breeding approaches, the latex yield has increased by sixfold in the last century. However, the underlying genetic basis of rubber yield improvement is largely unknown. Here, we present a high-quality, chromosome-level genome sequence of the wild rubber tree, the first report on selection signatures and a genome-wide association study (GWAS) of its yield traits. Population genomic analysis revealed a moderate population divergence between the Wickham clones and wild accessions. Interestingly, it is suggestive that H. brasiliensis and six relatives of the Hevea genus might belong to the same species. The selective sweep analysis found 361 obvious signatures in the domesticated clones associated with 245 genes. In a 15-year field trial, GWAS identified 155 marker-trait associations with latex yield, in which 326 candidate genes were found. Notably, six genes related to sugar transport and metabolism, and four genes related to ethylene biosynthesis and signalling are associated with latex yield. The homozygote frequencies of the causal nonsynonymous SNPs have been greatly increased under selection, which may have contributed to the fast latex yield improvement during the short domestication history. Our study provides insights into the genetic basis of the latex yield trait and has implications for genomic-assisted breeding by offering valuable resources in this new domesticated crop.

Genome-Wide Identification of the ERF Transcription Factor Family for Structure Analysis, Expression Pattern, and Response to Drought Stress in Populus alba × Populus glandulosa

The Ethylene Responsive Factor (ERF) transcription factor family is important for regulating plant growth and stress responses. Although the expression patterns of ERF family members have been reported in many plant species, their role in Populus alba × Populus glandulosa, an important model plant for forest research, remains unclear. In this study, we identified 209 PagERF transcription factors by analyzing the P. alba × P. glandulosa genome. We analyzed their amino acid sequences, molecular weight, theoretical pI (Isoelectric point), instability index, aliphatic index, grand average of hydropathicity, and subcellular localization. Most PagERFs were predicted to localize in the nucleus, with only a few PagERFs localized in the cytoplasm and nucleus. Phylogenetic analysis divided the PagERF proteins into ten groups, Class I to X, with those belonging to the same group containing similar motifs. Cis-acting elements associated with plant hormones, abiotic stress responses, and MYB binding sites were analyzed in the promoters of PagERF genes. We used transcriptome data to analyze the expression patterns of PagERF genes in different tissues of P. alba × P. glandulosa, including axillary buds, young leaves, functional leaves, cambium, xylem, and roots, and the results indicated that PagERF genes are expressed in all tissues of P. alba × P. glandulosa, especially in roots. Quantitative verification results were consistent with transcriptome data. When P. alba × P. glandulosa seedlings were treated with 6% polyethylene glycol 6000 (PEG6000), the results of RT-qRCR showed that nine PagERF genes responded to drought stress in various tissues. This study provides a new perspective on the roles of PagERF family members in regulating plant growth and development, and responses to stress in P. alba × P. glandulosa. Our study provides a theoretical basis for ERF family research in the future.

Wide-Range Portrayal of AP2/ERF Transcription Factor Family in Maize (Zea mays L.) Development and Stress Responses

The APETALA2/Ethylene-Responsive Transcriptional Factors containing conservative AP2/ERF domains constituted a plant-specific transcription factor (TF) superfamily, called AP2/ERF. The configuration of the AP2/ERF superfamily in maize has remained unresolved. In this study, we identified the 229 AP2/ERF genes in the latest (B73 RefGen_v5) maize reference genome. Phylogenetic classification of the ZmAP2/ERF family members categorized it into five clades, including 27 AP2 (APETALA2), 5 RAV (Related to ABI3/VP), 89 DREB (dehydration responsive element binding), 105 ERF (ethylene responsive factors), and a soloist. The duplication events of the paralogous genes occurred from 1.724-25.855 MYA, a key route to maize evolution. Structural analysis reveals that they have more introns and few exons. The results showed that 32 ZmAP2/ERFs regulate biotic stresses, and 24 ZmAP2/ERFs are involved in responses towards abiotic stresses. Additionally, the expression analysis showed that DREB family members are involved in plant sex determination. The real-time quantitative expression profiling of ZmAP2/ERFs in the leaves of the maize inbred line B73 under ABA, JA, salt, drought, heat, and wounding stress revealed their specific expression patterns. Conclusively, this study unveiled the evolutionary pathway of ZmAP2/ERFs and its essential role in stress and developmental processes. The generated information will be useful for stress resilience maize breeding programs.

Conventional and Omics Approaches for Understanding the Abiotic Stress Response in Cereal Crops-An Updated Overview

Cereals have evolved various tolerance mechanisms to cope with abiotic stress. Understanding the abiotic stress response mechanism of cereal crops at the molecular level offers a path to high-yielding and stress-tolerant cultivars to sustain food and nutritional security. In this regard, enormous progress has been made in the omics field in the areas of genomics, transcriptomics, and proteomics. Omics approaches generate a massive amount of data, and adequate advancements in computational tools have been achieved for effective analysis. The combination of integrated omics and bioinformatics approaches has been recognized as vital to generating insights into genome-wide stress-regulation mechanisms. In this review, we have described the self-driven drought, heat, and salt stress-responsive mechanisms that are highlighted by the integration of stress-manipulating components, including transcription factors, co-expressed genes, proteins, etc. This review also provides a comprehensive catalog of available online omics resources for cereal crops and their effective utilization. Thus, the details provided in the review will enable us to choose the appropriate tools and techniques to reduce the negative impacts and limit the failures in the intensive crop improvement study.

Genome-Wide Analysis of the Almond AP2/ERF Superfamily and Its Functional Prediction during Dormancy in Response to Freezing Stress

Simple Summary

The ethylene-responsive element (AP2/ERF) is one of the key and conserved transcription factors (TFs) in plants, and it plays a crucial role in regulating plant growth, development, and stress response. The cultivated almond in Xinjiang is often affected by short-term ultralow temperature freezing stress during the winter dormancy period, resulting in the death of large-scale almond plants. In this study, we conducted the first genome-wide analysis of the PdAP2/ERF family in almond, including protein physicochemical properties, phylogenetic relationships, motif types, gene structures, gene replication types, collinearity relationships, and cis-element types in promoter regions. We further analyzed the expression patterns of the PdAP2/ERF gene in different tissues of almond and under freezing stress at different temperatures in annual dormant branches using transcriptome data. In addition, we also analyzed the expression levels of 13 PdAP2/ERF genes in four tissues of almond and in annual dormant branches treated with freezing stress at different temperatures using fluorescence quantitative technology. This study laid the foundation for further exploring the function of the PdAP2/ERF gene in almond.

Abstract

The AP2/ERF transcription factor family is one of the largest transcription factor families in plants and plays an important role in regulating plant growth and development and the response to biotic and abiotic stresses. However, there is no report on the AP2/ERF gene family in almond (Prunus dulcis). In this study, a total of 136 PdAP2/ERF genes were identified from the almond genome, and their protein physicochemical properties were analyzed. The PdAP2/ERF members were divided into five subgroups: AP2, RAV, ERF, DREB, and Soloist. The PdAP2/ERF members in each subgroup had conserved motif types and exon/intron numbers. PdAP2/ERFS members are distributed on eight chromosomes, with 22 pairs of segmental duplications and 28 pairs of tandem duplications. We further explored the colinear relationship between almond and Arabidopsis thaliana, Oryza sativa, Malus domestica, and Prunus persicaAP2/ERF genes and their evolution. The results of cis-acting elements showed that PdAP2/ERF members are widely involved in various processes, such as growth and development, hormone regulation, and stress response. The results based on transcriptome expression patterns showed that PdAP2/ERF genes had significant tissue-specific expression characteristics and were involved in the response of annual dormant branches of almond to low-temperature freezing stress. In addition, the fluorescence quantitative relative expression results of 13 representative PdAP2/ERF genes in four tissues of ‘Wanfeng’ almond and under six low-temperature freezing treatments of annual dormant branches were consistent with the transcriptome results. It is worth noting that the fluorescence quantitative expression level showed that the PdERF24 gene was extremely significant at −30 °C, suggesting that this gene may play an important role in the response of almond dormancy to ultralow temperature freezing stress. Finally, we identified 7424 and 6971 target genes based on AP2 and ERF/DREB DNA-binding sites, respectively. The GO and KEGG enrichment results showed that these target genes play important roles in protein function and multiple pathways. In summary, we conducted bioinformatics and expression pattern studies on PdAP2/ERF genes, including 13 PdAP2/ERF genes, and performed fluorescence quantitative analysis of annual dormant shoots under different low-temperature freezing stress treatments to understand the tolerance of almond dormancy to freezing stress and suggest future improvements.

The moss-specific transcription factor PpERF24 positively modulates immunity against fungal pathogens in Physcomitrium patens

APETALA2/ethylene response factors (AP2/ERFs) transcription factors (TFs) have greatly expanded in land plants compared to algae. In angiosperms, AP2/ERFs play important regulatory functions in plant defenses against pathogens and abiotic stress by controlling the expression of target genes. In the moss Physcomitrium patens, a high number of members of the ERF family are induced during pathogen infection, suggesting that they are important regulators in bryophyte immunity. In the current study, we analyzed a P. patens pathogen-inducible ERF family member designated as PpERF24. Orthologs of PpERF24 were only found in other mosses, while they were absent in the bryophytes Marchantia polymorpha and Anthoceros agrestis, the vascular plant Selaginella moellendorffii, and angiosperms. We show that PpERF24 belongs to a moss-specific clade with distinctive amino acids features in the AP2 domain that binds to the DNA. Interestingly, all P. patens members of the PpERF24 subclade are induced by fungal pathogens. The function of PpERF24 during plant immunity was assessed by an overexpression approach and transcriptomic analysis. Overexpressing lines showed increased defenses to infection by the fungal pathogens Botrytis cinerea and Colletotrichum gloeosporioides evidenced by reduced cellular damage and fungal biomass compared to wild-type plants. Transcriptomic and RT-qPCR analysis revealed that PpERF24 positively regulates the expression levels of defense genes involved in transcriptional regulation, phenylpropanoid and jasmonate pathways, oxidative burst and pathogenesis-related (PR) genes. These findings give novel insights into potential mechanism by which PpERF24 increases plant defenses against several pathogens by regulating important players in plant immunity.

Reprogramming of the Hevea brasiliensis Epigenome and Transcriptome in Response to Cold Stress

Low temperature is a key factor limiting the rubber plantation extending to high latitude area. Previous work has shown that cold-induced DNA demethylation was coordinated with the expression of cold-responsive (COR) genes in Hevea brasiliensis. In this work, reduced representation bisulphite sequencing analysis of H. brasiliensis showed that cold treatment induced global genomic DNA demethylation and altered the sequence contexts of methylated cytosines, but the levels of mCG methylation in transposable elements were slightly enhanced by cold treatment. Integrated analysis of the DNA methylome and transcriptome revealed 400 genes whose expression correlated with altered DNA methylation. DNA demethylation in the upstream region of gene seems to correlate with higher gene expression, whereas demethylation in the gene body has less association. Our results suggest that cold treatment globally change the genomic DNA methylation status of the rubber tree, which might coordinate reprogramming of the transcriptome.

Transcriptome profiling in susceptible and tolerant rubber tree clones in response to cassiicolin Cas1, a necrotrophic effector from Corynespora cassiicola

Corynespora cassiicola, a fungal plant pathogen with a large host range, causes important damages in rubber tree (Hevea brasiliensis), in Asia and Africa. A small secreted protein named cassiicolin was previously identified as a necrotrophic effector required for the virulence of C. cassiicola in specific rubber tree clones. The objective of this study was to decipher the cassiicolin-mediated molecular mechanisms involved in this compatible interaction. We comparatively analyzed the RNA-Seq transcriptomic profiles of leaves treated or not with the purified cassiicolin Cas1, in two rubber clones: PB260 (susceptible) and RRIM600 (tolerant). The reads were mapped against a synthetic transcriptome composed of all available transcriptomic references from the two clones. Genes differentially expressed in response to cassiicolin Cas1 were identified, in each clone, at two different time-points. After de novo annotation of the synthetic transcriptome, we analyzed GO enrichment of the differentially expressed genes in order to elucidate the main functional pathways impacted by cassiicolin. Cassiicolin induced qualitatively similar transcriptional modifications in both the susceptible and the tolerant clones, with a strong negative impact on photosynthesis, and the activation of defense responses via redox signaling, production of pathogenesis-related protein, or activation of the secondary metabolism. In the tolerant clone, transcriptional reprogramming occurred earlier but remained moderate. By contrast, the susceptible clone displayed a late but huge transcriptional burst, characterized by massive induction of phosphorylation events and all the features of a hypersensitive response. These results confirm that cassiicolin Cas1 is a necrotrophic effector triggering a hypersensitive response in susceptible rubber clones, in agreement with the necrotrophic-effector-triggered susceptibility model.

Ethylene response factors regulate expression of HbSUT3, the sucrose influx carrier in laticifers of Hevea brasiliensis

Natural rubber is an important industrial raw material and is commercially produced by rubber trees (Hevea brasiliensis). The sucrose transporter HbSUT3 plays an essential role in rubber production. Its expression in latex (cytoplasm of rubber-producing laticifers) is induced by bark treatment with Ethrel, an ethylene releaser, and the inducing effect correlates well with Ethrel-stimulated rubber yield increase. However, the mechanisms of ethylene induction on HbSUT3 expression are not known. Here, five Ethylene Response Factor (ERF) genes were identified from the cDNA library of Hevea latex by yeast one-hybrid screening with the promoter of HbSUT3 gene as bait. As revealed in a tobacco (Nicotiana tabacum) protoplast transient expression system, these HbERFs were mainly localized in the nucleus and four of them exhibited apparent transactivation activity. Of the five HbERF genes, HbERF-IXc4 was the most frequently screened in yeast one-hybrid, accounting for 65% of the ERF clones obtained. Moreover, among the five HbERFs, HbERF-IXc4 showed the strongest transactivation capacity when expressed in tobacco protoplast, the highest transcript abundance in latex and a close expressional correlation with its target gene, HbSUT3, in response to the Ethrel treatment. Taken together, our results indicate that ERFs, especially HbERF-IXc4, are critically involved in the activation of HbSUT3 expression in latex after Ethrel treatment on Hevea bark, and thus the stimulated latex yield.

Genome-wide analysis of AP2/ERF superfamily in lotus (Nelumbo nucifera) and the association between NnADAP and rhizome morphology

Background

The AP2/ERF family is widely present in plants and plays a crucial regulatory role in plant growth and development. As an essential aquatic horticultural model plant, lotus has an increasingly prominent economic and research value.

Results

We have identified and analysed the AP2/ERF gene family in the lotus. Initially, 121 AP2/ERF family genes were identified. By analysing their gene distribution and protein structure, and their expression patterns during the development of lotus rhizome, combined with previous studies, we obtained an SNP (megascaffold_20:3578539) associated with lotus rhizome phenotype. This SNP was in the NnADAP gene of the AP2 subfamily, and the changes in SNP (C/T) caused amino acid conversion (proline/leucine). We constructed a population of 95 lotus varieties for SNP verification. Through population typing experiments, we found that the group with SNP CC had significantly larger lotus rhizome and higher soluble sugar content among the population.

Conclusions

In conclusion, we speculate that the alteration of the SNP in the NnADAP can affect the size and sugar content of the lotus rhizome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07473-w.

Identification, classification, and characterization of AP2/ERF superfamily genes in Masson pine (Pinus massoniana Lamb.)

Transcription factors (TFs) play crucial regulatory roles in controlling the expression of the target genes in plants. APETALA2/Ethylene-responsive factors (AP2/ERF) are part of a large superfamily of plant-specific TFs whose members are involved in the control of plant metabolism, development and responses to various biotic and abiotic stresses. However, the AP2/ERF superfamily has not been identified systematically in Masson pine (Pinus massoniana), which is one of the most important conifer in southern China. Therefore, we performed systematic identification of the AP2/ERF superfamily using transcriptome sequencing data from Masson pine. In the current study, we obtained 88 members of the AP2/ERF superfamily. All PmAP2/ERF members could be classified into 3 main families, AP2 (7 members), RAV (7 members), ERF (73 members) families, and a soloist protein. Subcellular localization assays suggested that two members of PmAP2/ERF were nuclear proteins. Based on pine wood nematode (PWN) inoculated transcriptome and qPCR analysis, we found that many members of PmAP2/ERF could respond to PWN inoculation and PWN related treatment conditions in vitro. In general, members of the AP2/ERF superfamily play an important role in the response of Masson pine responds to PWN. Furthermore, the roles of the AP2/ERF superfamily in other physiological activities of Masson pine remain to be further studied.

Structural and Functional Annotation of Transposable Elements Revealed a Potential Regulation of Genes Involved in Rubber Biosynthesis by TE-Derived siRNA Interference in Hevea brasiliensis

The natural rubber biosynthetic pathway is well described in Hevea, although the final stages of rubber elongation are still poorly understood. Small Rubber Particle Proteins and Rubber Elongation Factors (SRPPs and REFs) are proteins with major function in rubber particle formation and stabilization. Their corresponding genes are clustered on a scaffold1222 of the reference genomic sequence of the Hevea brasiliensis genome. Apart from gene expression by transcriptomic analyses, to date, no deep analyses have been carried out for the genomic environment of SRPPs and REFs loci. By integrative analyses on transposable element annotation, small RNAs production and gene expression, we analysed their role in the control of the transcription of rubber biosynthetic genes. The first in-depth annotation of TEs (Transposable Elements) and their capacity to produce TE-derived siRNAs (small interfering RNAs) is presented, only possible in the Hevea brasiliensis clone PB 260 for which all data are available. We observed that 11% of genes are located near TEs and their presence may interfere in their transcription at both genetic and epigenetic level. We hypothesized that the genomic environment of rubber biosynthesis genes has been shaped by TE and TE-derived siRNAs with possible transcriptional interference on their gene expression. We discussed possible functionalization of TEs as enhancers and as donors of alternative transcription start sites in promoter sequences, possibly through the modelling of genetic and epigenetic landscapes.

Post-transcriptional regulation of several biological processes involved in latex production in Hevea brasiliensis

Background

Small RNAs modulate plant gene expression at both the transcriptional and post-transcriptional level, mostly through the induction of either targeted DNA methylation or transcript cleavage, respectively. Small RNA networks are involved in specific plant developmental processes, in signaling pathways triggered by various abiotic stresses and in interactions between the plant and viral and non-viral pathogens. They are also involved in silencing maintenance of transposable elements and endogenous viral elements. Alteration in small RNA production in response to various environmental stresses can affect all the above-mentioned processes. In rubber trees, changes observed in small RNA populations in response to trees affected by tapping panel dryness, in comparison to healthy ones, suggest a shift from a transcriptional to a post-transcriptional regulatory pathway. This is the first attempt to characterise small RNAs involved in post-transcriptional silencing and their target transcripts in Hevea.

Methods

Genes producing microRNAs (MIR genes) and loci producing trans-activated small interfering RNA (ta-siRNA) were identified in the clone PB 260 re-sequenced genome. Degradome libraries were constructed with a pool of total RNA from six different Hevea tissues in stressed and non-stressed plants. The analysis of cleaved RNA data, associated with genomics and transcriptomics data, led to the identification of transcripts that are affected by 20–22 nt small RNA-mediated post-transcriptional regulation. A detailed analysis was carried out on gene families related to latex production and in response to growth regulators.

Results

Compared to other tissues, latex cells had a higher proportion of transcript cleavage activity mediated by miRNAs and ta-siRNAs. Post-transcriptional regulation was also observed at each step of the natural rubber biosynthesis pathway. Among the genes involved in the miRNA biogenesis pathway, our analyses showed that all of them are expressed in latex. Using phylogenetic analyses, we show that both the Argonaute and Dicer-like gene families recently underwent expansion. Overall, our study underlines the fact that important biological pathways, including hormonal signalling and rubber biosynthesis, are subject to post-transcriptional silencing in laticifers.

In silico mining and functional analysis of AP2/ERF gene in Withania somnifera

Withania somnifera owing to its strong and remarkable stress tolerance property is a reliable candidate for the determination of genes involved in mechanism of adaption/tolerance of various stress conditions. 187 AP2/ERF gene related transcripts (GRTs) were identified during comprehensive search in W. somnifera transcriptome repertoire. Major hits in homology search were observed from the model plant Arabidopsis and members of Solanaceae family. Cloning, expression analysis of the gene and genetic transient transformation with the gene (WsAP2) were performed to predict its functional role in planta. Enhanced expression of some of the pathway genes for terpenoid biosynthesis was observed in transformed tissues in comparison to the control tissues. It is speculated that WsAP2 gene crucially regulates the expression of GGPPS gene in addition to the regulation of other important genes of terpenoid pathway via induction of expression of other genes such as HMGR, CAS, DXS and DXR. To the best of our knowledge, this is the first report representing detailed study of AP2/ERF gene family in W. somnifera. It is also suggested from the study that gene might have role in eliciting responses to combat stress and attribute the strong stress tolerant property associated with the plant.

Structural variation, functional differentiation and expression characteristics of the AP2/ERF gene family and its response to cold stress and methyl jasmonate in Panax ginseng C.A. Meyer

The APETALA2/Ethylene Responsive Factor (AP2/ERF) gene family has been shown to play a crucial role in plant growth and development, stress responses and secondary metabolite biosynthesis. Nevertheless, little is known about the gene family in ginseng (Panax ginseng C.A. Meyer), an important medicinal herb in Asia and North America. Here, we report the systematic analysis of the gene family in ginseng using several transcriptomic databases. A total of 189 putative AP2/ERF genes, defined as PgERF001 through PgERF189, were identified and these PgERF genes were spliced into 397 transcripts. The 93 PgERF genes that have complete AP2 domains in open reading frame were classified into five subfamilies, DREB, ERF, AP2, RAV and Soloist. The DREB subfamily and ERF subfamily were further clustered into four and six groups, respectively, compared to the 12 groups of these subfamilies found in Arabidopsis thaliana. Gene ontology categorized these 397 transcripts of the 189 PgERF genes into eight functional subcategories, suggesting their functional differentiation, and they have been especially enriched for the subcategory of nucleic acid binding transcription factor activity. The expression activity and networks of the 397 PgERF transcripts have substantially diversified across tissues, developmental stages and genotypes. The expressions of the PgERF genes also significantly varied, when ginseng was subjected to cold stress, as tested using six PgERF genes, PgERF073, PgERF079, PgERF110, PgERF115, PgERF120 and PgERF128, randomly selected from the DREB subfamily. This result suggests that the DREB subfamily genes play an important role in plant response to cold stress. Finally, we studied the responses of the PgERF genes to methyl jasmonate (MeJA). We found that 288 (72.5%) of the 397 PgERF gene transcripts responded to the MeJA treatment, with 136 up-regulated and 152 down-regulated, indicating that most members of the PgERF gene family are responsive to MeJA. These results, therefore, provide new resources and knowledge necessary for family-wide functional analysis of the PgERF genes in ginseng and related species.

Single-nucleotide polymorphism markers within MVA and MEP pathways among Hevea brasiliensis clones through transcriptomic analysis

In this study, RNA sequencing of several Hevea brasiliensis clones grown in Malaysia with different annual rubber production yields and disease resistance was performed on the Illumina platform. A total of 29,862,548 reads were generated, resulting in 101,269 assembled transcripts that were used as the reference transcripts. A similarity search against the non-redundant (nr) protein databases presented 83,771 (83%) positive BLASTx hits. The transcriptome was annotated using gene ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the Pfam database. A search for putative molecular markers was performed to identify single-nucleotide polymorphisms (SNPs). Overall, 3,210,629 SNPs were detected and a total of 1314 SNPs associated with the genes involved in MVA and MEP pathways were identified. A total of 176 SNP primer pairs were designed from sequences that were related to the MVA and MEP pathways. The transcriptome of RRIM 3001 and RRIM 712 were subjected to pairwise comparison and the results revealed that there were 1262 significantly differentially expressed genes unique to RRIM 3001, 1499 significantly differentially expressed genes unique to RRIM 712 and several genes related to the MVA and MEP pathways such as AACT, HMGS, PMK, MVD, DXS and HDS were included. The results will facilitate the characterization of H. brasiliensis transcriptomes and the development of a new set of molecular markers in the form of SNPs from transcriptome assembly for the genotype identification of various rubber varieties with superior traits in Malaysia.

The expanding roles of APETALA2/Ethylene Responsive Factors and their potential applications in crop improvement

Understanding the molecular basis of the gene-regulatory networks underlying agronomic traits or plant responses to abiotic/biotic stresses is very important for crop improvement. In this context, transcription factors, which either singularly or in conjugation directly control the expression of many target genes, are suitable candidates for improving agronomic traits via genetic engineering. In this regard, members of one of the largest class of plant-specific APETALA2/Ethylene Response Factor (AP2/ERF) superfamily, which is implicated in various aspects of development and plant stress adaptation responses, are considered high-value targets for crop improvement. Besides their long-known regulatory roles in mediating plant responses to abiotic stresses such as drought and submergence, the novel roles of AP2/ERFs during fruit ripening or secondary metabolites production have also recently emerged. The astounding functional plasticity of AP2/ERF members is considered to be achieved by their interplay with other regulatory networks and signalling pathways. In this review, we have integrated the recently accumulated evidence from functional genomics studies and described their newly emerged functions in plants. The key structural features of AP2/ERF proteins and the modes of their action are briefly summarized. The importance of AP2/ERFs in plant development and stress responses and a summary of the event of their successful applications in crop improvement programs are also provided. Altogether, we envisage that the synthesized information presented in this review will be useful to design effective strategies for improving agronomic traits in crop plants.

Identification of the targets of HbEIN3/EILs in genomic wide in Hevea brasiliensis

EIN3/EILs are key regulators in ET signaling pathway. In this work, 4 members of EIN3/EILs of Hevea brasiliensis (HbEIN3/EILs) showed interaction with two F box proteins, HbEBF1 and HbEBF2. HbEIN3 located in nucleus and exhibited strong transcriptional activity. HbEIN3 was induced by ET treatment in C-serum, but not in B-serum of latex. HbEIN3/EILs bound to G-box cis-element. To globally search the potential targets of HbEIN3/EILs, genomic sequences of H. brasiliensis was re-annotated and an HCES (Hevea Cis-Elements Scanning) program was developed ( www.h-brasiliensis.com ). HCES scanning results showed that ET- and JA- responsive cis-elements distribute overlapping in gene promoters. 3146 genes containing G-box in promoters are potential targets of HbEIN3, including 41 genes involved in biosynthesis and drainage of latex, of which 7 rate-limiting genes of latex production were regulated by both ET and JA, suggesting that ET and JA signaling pathways coordinated the latex biosynthesis and drainage in H. brasiliensis. Abbreviations: ABRE: ABA responsive elements; bHLH: basic helix-loop-helix; COG: Orthologous Groups; DRE: dehydration response element; ERE: ethylene responsive element; ET: Ethylene; GO: Gene Ontology; HCES: Hevea Cis-Elements Scanning; JA: jasmonates; JRE: Jasmonate-responsive element; KEGG: Kyoto Encyclopedia of Genes and Genomes; NR: non-redundant database; PLACE: Plant Cis-acting Regulatory DNA Elements; qRT-PCR: quantitative real-time RT-PCR.

Identification, Classification, and Functional Analysis of AP2/ERF Family Genes in the Desert Moss Bryum argenteum

Bryum argenteum is a desert moss which shows tolerance to the desert environment and is emerging as a good plant material for identification of stress-related genes. AP2/ERF transcription factor family plays important roles in plant responses to biotic and abiotic stresses. AP2/ERF genes have been identified and extensively studied in many plants, while they are rarely studied in moss. In the present study, we identified 83 AP2/ERF genes based on the comprehensive dehydrationrehydration transcriptomic atlas of B. argenteum. BaAP2/ERF genes can be classified into five families, including 11 AP2s, 43 DREBs, 26 ERFs, 1 RAV, and 2 Soloists. RNA-seq data showed that 83 BaAP2/ERFs exhibited elevated transcript abundances during dehydration⁻rehydration process. We used RT-qPCR to validate the expression profiles of 12 representative BaAP2/ERFs and confirmed the expression trends using RNA-seq data. Eight out of 12 BaAP2/ERFs demonstrated transactivation activities. Seven BaAP2/ERFs enhanced salt and osmotic stress tolerances of yeast. This is the first study to provide detailed information on the identification, classification, and functional analysis of the AP2/ERFs in B. argenteum. This study will lay the foundation for the further functional analysis of these genes in plants, as well as provide greater insights into the molecular mechanisms of abiotic stress tolerance of B. argenteum.

Identification of browning-related microRNAs and their targets reveals complex miRNA-mediated browning regulatory networks in Luffa cylindrica

As a non-coding and endogenous small RNA, MicroRNA (miRNA) takes a vital regulatory role in plant growth and development. Long-term storage and processing of many fruits and vegetables, including Luffa, are subject to influences from browning, a common post-harvest problem that adversely affects flavor, aroma, and nutritional value. The browning regulatory networks mediated by miRNA, however, remain largely unexplored. For a systematic identification of browning-responsive miRNAs and the targets, we built two RNA libraries from Luffa pulps of near-isogenic line, with resistant and sensitive browning characteristics respectively, and then sequenced them using Solexa high-throughput technology. We consequently identified 179 known miRNAs that represent 17 non-conserved miRNA families and 24 conserved families, as well as 84 potential novel miRNAs, among which 16 miRNAs (eight known and eight novel miRNAs) were found to exhibit significant differential expressions and were thus identified as browning-related miRNAs. We then studied those browning-responsive miRNAs and the corresponding targets with RT-qPCR and finally validated their expression patterns. The results revealed that the expression patterns are specific to plant development stages and the miRNAs are identified with 39 target transcripts, which involve in plant development, defense response, transcriptional regulation, and signal transduction. After characterizing these miRNAs and their targets, we propose a browning regulatory network model of miRNA-mediatation in this paper. The findings of the work are helpful for the understanding of miRNA-mediated regulatory mechanisms of browning in Luffa, and will facilitate genetic improvement of pulp characteristics in Luffa.

Transcriptome analysis in Hevea brasiliensis latex revealed changes in hormone signalling pathways during ethephon stimulation and consequent Tapping Panel Dryness

Tapping Panel Dryness (TPD) affects latex production in Hevea brasiliensis. This physiological syndrome involves the agglutination of rubber particles, which leads to partial or complete cessation of latex flow. Latex harvesting consists in tapping soft bark. Ethephon can be applied to stimulate latex flow and its regeneration in laticifers. Several studies have reported transcriptome changes in bark tissues. This study is the first report on deep RNA sequencing of latex to compare the effect of ethephon stimulation and TPD severity. Trees were carefully selected for paired-end sequencing using an Illumina HiSeq 2000. In all, 43 to 60 million reads were sequenced for each treatment in three biological replicates (slight TPD trees without ethephon stimulation, and slight and severe TPD trees with ethephon treatment). Differentially expressed genes were identified and annotated, giving 8,111 and 728 in response to ethephon in slight TPD trees and in ethephon-induced severe TPD trees, respectively. A biological network of responses to ethephon and TPD highlighted the major influence of metabolic processes and the response to stimulus, especially wounding and jasmonate depression in TPD-affected trees induced by ethephon stimulation.

RAV transcription factors are essential for disease resistance against cassava bacterial blight via activation of melatonin biosynthesis genes

With 1 AP2 domain and 1 B3 domain, 7 MeRAVs in apetala2/ethylene response factor (AP2/ERF) gene family have been identified in cassava. However, the in vivo roles of these remain unknown. Gene expression assays showed that the transcripts of MeRAVs were commonly regulated after Xanthomonas axonopodis pv manihotis (Xam) and MeRAVs were specifically located in plant cell nuclei. Through virus-induced gene silencing (VIGS) in cassava, we found that MeRAV1 and MeRAV2 are essential for plant disease resistance against cassava bacterial blight, as shown by the bacterial propagation of Xam in plant leaves. Through VIGS in cassava leaves and overexpression in cassava leave protoplasts, we found that MeRAV1 and MeRAV2 positively regulated melatonin biosynthesis genes and the endogenous melatonin level. Further investigation showed that MeRAV1 and MeRAV2 are direct transcriptional activators of 3 melatonin biosynthesis genes in cassava, as evidenced by chromatin immunoprecipitation-PCR in cassava leaf protoplasts and electrophoretic mobility shift assay. Moreover, cassava melatonin biosynthesis genes also positively regulated plant disease resistance. Taken together, this study identified MeRAV1 and MeRAV2 as common and upstream transcription factors of melatonin synthesis genes in cassava and revealed a model of MeRAV1 and MeRAV2-melatonin biosynthesis genes-melatonin level in plant disease resistance against cassava bacterial blight.

Overexpression of Hevea brasiliensis ethylene response factor HbERF-IXc5 enhances growth and tolerance to abiotic stress and affects laticifer differentiation

Ethylene response factor 1 (ERF1) is an essential integrator of the jasmonate and ethylene signalling pathways coordinating a large number of genes involved in plant defences. Its orthologue in Hevea brasiliensis, HbERF-IXc5, has been assumed to play a major role in laticifer metabolism and tolerance to harvesting stress for better latex production. This study sets out to establish and characterize rubber transgenic lines overexpressing HbERF-IXc5. Overexpression of HbERF-IXc5 dramatically enhanced plant growth and enabled plants to maintain some ecophysiological parameters in response to abiotic stress such as water deficit, cold and salt treatments. This study revealed that HbERF-IXc5 has rubber-specific functions compared to Arabidopsis ERF1 as transgenic plants overexpressing HbERF-IXc5 accumulated more starch and differentiated more latex cells at the histological level. The role of HbERF-IXc5 in driving the expression of some target genes involved in laticifer differentiation is discussed.

Papain-like cysteine protease encoding genes in rubber (Hevea brasiliensis): comparative genomics, phylogenetic, and transcriptional profiling analysis

43 HbPLCPs representing nine subfamilies or 20 orthologous groups were found in rubber, where paralogs were resulted from the recent WGD and local duplication. Several senescence-associated genes were also identified. Papain-like cysteine proteases (PLCPs) comprise a large family of proteolytic enzymes involved in plant growth and development, seed germination, organ senescence, immunity, and stress response. Despite their importance and the extensive research in the model plant Arabidopsis thaliana, little information is available on rubber tree (Hevea brasiliensis), a rubber-producing plant of the Euphorbiaceae family. This study performed a genome-wide identification of PLCP family genes in rubber, resulting in a relatively high number of 43 members. The phylogenetic analysis assigned these genes into nine subfamilies, i.e., RD21 (6), CEP (4), XCP (4), XBCP3 (2), THI (1), SAG12 (18), RD19 (4), ALP (2), and CTB (2). Most of them were shown to have orthologs in Arabidopsis; however, several members in SAG12, CEP and XBCP3 subfamilies form new groups as observed in other core eudicots such as Manihot esculenta, Ricinus communis, Populus trichocarpa, and Vitis vinifera. Based on an expert sequence comparison, 20 orthologous groups (OGs) were proposed for core eudicots, and rubber paralogs were shown to be resulted from the recent whole-genome duplication (WGD) as well as local duplication. Transcriptional profiling showed distinct expression pattern of different members across various tissues, e.g., root, leaf, bark, laticifer, flower, and seed. By using the senescence-specific HbSAG12H1 as the indicator, the transcriptome of senescent rubber leaves was deeply sequenced and several senescence-associated PLCP genes were identified. Results obtained from this study provide valuable information for future functional analysis and utilization of PLCP genes in Hevea and other species.

Identification of novel miRNAs and miRNA expression profiling in embryogenic tissues of Picea balfouriana treated by 6-benzylaminopurine

Here, we compared miRNA expression profiles in embryonic cell cultures of the conifer Picea balfouriana following application of the synthetic cytokinin 6-benzylaminopurine (6-BAP). We used next-generation sequencing to analyze three libraries of small RNAs from the treated embryogenic cell cultures and generated 24,000,000 raw reads from each of the libraries. Over 70 differentially regulated micro RNA (miRNA) families (≥2 fold change in expression) were identified between pairs of treatments. A quantitative analysis showed that miR3633 and miR1026 were upregulated in tissues with the highest embryogenic ability. These two miRNAs were predicted to target genes encoding receptor-like protein kinase and GAMYB transcription factors, respectively. In one library, miR1160, miR5638, miR1315, and miR5225 were downregulated. These four miRNAs were predicted to target genes encoding APETALA2, calmodulin-binding protein, and calcium-dependent protein kinase transcription factors. The expression patterns of the miRNAs and their targets were negatively correlated. Approximately 181 potentially novel P. balfouriana miRNAs were predicted from the three libraries, and seven were validated during the quantitative analysis. This study is the first report of differential miRNA regulation in tissues treated with 6-BAP during somatic embryogenesis. The differentially expressed miRNAs will be of value for investigating the mechanisms of embryogenic processes that are responsive to 6-BAP in P. balfouriana.

Transcriptome-Wide Identification, Classification, and Characterization of AP2/ERF Family Genes in the Desert Moss Syntrichia caninervis

APETALA2/Ethylene Responsive Factor (AP2/ERF) is a large family of plant transcription factors which play important roles in the control of plant metabolism and development as well as responses to various biotic and abiotic stresses. The desert moss Syntrichia caninervis, due to its robust and comprehensive stress tolerance, is a promising organism for the identification of stress-related genes. Using S. caninervis transcriptome data, 80 AP2/ERF unigenes were identified by HMM modeling and BLASTP searching. Based on the number of AP2 domains, multiple sequence alignment, motif analysis, and gene tree construction, ScAP2/ERF genes were classified into three main subfamilies (including 5 AP2 gene members, 72 ERF gene members, and 1 RAV member) and two Soloist members. We found that the ratio for each subfamily was constant between S. caninervis and the model moss Physcomitrella patens, however, as compared to the angiosperm Arabidopsis, the percentage of ERF subfamily members in both moss species were greatly expanded, while the members of the AP2 and RAV subfamilies were reduced accordingly. The amino acid composition of the AP2 domain of ScAP2/ERFs was conserved as compared with Arabidopsis. Interestingly, most of the identified DREB genes in S. caninervis belonged to the A-5 group which play important roles in stress responses and are rarely reported in the literature. Expression profile analysis of ScDREB genes showed different gene expression patterns under dehydration and rehydration; the majority of ScDREB genes demonstrated a stronger response to dehydration relative to rehydration indicating that ScDREB may play an important role in dehydrated moss tissues. To our knowledge, this is the first study to detail the identification and characterization of the AP2/ERF gene family in a desert moss. Further, this study will lay the foundation for further functional analysis of these genes, provide greater insight to the stress tolerance mechanisms in S. caninervis and provide a reference for AP2/ERF gene family classification in other moss species.

Survey of the rubber tree genome reveals a high number of cysteine protease-encoding genes homologous to Arabidopsis SAG12

Arabidopsis thaliana SAG12, a senescence-specific gene encoding a cysteine protease, is widely used as a molecular marker for the study of leaf senescence. To date, its potential orthologues have been isolated from several plant species such as Brassica napus and Nicotiana tabacum. However, little information is available in rubber tree (Hevea brasiliensis), a rubber-producing plant of the Euphorbiaceae family. This study presents the identification of SAG12-like genes from the rubber tree genome. Results showed that an unexpected high number of 17 rubber orthologues with a single intron were found, contrasting the single copy with two introns in Arabidopsis. The gene expansion was also observed in another two Euphorbiaceae plants, castor bean (Ricinus communis) and physic nut (Jatropha curcas), both of which contain 8 orthologues. In accordance with no occurrence of recent whole-genome duplication (WGD) events, most duplicates in castor and physic nut were resulted from tandem duplications. In contrast, the duplicated HbSAG12H genes were derived from tandem duplications as well as the recent WGD. Expression analysis showed that most HbSAG12H genes were lowly expressed in examined tissues except for root and male flower. Furthermore, HbSAG12H1 exhibits a strictly senescence-associated expression pattern in rubber tree leaves, and thus can be used as a marker gene for the study of senescence mechanism in Hevea.

Divergent Expression Patterns in Two Vernicia Species Revealed the Potential Role of the Hub Gene VmAP2/ERF036 in Resistance to Fusarium oxysporum in Vernicia montana

Tung oil tree (Vernicia fordii) is a promising industrial oil crop; however, this tree is highly susceptible to Fusarium wilt disease. Conversely, Vernicia montana is resistant to the pathogen. The APETALA2/ethylene-responsive element binding factor (AP2/ERF) transcription factor superfamily has been reported to play a significant role in resistance to Fusarium oxysporum. In this study, comprehensive analysis identified 75 and 81 putative Vf/VmAP2/ERF transcription factor-encoding genes in V. fordii and V. montana, respectively, which were divided into AP2, ERF, related to ABI3 and VP1 (RAV) and Soloist families. After F. oxysporum infection, a majority of AP2/ERF superfamily genes showed strong patterns of repression in both V. fordii and V. montana. We then identified 53 pairs of one-to-one orthologs in V. fordii and V. montana, with most pairs of orthologous genes exhibiting similar expression in response to the pathogen. Further investigation of Vf/VmAP2/ERF gene expression in plant tissues indicated that the pairs of genes with different expression patterns in response to F. oxysporum tended to exhibit different tissue profiles in the two species. In addition, VmAP2/ERF036, showing the strongest interactions with 666 genes, was identified as a core hub gene mediating resistance. Moreover, qRT-PCR results indicated VmAP2/ERF036 showed repressed expression while its orthologous gene VfAP2/ERF036 had the opposite expression pattern during pathogen infection. Overall, comparative analysis of the Vf/VmAP2/ERF superfamily and indication of a potential hub resistance gene in resistant and susceptible Vernicia species provides valuable information for understanding the molecular basis and selection of essential functional genes for V. fordii genetic engineering to control Fusarium wilt disease.

Laser microdissection of tomato fruit cell and tissue types for transcriptome profiling

This protocol enables transcriptome profiling of specific cell or tissue types that are isolated from tomato using laser microdissection (LM). To prepare tissue for LM, fruit samples are first fixed in optimal cutting temperature (OCT) medium and frozen in molds. The tissue is then sectioned using a cryostat before being dissected using an LM instrument. The RNAs contained in the harvested cells are purified and subjected to two rounds of amplification to yield sufficient quantities of RNA to generate cDNA libraries. Unlike several other techniques that are used to isolate specific cell types, LM has the advantage of being readily applied to any plant species without having to generate transgenic plants. Using the protocols described here, LM-mediated cell-type transcriptomic analysis of two samples requires ∼8 d from tissue harvest to RNA sequencing (RNA-seq), whereas each additional sample, up to a total of 12 samples, requires ∼1 additional day for the LM step. RNA obtained using this method has been successfully used for deep-coverage transcriptome profiling, which is a particularly effective strategy for identifying genes that are differentially expressed between cell or tissue types.

RAP2.4a Is Transported through the Phloem to Regulate Cold and Heat Tolerance in Papaya Tree (Carica papaya cv. Maradol): Implications for Protection Against Abiotic Stress

Plants respond to stress through metabolic and morphological changes that increase their ability to survive and grow. To this end, several transcription factor families are responsible for transmitting the signals that are required for these changes. Here, we studied the transcription factor superfamily AP2/ERF, particularly, RAP2.4 from Carica papaya cv. Maradol. We isolated four genes (CpRap2.4a, CpRAap2.4b, CpRap2.1 and CpRap2.10), and an in silico analysis showed that the four genes encode proteins that contain a conserved APETALA2 (AP2) domain located within group I and II transcription factors of the AP2/ERF superfamily. Semiquantitative PCR experiments indicated that each CpRap2 gene is differentially expressed under stress conditions, such as extreme temperatures. Moreover, genetic transformants of tobacco plants overexpressing CpRap2.4a and CpRap2.4b genes show a high level of tolerance to cold and heat stress compared to non-transformed plants. Confocal microscopy analysis of tobacco transgenic plants showed that CpRAP2.4a and CpRAP2.4b proteins were mainly localized to the nuclei of cells from the leaves and roots and also in the sieve elements. Moreover, the movement of CpRap2.4a RNA in tobacco grafting was analyzed. Our results indicate that CpRap2.4a and CpRap2.4b RNA in the papaya tree have a functional role in the response to stress conditions such as exposure to extreme temperatures via direct translation outside the parental RNA cell.

Expansion and Functional Divergence of AP2 Group Genes in Spermatophytes Determined by Molecular Evolution and Arabidopsis Mutant Analysis

The APETALA2 (AP2) genes represent the AP2 group within a large group of DNA-binding proteins called AP2/EREBP. The AP2 gene is functional and necessary for flower development, stem cell maintenance, and seed development, whereas the other members of AP2 group redundantly affect flowering time. Here we study the phylogeny of AP2 group genes in spermatophytes. Spermatophyte AP2 group genes can be classified into AP2 and TOE types, six clades, and we found that the AP2 group homologs in gymnosperms belong to the AP2 type, whereas TOE types are absent, which indicates the AP2 type gene are more ancient and TOE type was split out of AP2 type and losing the major function. In Brassicaceae, the expansion of AP2 and TOE type lead to the gene number of AP2 group were up to six. Purifying selection appears to have been the primary driving force of spermatophyte AP2 group evolution, although positive selection occurred in the AP2 clade. The transition from exon to intron of AtAP2 in Arabidopsis mutant leads to the loss of gene function and the same situation was found in AtTOE2. Combining this evolutionary analysis and published research, the results suggest that typical AP2 group genes may first appear in gymnosperms and diverged in angiosperms, following expansion of group members and functional differentiation. In angiosperms, AP2 genes (AP2 clade) inherited key functions from ancestors and other genes of AP2 group lost most function but just remained flowering time controlling in gene formation. In this study, the phylogenies of AP2 group genes in spermatophytes was analyzed, which supported the evidence for the research of gene functional evolution of AP2 group.

Insight into the AP2/ERF transcription factor superfamily in sesame and expression profiling of DREB subfamily under drought stress

BACKGROUND

Sesame is an important oilseed crop mainly grown in inclement areas with high temperatures and frequent drought. Thus, drought constitutes one of the major constraints of its production. The AP2/ERF is a large family of transcription factors known to play significant roles in various plant processes including biotic and abiotic stress responses. Despite their importance, little is known about sesame AP2/ERF genes. This constitutes a limitation for drought-tolerance candidate genes discovery and breeding for tolerance to water deficit.

RESULTS

One hundred thirty-two AP2/ERF genes were identified in the sesame genome. Based on the number of domains, conserved motifs, genes structure and phylogenetic analysis including 5 relatives species, they were classified into 24 AP2, 41 DREB, 61 ERF, 4 RAV and 2 Soloist. The number of sesame AP2/ERF genes was relatively few compared to that of other relatives, probably due to gene loss in ERF and DREB subfamilies during evolutionary process. In general, the AP2/ERF genes were expressed differently in different tissues but exhibited the highest expression levels in the root. Mostly all DREB genes were responsive to drought stress. Regulation by drought is not specific to one DREB group but depends on the genes and the group A6 and A1 appeared to be more actively expressed to cope with drought.

CONCLUSIONS

This study provides insights into the classification, evolution and basic functional analysis of AP2/ERF genes in sesame which revealed their putative involvement in multiple tissue-/developmental stages. Out of 20 genes which were significantly up- /down-regulated under drought stress, the gene AP2si16 may be considered as potential candidate gene for further functional validation as well for utilization in sesame improvement programs for drought stress tolerance.

Genome-wide dissection of AP2/ERF and HSP90 gene families in five legumes and expression profiles in chickpea and pigeonpea

APETALA2/ethylene response factor (AP2/ERF) and heat-shock protein 90 (HSP90) are two significant classes of transcription factor and molecular chaperone proteins which are known to be implicated under abiotic and biotic stresses. Comprehensive survey identified a total of 147 AP2/ERF genes in chickpea, 176 in pigeonpea, 131 in Medicago, 179 in common bean and 140 in Lotus, whereas the number of HSP90 genes ranged from 5 to 7 in five legumes. Sequence alignment and phylogenetic analyses distinguished AP2, ERF, DREB, RAV and soloist proteins, while HSP90 proteins segregated on the basis of their cellular localization. Deeper insights into the gene structure allowed ERF proteins to be classified into AP2s based on DNA-binding domains, intron arrangements and phylogenetic grouping. RNA-seq and quantitative real-time PCR (qRT-PCR) analyses in heat-stressed chickpea as well as Fusarium wilt (FW)- and sterility mosaic disease (SMD)-stressed pigeonpea provided insights into the modus operandi of AP2/ERF and HSP90 genes. This study identified potential candidate genes in response to heat stress in chickpea while for FW and SMD stresses in pigeonpea. For instance, two DREB genes (Ca_02170 and Ca_16631) and three HSP90 genes (Ca_23016, Ca_09743 and Ca_25602) in chickpea can be targeted as potential candidate genes. Similarly, in pigeonpea, a HSP90 gene, C.cajan_27949, was highly responsive to SMD in the resistant genotype ICPL 20096, can be recommended for further functional validation. Also, two DREB genes, C.cajan_41905 and C.cajan_41951, were identified as leads for further investigation in response to FW stress in pigeonpea.

Gene expression analysis and SNP/InDel discovery to investigate yield heterosis of two rubber tree F1 hybrids

As an important industrial material, natural rubber is mainly harvested from the rubber tree. Rubber tree breeding is inefficient, expensive and time-consuming, whereas marker-assisted selection is a feasible method for early selection of high-yield hybrids. We thus sequenced and analyzed the transcriptomes of two parent rubber trees (RRIM 600 and PR 107) and their most productive hybrids (RY 7-33-97 and RY 7-20-59) to understand their gene expression patterns and genetic variations including single nucleotide polymorphisms (SNPs) and small insertions/deletions (InDels). We discovered >31,000 genetic variations in 112,702 assembled unigenes. Our results showed that the higher yield in F₁ hybrids was positively associated with their higher genome heterozygosity, which was further confirmed by genotyping 10 SNPs in 20 other varieties. We also showed that RY 7-33-97 and RY 7-20-59 were genetically closer to RRIM 600 and PR 107, respectively, in agreement with both their phenotypic similarities and gene expression profiles. After identifying ethylene- and jasmonic acid–responsive genes at the transcription level, we compared and analyzed the genetic variations underlying rubber biosynthesis and the jasmonic acid and ethylene pathways in detail. Our results suggest that genome-wide genetic variations play a substantive role in maintaining rubber tree heterosis.

Molecular mechanism of ethylene stimulation of latex yield in rubber tree (Hevea brasiliensis) revealed by de novo sequencing and transcriptome analysis

Background

Rubber tree (Hevea brasiliensis) is an important industrial crop cultivated in tropical areas for natural rubber production. Treatment of the bark of rubber trees with ehephon (an ethylene releaser) has been a routine measure to increase latex yield, but the molecular mechanism behind the stimulation of rubber production by ethylene still remains a puzzle. Deciphering the enigma is of great importance for improvement of rubber tree for high yield.

Results

De novo sequencing and assembly of the bark transciptomes of Hevea brasiliensis induced with ethephon for 8 h (E8) and 24 h (E24) were performed. 51,965,770, 52,303,714 and 53,177,976 high-quality clean reads from E8, E24 and C (control) samples were assembled into 81,335, 80,048 and 80,800 unigenes respectively, with a total of 84,425 unigenes and an average length of 1,101 bp generated. 10,216 and 9,374 differentially expressed genes (DEGs) in E8 and E24 compared with C were respectively detected. The expression of several enzymes in crucial points of regulation in glycolysis were up-regulated and DEGs were not significantly enriched in isopentenyl diphosphate (IPP) biosynthesis pathway. In addition, up-regulated genes of great regulatory importance in carbon fixation (Calvin cycle) were identified.

Conclusions

The rapid acceleration of glycolytic pathway supplying precursors for the biosynthesis of IPP and natural rubber, instead of rubber biosynthesis per se, may be responsible for ethylene stimulation of latex yield in rubber tree. The elevated rate of flux throughout the Calvin cycle may account for some durability of ethylene-induced stimulation. Our finding lays the foundations for molecular diagnostic and genetic engineering for high-yielding improvement of rubber tree.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2587-4) contains supplementary material, which is available to authorized users.

Profiling Ethylene-Responsive Genes Expressed in the Latex of the Mature Virgin Rubber Trees Using cDNA Microarray

Ethylene is commonly used as a latex stimulant of Hevea brasiliensis by application of ethephon (chloro-2-ethylphosphonic acid); however, the molecular mechanism by which ethylene increases latex production is not clear. To better understand the effects of ethylene stimulation on the laticiferous cells of rubber trees, a latex expressed sequence tag (EST)-based complementary DNA microarray containing 2,973 unique genes (probes) was first developed and used to analyze the gene expression changes in the latex of the mature virgin rubber trees after ethephon treatment at three different time-points: 8, 24 and 48 h. Transcript levels of 163 genes were significantly altered with fold-change values ≥ 2 or ≤ –2 (q-value < 0.05) in ethephon-treated rubber trees compared with control trees. Of the 163 genes, 92 were up-regulated and 71 down-regulated. The microarray results were further confirmed using real-time quantitative reverse transcript-PCR for 20 selected genes. The 163 ethylene-responsive genes were involved in several biological processes including organic substance metabolism, cellular metabolism, primary metabolism, biosynthetic process, cellular response to stimulus and stress. The presented data suggest that the laticifer water circulation, production and scavenging of reactive oxygen species, sugar metabolism, and assembly and depolymerization of the latex actin cytoskeleton might play important roles in ethylene-induced increase of latex production. The results may provide useful insights into understanding the molecular mechanism underlying the effect of ethylene on latex metabolism of H. brasiliensis.

Transcriptional identification and characterization of differentially expressed genes associated with embryogenesis in radish (Raphanus sativus L.)

Embryogenesis is an important component in the life cycle of most plant species. Due to the difficulty in embryo isolation, the global gene expression involved in plant embryogenesis, especially the early events following fertilization are largely unknown in radish. In this study, three cDNA libraries from ovules of radish before and after fertilization were sequenced using the Digital Gene Expression (DGE) tag profiling strategy. A total of 5,777 differentially expressed transcripts were detected based on pairwise comparison in the three libraries (0_DAP, 7_DAP and 15_DAP). Results from Gene Ontology (GO) and pathway enrichment analysis revealed that these differentially expressed genes (DEGs) were implicated in numerous life processes including embryo development and phytohormones biosynthesis. Notably, some genes encoding auxin response factor (ARF ), Leafy cotyledon1 (LEC1) and somatic embryogenesis receptor-like kinase (SERK ) known to be involved in radish embryogenesis were differentially expressed. The expression patterns of 30 genes including LEC1-2, AGL9, LRR, PKL and ARF8-1 were validated by qRT-PCR. Furthermore, the cooperation between miRNA and mRNA may play a pivotal role in the radish embryogenesis process. This is the first report on identification of DEGs profiles related to radish embryogenesis and seed development. These results could facilitate further dissection of the molecular mechanisms underlying embryogenesis and seed development in radish.

Origination, Expansion, Evolutionary Trajectory, and Expression Bias of AP2/ERF Superfamily in Brassica napus

The AP2/ERF superfamily, one of the most important transcription factor families, plays crucial roles in response to biotic and abiotic stresses. So far, a comprehensive evolutionary inference of its origination and expansion has not been available. Here, we identified 515 AP2/ERF genes in B. napus, a neo-tetraploid forming ~7500 years ago, and found that 82.14% of them were duplicated in the tetraploidization. A prominent subgenome bias was revealed in gene expression, tissue-specific, and gene conversion. Moreover, a large-scale analysis across plants and alga suggested that this superfamily could have been originated from AP2 family, expanding to form other families (ERF, and RAV). This process was accompanied by duplicating and/or alternative deleting AP2 domain, intragenic domain sequence conversion, and/or by acquiring other domains, resulting in copy number variations, alternatively contributing to functional innovation. We found that significant positive selection occurred at certain critical nodes during the evolution of land plants, possibly responding to changing environment. In conclusion, the present research revealed origination, functional innovation, and evolutionary trajectory of the AP2/ERF superfamily, contributing to understanding their roles in plant stress tolerance.

Genome-wide identification of rubber tree (Hevea brasiliensis Muell. Arg.) aquaporin genes and their response to ethephon stimulation in the laticifer, a rubber-producing tissue

BACKGROUND

Natural rubber, an important industrial raw material, is specifically synthesized in laticifers located inside the rubber tree (Hevea brasiliensis Muell. Arg.) trunk. Due to the absence of plasmodesmata, the laticifer water balance is mediated by aquaporins (AQPs). However, to date, the characterization of H. brasiliensis AQPs (HbAQPs) is still in its infancy.

RESULTS

In this study, 51 full-length AQP genes were identified from the rubber tree genome. The phylogenetic analysis assigned these AQPs to five subfamilies, including 15 plasma membrane intrinsic proteins (PIPs), 17 tonoplast intrinsic proteins (TIPs), 9 NOD26-like intrinsic proteins (NIPs), 4 small basic intrinsic proteins (SIPs) and 6 X intrinsic proteins (XIPs). Functional prediction based on the analysis of the aromatic/arginine (ar/R) selectivity filter, Froger's positions and specificity-determining positions (SDPs) showed a remarkable difference in substrate specificity among subfamilies. Homology analysis supported the expression of 44 HbAQP genes in at least one of the examined tissues. Furthermore, deep sequencing of the laticifer transcriptome in the form of latex revealed a key role of several PIP subfamily members in the laticifer water balance, and qRT-PCR analysis showed diverse expression patterns of laticifer-expressed HbAQP genes upon ethephon treatment, a widely-used practice for the stimulation of latex yield.

CONCLUSIONS

This study provides an important genetic resource of HbAQP genes, which will be useful to improve the water use efficiency and latex yield of Hevea.

Molecular cloning, expression profiles and characterization of a glutathione reductase in Hevea brasiliensis

Glutathione reductase (GR; EC 1.8.1.7) is an important oxidoreductase that can protect organisms against various oxidative stresses. In this study, a new GR gene, named as HbGR2, was isolated from Hevea brasiliensis. The HbGR2 cDNA contained a 1674-bp open reading frame encoding 557 amino acids and the deduced HbGR2 protein showed high identities to the chloroplastic GRs from other plant species. HbGR2 was localized in the chloroplasts of tobacco mesophyll protoplasts. The cis-acting regulatory elements related to stress or hormone responses were predicted in the promoter region of HbGR2. The results from real-time RT-PCR analyses showed that HbGR2 was expressed throughout different tissues and developmental stages of leaves. Besides being related to tapping panel dryness (TPD), HbGR2 was regulated by several treatments including ethephon (ET), methyl jasmonate (MeJA), drought, low temperature, high salt, wounding and hydrogen peroxide (H2O2). The Escherichia coli (E. coli) cells overexpressing HbGR2 markedly increased their tolerance and survival at high concentrations of H2O2, suggesting that HbGR2 might play an important role in oxidative stress response in Hevea brasiliensis.

Functional Characterization of Hevea brasiliensis CRT/DRE Binding Factor 1 Gene Revealed Regulation Potential in the CBF Pathway of Tropical Perennial Tree

Rubber trees (Hevea brasiliensis) are susceptible to low temperature and therefore are only planted in the tropical regions. In the past few decades, although rubber trees have been successfully planted in the northern margin of tropical area in China, they suffered from cold injury during the winter. To understand the physiological response under cold stress, we isolated a C-repeat binding factor 1 (CBF1) gene from the rubber tree. This gene (HbCBF1) was found to respond to cold stress but not drought or ABA stress. The corresponding HbCBF1 protein showed CRT/DRE binding activity in gel shift experiment. To further characterize its molecular function, the HbCBF1 gene was overexpressed in Arabidopsis. The HbCBF1 over expression (OE) line showed enhanced cold resistance and relatively slow dehydration, and the expression of Arabidopsis CBF pathway downstream target genes, e.g. AtCOR15a and AtRD29a, were significantly activated under non-acclimation condition. These data suggest HbCBF1 gene is a functional member of the CBF gene family, and may play important regulation function in rubber tree.

Involvement of Ethylene in the Latex Metabolism and Tapping Panel Dryness of Hevea brasiliensis

Ethephon, an ethylene releaser, is used to stimulate latex production in Hevea brasiliensis. Ethylene induces many functions in latex cells including the production of reactive oxygen species (ROS). The accumulation of ROS is responsible for the coagulation of rubber particles in latex cells, resulting in the partial or complete stoppage of latex flow. This study set out to assess biochemical and histological changes as well as changes in gene expression in latex and phloem tissues from trees grown under various harvesting systems. The Tapping Panel Dryness (TPD) susceptibility of Hevea clones was found to be related to some biochemical parameters, such as low sucrose and high inorganic phosphorus contents. A high tapping frequency and ethephon stimulation induced early TPD occurrence in a high latex metabolism clone and late occurrence in a low latex metabolism clone. TPD-affected trees had smaller number of laticifer vessels compared to healthy trees, suggesting a modification of cambial activity. The differential transcript abundance was observed for twenty-seven candidate genes related to TPD occurrence in latex and phloem tissues for ROS-scavenging, ethylene biosynthesis and signalling genes. The predicted function for some Ethylene Response Factor genes suggested that these candidate genes should play an important role in regulating susceptibility to TPD.

Transcriptome-based discovery of AP2/ERF transcription factors related to temperature stress in tea plant (Camellia sinensis)

Tea plant (Camellia sinensis) is an important natural resource for the global supply of non-alcoholic beverage production. The extension of tea plant cultivation is challenged by biotic and abiotic stresses. Transcription factors (TFs) of the APETALA 2 (AP2)/ethylene-responsive factor (ERF) family are involved in growth and anti-stresses through multifaceted transcriptional regulation in plants. This study comprehensively analyzed AP2/ERF family TFs from C. sinensis on the basis of the transcriptome sequencing data of four tea plant cultivars, namely, 'Yunnanshilixiang', 'Chawansanhao', 'Ruchengmaoyecha', and 'Anjibaicha'. A total of 89 putative AP2/ERF transcription factors with full-length AP2 domain were identified from C. sinensis and classified into five subfamilies, namely, AP2, dehydration-responsive-element-binding (DREB), ERF, related to ABI3/VP (RAV), and Soloist. All identified CsAP2/ERF genes presented relatively stable expression levels in the four tea plant cultivars. Many groups also showed cultivar specificity. Five CsAP2/ERF genes from each AP2/ERF subfamily (DREB, ERF, AP2, and RAV) were related to temperature stresses; these results indicated that AP2/ERF TFs may play important roles in abnormal temperature stress response in C. sinensis.

MicroRNAs and targets in senescent litchi fruit during ambient storage and post-cold storage shelf life

BACKGROUND

Litchi has a high commercial value due to its bright color and rich nutrients. However, it deteriorates with the pericarp turning brown within 1-2 days after harvest. The factors that mediate litchi fruit senescence are complicated. MicroRNAs act as negative regulators involved in almost every physiological process. To understand the mechanism of litchi fruit senescence and pericarp browning at the miRNA level, five small RNA libraries and a degradome library prepared from the pericarp of litchi fruit subjected to ambient storage and post-cold storage shelf life were sequenced.

RESULTS

By aligning the sRNA reads onto the litchi unigene assembly, 296 miRNAs belonging to 49 known miRNA families were first identified from litchi. In addition, 11 litchi-specific miRNAs were identified. Among these, 167 known miRNAs were identified to cleave 197 targets, and three litchi-specific miRNAs were found to have five targets. Through combined analysis of stem-loop quantitative real-time polymerase chain reaction (qRT-PCR) and transcriptome profiling, 14 miRNA-target pairs were found to be actively involved in litchi fruit senescence-related processes, including energy regulation, anthocyanin metabolism, hormone signaling, and pathogen-infection defense.

CONCLUSIONS

A network of miRNA-targets that regulate litchi fruit senescence has been proposed, revealing the miRNA-mediated regulation in senescent litchi fruit. This will aid in developing new strategies to postpone the senescence of litchi fruit and other horticultural products.

Transcriptome sequencing and analysis of rubber tree (Hevea brasiliensis Muell.) to discover putative genes associated with tapping panel dryness (TPD)

Background

Tapping panel dryness (TPD) involves in the partial or complete cessation of latex flow thus seriously affect latex production in rubber tree (Hevea brasiliensis). Numerous studies have been conducted to define its origin and nature, but the molecular nature and mechanism of TPD occurrence remains unknown. This study is committed to de novo sequencing and comparative analysis of the transcriptomes of healthy (H) and TPD-affected (T) rubber trees to identify the genes and pathways related to the TPD.

Results

Total raw reads of 34,632,012 and 35,913,020 bp were obtained from H and T library, respectively using Illumina Hiseq 2000 sequencing technology. De novo assemblies yielded 141,456 and 169,285 contigs, and 96,070 and 112,243 unigenes from H and T library, respectively.

Among 73597 genes, 22577 genes were identified as differential expressed genes between H and T library via comparative transcript profiling. A majority of genes involved in natural rubber biosynthesis and jasmonate synthesis with most potential relevance in TPD occurrence were found to be differentially expressed.

Conclusions

In TPD-affected trees, the expression of most genes related to the latex biosynthesis and jasmonate synthesis was severely inhibited and is probably the direct cause of the TPD. These new de novo transcriptome data sets provide a significant resource for the discovery of genes related to TPD and improve our understanding of the occurrence and maintainace of TPD.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1562-9) contains supplementary material, which is available to authorized users.

Ethylene Response Factors Are Controlled by Multiple Harvesting Stresses in Hevea brasiliensis

Tolerance of recurrent mechanical wounding and exogenous ethylene is a feature of the rubber tree. Latex harvesting involves tapping of the tree bark and ethephon is applied to increase latex flow. Ethylene is an essential element in controlling latex production. The ethylene signalling pathway leads to the activation of Ethylene Response Factor (ERF) transcription factors. This family has been identified in Hevea brasiliensis. This study set out to understand the regulation of ERF genes during latex harvesting in relation to abiotic stress and hormonal treatments. Analyses of the relative transcript abundance were carried out for 35 HbERF genes in latex, in bark from mature trees and in leaves from juvenile plants under multiple abiotic stresses. Twenty-one HbERF genes were regulated by harvesting stress in laticifers, revealing an overrepresentation of genes in group IX. Transcripts of three HbERF-IX genes from HbERF-IXc4, HbERF-IXc5 and HbERF-IXc6 were dramatically accumulated by combining wounding, methyl jasmonate and ethylene treatments. When an ethylene inhibitor was used, the transcript accumulation for these three genes was halted, showing ethylene-dependent induction. Subcellular localization and transactivation experiments confirmed that several members of HbERF-IX are activator-type transcription factors. This study suggested that latex harvesting induces mechanisms developed for the response to abiotic stress. These mechanisms probably depend on various hormonal signalling pathways. Several members of HbERF-IX could be essential integrators of complex hormonal signalling pathways in Hevea.

Isolation and molecular characterization of 1-aminocyclopropane-1-carboxylic acid synthase genes in Hevea brasiliensis

Ethylene is an important factor that stimulates Hevea brasiliensis to produce natural rubber. 1-Aminocyclopropane-1-carboxylic acid synthase (ACS) is a rate-limiting enzyme in ethylene biosynthesis. However, knowledge of the ACS gene family of H. brasiliensis is limited. In this study, nine ACS-like genes were identified in H. brasiliensis. Sequence and phylogenetic analysis results confirmed that seven isozymes (HbACS1–7) of these nine ACS-like genes were similar to ACS isozymes with ACS activity in other plants. Expression analysis results showed that seven ACS genes were differentially expressed in roots, barks, flowers, and leaves of H. brasiliensis. However, no or low ACS gene expression was detected in the latex of H. brasiliensis. Moreover, seven genes were differentially up-regulated by ethylene treatment.These results provided relevant information to help determine the functions of the ACS gene in H. brasiliensis, particularly the functions in regulating ethylene stimulation of latex production.