BackgroundIn plants, ERF genes participate in a variety of regulatory pathways, such as plant growth and biotic and/or abiotic stress responses. Although the genome of Chinese white pear (‘Dangshansuli’) has been released, knowledge regarding the ERF family in pear, such as gene functions, evolutionary history and expression patterns, remains limited.ResultsIn our study, a total of 155 members of ERF families were identified in pear (Pyrus bretschneideri). The Ka and Ks values suggested that whole-genome duplication (WGD) and dispersed duplication have effectively contributed to the expansion of the pear ERF family. Gene structure and phylogeny analysis divided the PbrERF family into 12 groups, and their gene functions were predicted by comparative analysis. qRT-PCR was carried out to verify the relative expression levels of 7 genes in group III using wild and cultivated pear fruits at three key developmental stages. Wild samples had higher expression of these genes than cultivated samples, especially at the enlarged fruit stage. The transcriptome data of pear seedlings subjected to dehydration treatment further revealed that 4 of the 7 genes responded to drought conditions.ConclusionThe AP2/ERF gene family is greatly expanded in pear. Comparative analysis revealed the probability of ERF genes performing functional roles in multiple pathways. Expression analysis at different stages of pear fruit development in wild and cultivated samples indicated that genes in group III might be involved in abiotic and/or biotic stresses. Further transcriptome data on seedlings subjected to drought treatment verified the potential role of ERF genes in stress response. These results will provide a valuable reference for understanding the function and evolution of the ERF family in higher plants.Electronic supplementary materialThe online version of this article (10.1186/s12870-018-1265-x) contains supplementary material, which is available to authorized users.

Bread wheat is an essential crop with the second-highest global production after maize. Currently, wheat diseases are a serious threat to wheat production. Therefore, efficient breeding for disease resistance is extremely urgent in modern wheat. Here, we identified 2012 NLR genes from hexaploid wheat, and Ks values of paired syntenic NLRs showed a significant peak at 3.1–6.3 MYA, which exactly coincided with the first hybridization event between A and B genome lineages at ~5.5 MYA. We provided a landscape of dynamic diversity of NLRs from Triticum and Aegilops and found that NLR genes have higher diversity in wild progenitors and relatives. Further, most NLRs had opposite diversity patterns between genic and 2 Kb-promoter regions, which might respectively link sub/neofunctionalization and loss of duplicated NLR genes. Additionally, we identified an alien introgression of chromosome 4A in tetraploid emmer wheat, which was similar to that in hexaploid wheat. Transcriptome data from four experiments of wheat disease resistance helped to profile the expression pattern of NLR genes and identified promising NLRs involved in broad-spectrum disease resistance. Our study provided insights into the diversity evolution of NLR genes and identified beneficial NLRs to deploy into modern wheat in future wheat disease-resistance breeding.

The MYB transcription factor superfamily includes key regulators of plant development and responses to environmental changes. The diversity of lifestyles and morphological characteristics exhibited by plants are potentially associated with the genomic dynamics of the MYB superfamily. With the release of the plant genomes, a comprehensive phylogenomic analysis of the MYB superfamily across Viridiplantae is allowed. The present study performed phylogenetic, phylogenomic, syntenic, horizontal gene transfer, and neo/sub-functionalization analysis of the MYB superfamily to explore the evolutionary contributions of MYB members to species diversification, trait formation, and environmental adaptation in 437 different plant species. We identified major changes in copy number variation and genomic context within subclades across lineages. Multiple MYB subclades showed highly conserved copy number patterns and synteny across flowering plants, whereas others were more dynamic and showed lineage-specific patterns. As examples of lineage-specific morphological divergence, we hypothesize that the gain of a MYB orthogroup associated with flower development and environmental responses and an orthogroup associated with auxin and wax biosynthesis in angiosperms were correlated with the emergence of flowering plants, unbiased neo-/sub-functionalization of gene duplicates contributed to environmental adaptation, and species-specific neo-/sub-functionalization contributed to phenotype divergence between species. Transposable element insertion in promoter regions may have facilitated the sub-/neo-functionalization of MYB genes and likely played a tissue-specific role contributing to sub-/neo-functionalization in plant root tissues. This study provides new insights into the evolutionary divergence of the MYB superfamily across major flowering and non-flowering lineages and emphasizes the need for lineage-/tissue-specific characterization to further understand trait variability and environmental adaptation.