Background The sharp increase of plant genome and transcriptome data provide valuable resources to investigate evolutionary consequences of gene duplication in a range of taxa, and unravel common principles underlying duplicate gene retention. Results We survey 141 sequenced plant genomes to elucidate consequences of gene and genome duplication, processes central to the evolution of biodiversity. We develop a pipeline named DupGen_finder to identify different modes of gene duplication in plants. Genes derived from whole-genome, tandem, proximal, transposed, or dispersed duplication differ in abundance, selection pressure, expression divergence, and gene conversion rate among genomes. The number of WGD-derived duplicate genes decreases exponentially with increasing age of duplication events—transposed duplication- and dispersed duplication-derived genes declined in parallel. In contrast, the frequency of tandem and proximal duplications showed no significant decrease over time, providing a continuous supply of variants available for adaptation to continuously changing environments. Moreover, tandem and proximal duplicates experienced stronger selective pressure than genes formed by other modes and evolved toward biased functional roles involved in plant self-defense. The rate of gene conversion among WGD-derived gene pairs declined over time, peaking shortly after polyploidization. To provide a platform for accessing duplicated gene pairs in different plants, we constructed the Plant Duplicate Gene Database. Conclusions We identify a comprehensive landscape of different modes of gene duplication across the plant kingdom by comparing 141 genomes, which provides a solid foundation for further investigation of the dynamic evolution of duplicate genes. Electronic supplementary material The online version of this article (10.1186/s13059-019-1650-2) contains supplementary material, which is available to authorized users.

BackgroundHeat shock transcription factors (Hsfs), which act as important transcriptional regulatory proteins in eukaryotes, play a central role in controlling the expression of heat-responsive genes. At present, the genomes of Chinese white pear (‘Dangshansuli’) and five other Rosaceae fruit crops have been fully sequenced. However, information about the Hsfs gene family in these Rosaceae species is limited, and the evolutionary history of the Hsfs gene family also remains unresolved.ResultsIn this study, 137 Hsf genes were identified from six Rosaceae species (Pyrus bretschneideri, Malus × domestica, Prunus persica, Fragaria vesca, Prunus mume, and Pyrus communis), 29 of which came from Chinese white pear, designated as PbHsf. Based on the structural characteristics and phylogenetic analysis of these sequences, the Hsf family genes could be classified into three main groups (classes A, B, and C). Segmental and dispersed duplications were the primary forces underlying Hsf gene family expansion in the Rosaceae. Most of the PbHsf duplicated gene pairs were dated back to the recent whole-genome duplication (WGD, 30–45 million years ago (MYA)). Purifying selection also played a critical role in the evolution of Hsf genes. Transcriptome data demonstrated that the expression levels of the PbHsf genes were widely different. Six PbHsf genes were upregulated in fruit under naturally increased temperature.ConclusionA comprehensive analysis of Hsf genes was performed in six Rosaceae species, and 137 full length Hsf genes were identified. The results presented here will undoubtedly be useful for better understanding the complexity of the Hsf gene family and will facilitate functional characterization in future studies.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-014-0401-5) contains supplementary material, which is available to authorized users.

The MYB superfamily is large and functionally diverse in plants. To date, MYB family genes have not yet been identified in Chinese white pear (Pyrus bretschneideri), and their functions remain unclear. In this study, we identified 231 genes as candidate MYB genes and divided them into four subfamilies. The R2R3-MYB (PbrMYB) family shared an R2R3 domain with 104 amino acid residues, including five conserved tryptophan residues. The Pbr MYB family was divided into 37 functional subgroups including 33 subgroups which contained both MYB genes of Rosaceae plants and AtMYB genes, and four subgroups which included only Rosaceae MYB genes or AtMYB genes. PbrMYB genes with similar functions clustered into the same subgroup, indicating functional conservation. We also found that whole-genome duplication (WGD) and dispersed duplications played critical roles in the expansion of the MYB family. The 87 Pbr MYB duplicated gene pairs dated back to the two WGD events. Purifying selection was the primary force driving Pbr MYB gene evolution. The 15 gene pairs presented 1-7 codon sites under positive selection. A total of 147 expressed genes were identified from RNA-sequencing data of fruit, and six Pbr MYB members in subgroup C1 were identified as important candidate genes in the regulation of lignin synthesis by quantitative real-time PCR analysis. Further correlation analysis revealed that six PbrMYBs were significantly correlated with five structural gene families (F5H, HCT, CCR, POD and C3'H) in the lignin pathway. The phylogenetic, evolution and expression analyses of the MYB gene family in Chinese white pear establish a solid foundation for future comprehensive functional analysis of Pbr MYB genes.

Pear is an important fruit crop of the Rosaceae family and has experienced two rounds of ancient whole-genome duplications (WGDs). However, whether different types of gene duplications evolved differently after duplication remains unclear in the pear genome. In this study, we identified the different modes of gene duplication in pear. Duplicate genes derived from WGD, tandem, proximal, retrotransposed, DNA-based transposed or dispersed duplications differ in genomic distribution, gene features, selection pressure, expression divergence, regulatory divergence and biological roles. Widespread sequence, expression and regulatory divergence have occurred between duplicate genes over the 30–45 million years of evolution after the recent genome duplication in pear. The retrotransposed genes show relatively higher expression and regulatory divergence than other gene duplication modes. In contrast, WGD genes underwent a slower sequence divergence and may be influenced by abundant gene conversion events. Moreover, the different classes of duplicate genes exhibited biased functional roles. We also investigated the evolution and expansion patterns of the gene families involved in sugar and organic acid metabolism pathways, which are closely related to the fruit quality and taste in pear. Single-gene duplications largely account for the extensive expansion of gene families involved in the sorbitol metabolism pathway in pear. Gene family expansion was also detected in the sucrose metabolism pathway and tricarboxylic acid cycle pathways. Thus, this study provides insights into the evolutionary fates of duplicated genes.

BackgroundMicroRNAs (miRNAs) are a class of small, endogenous RNAs that take part in regulating genes through mediating gene expressions at the post-transcriptional level in plants. Previous studies have reported miRNA identification in various plants ranging from model plants to perennial fruit trees. However, the role of miRNAs in pear (Pyrus bretschneideri) fruit development is not clear. Here, we investigated the miRNA profiles of pear fruits from different time stages during development with Illumina HiSeq 2000 platform and bioinformatics analysis. Quantitative real-time PCR was used to validate the expression levels of miRNAs.ResultsBoth conserved and species-specific miRNAs in pear have been identified in this study. Total reads, ranging from 19,030,925 to 25,576,773, were obtained from six small RNA libraries constructed for different stages of fruit development after flowering. Comparative profiling showed that an average of 90 miRNAs was expressed with significant differences between various developmental stages. KEGG pathway analysis on 2,216 target genes of 188 known miRNAs and 1,127 target genes of 184 novel miRNAs showed that miRNAs are widely involved in the regulation of fruit development. Among these, a total of eleven miRNAs putatively participate in the pathway of lignin biosynthesis, nine miRNAs were identified to take part in sugar and acid metabolism, and MiR160 was identified to regulate auxin response factor.ConclusionComparative analysis of miRNAomes during pear fruit development is presented, and miRNAs were proved to be widely involved in the regulation of fruit development and formation of fruit quality, for example through lignin synthesis, sugar and acid metabolism, and hormone signaling. Combined with computational analysis and experimental confirmation, the research contributes valuable information for further functional research of microRNA in fruit development for pear and other species.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-953) contains supplementary material, which is available to authorized users.