Search

A comparative investigation was conducted to evaluate transcriptional changes in guard cells (GCs) of closely related halophytic (Chenopodium quinoa) and glycophytic (Spinacia oleracea) species. Plants were exposed to 3 weeks of 250 mM sodium chloride treatment, and GC-enriched epidermal fragments were mechanically prepared. In both species, salt-responsive genes were mainly related to categories of protein metabolism, secondary metabolites, signal transduction and transport systems. Genes related to abscisic acid (ABA) signaling and ABA biosynthesis were strongly induced in quinoa but not in spinach GCs. Also, expression of the genes encoding transporters of amino acids, proline, sugars, sucrose and potassium increased in quinoa GCs under salinity stress. Analysis of cell-wall-related genes suggests that genes involved in lignin synthesis (e.g. lignin biosynthesis LACCASE 4) were highly upregulated by salt in spinach GCs. In contrast, transcripts related to cell wall plasticity Pectin methylesterase3 (PME3) were highly induced in quinoa. Faster stomatal response to light and dark measured by observing kinetics of changes in stomatal conductance in quinoa might be associated with higher plasticity of the cell wall regulated by PME3 Furthermore, genes involved in the inhibition of stomatal development and differentiation were highly expressed by salt in quinoa, but not in spinach. These changes correlated with reduced stomatal density and index in quinoa, thus improving its water use efficiency. The fine modulation of transporters, cell wall modification and controlling stomatal development in GCs of quinoa may have resulted in high K+/Na+ ratio, lower stomatal conductance and higher stomatal speed for better adaptation to salinity stress in quinoa.

Broomcorn millet (Panicum miliaceum L.) was domesticated in northern China at least 7,000 years ago and was subsequentially adopted as a cereal in many areas throughout Eurasia. One such locale is Areni-1 an archaeological cave site in Southern Armenia, a region that has an important history in crop domestication. The rich botanical material found at Areni-1 includes grains identified by morphology as broomcorn millet that were 14C dated to the medieval era (873 ± 36 CE and 1118 ± 35 CE). To retrace the demographic history of these broomcorn millet samples, we used ancient DNA extraction and hybridization capture enrichment to sequence and assemble three chloroplast genomes from the Areni-1 grains and then compared these sequences to 50 modern chloroplast genomes. Overall, the chloroplast genomes contained a low amount of diversity and little inference on broomcorn demography could be made. However, in a phylogeny the chloroplast genomes separated into two clades with strong bootstrap support, similar to what has been reported for nuclear DNA from broomcorn millet. In a haplotype network, the chloroplast genomes of two accessions of wild (undomesticated) broomcorn millet contained a relatively large number of variants, 11 SNPs. These SNPs were not present in the domesticated varieties, suggesting these wild accessions may not be directly related to the lineages that underwent domestication or that broomcorn millet may have undergone a domestication bottleneck resulting in lost diversity in the chloroplast genome. These results demonstrate that broomcorn millet from archaeological sites can preserve DNA for at least 1000 years and serve as a genetic resource to study the domestication of this cereal crop.

Genome fractionation (also known as diploidization) frequently occurs following paleopolyploidization events. Biased fractionation between subgenomes has been found in some paleo-allopolyploids, while this phenomenon is absent in paleo-autopolyploids. Pear (Pyrus bretschneideri Rehd.) experienced a recent whole-genome duplication (WGD, ~30 million years ago); however, the evolutionary fate of the two subgenomes derived from this WGD event is not clear. In this study, we identified the two paleo-subgenomes in pear using peach (Prunus persica) as an outgroup and investigated differences in the gene loss rate, evolutionary rate, gene expression level, and DNA methylation level between these two subgenomes. Fractionation bias was not found between the two pear subgenomes, which evolved at similar evolutionary rates. The DNA methylation level of the two subgenomes showed little bias, and we found no expression dominance between the subgenomes. However, we found that singleton genes and homeologous genes within each subgenome showed divergent evolutionary patterns of selective constraints, expression and epigenetic modification. These results provide insights into subgenome evolution following paleopolyploidization in pear.

Broomcorn millet (Panicum miliaceum L.) is the most water-efficient cereal and one of the earliest domesticated plants. Here we report its high-quality, chromosome-scale genome assembly using a combination of short-read sequencing, single-molecule real-time sequencing, Hi-C, and a high-density genetic map. Phylogenetic analyses reveal two sets of homologous chromosomes that may have merged ~5.6 million years ago, both of which exhibit strong synteny with other grass species. Broomcorn millet contains 55,930 protein-coding genes and 339 microRNA genes. We find Paniceae-specific expansion in several subfamilies of the BTB (broad complex/tramtrack/bric-a-brac) subunit of ubiquitin E3 ligases, suggesting enhanced regulation of protein dynamics may have contributed to the evolution of broomcorn millet. In addition, we identify the coexistence of all three C4 subtypes of carbon fixation candidate genes. The genome sequence is a valuable resource for breeders and will provide the foundation for studying the exceptional stress tolerance as well as C4 biology., Broomcorn millet is one of the earliest domesticated plants and has the highest water use efficiency among cereals. Here, the authors report its genome assembly and annotation, which provides a valuable resource for breeders and paves the way for studying plant drought tolerance and C4 photosynthesis.

Aiming at the need of video abnormal events to be located in pixel-level regions, a video abnormal event detection method based on CNN (Convolutional Neural Networks) and multiple instance learning is proposed. Firstly, the Gaussian background model is used to extract the moving targets in the video, and the connected regions of the moving targets are obtained by the image processing method. Secondly, the pre-trained VGG16 model is used to extract the features of the connected regions what construct multiple instance learning packages. Finally, the multiple instance learning model is trained using MISVM (Multiple-Instance Support Vector Machines) and NSK (Normalized Set Kernel) algorithms and predicted at the pixel-level. The experimental results show that the video anomaly detection method based on CNN and multiple instance learning can accurately locate the abnormal events in the pixel-level region.

Epidermal fragments enriched in guard cells (GCs) were isolated from the halophyte quinoa (Chenopodium quinoa Wild.) species, and the response at the proteome level was studied after salinity treatment of 300 mM NaCl for 3 weeks. In total, 2147 proteins were identified, of which 36% were differentially expressed in response to salinity stress in GCs. Up and downregulated proteins included signaling molecules, enzyme modulators, transcription factors and oxidoreductases. The most abundant proteins induced by salt treatment were desiccation-responsive protein 29B (50-fold), osmotin-like protein OSML13 (13-fold), polycystin-1, lipoxygenase, alpha-toxin, and triacylglycerol lipase (PLAT) domain-containing protein 3-like (eight-fold), and dehydrin early responsive to dehydration (ERD14) (eight-fold). Ten proteins related to the gene ontology term &ldquo, response to ABA&rdquo, were upregulated in quinoa GC, this included aspartic protease, phospholipase D and plastid-lipid-associated protein. Additionally, seven proteins in the sucrose&ndash, starch pathway were upregulated in the GC in response to salinity stress, and accumulation of tryptophan synthase and L-methionine synthase (enzymes involved in the amino acid biosynthesis) was observed. Exogenous application of sucrose and tryptophan, L-methionine resulted in reduction in stomatal aperture and conductance, which could be advantageous for plants under salt stress. Eight aspartic proteinase proteins were highly upregulated in GCs of quinoa, and exogenous application of pepstatin A (an inhibitor of aspartic proteinase) was accompanied by higher oxidative stress and extremely low stomatal aperture and conductance, suggesting a possible role of aspartic proteinase in mitigating oxidative stress induced by saline conditions.

Nowadays, people pay more and more attention to security issues. If we can warn the occurrence of abnormal events, we can greatly reduce the occurrence of public security issues. Therefore, the detection of abnormal events in videos is an important research topic. In this paper, an improved deep learning model TCNN-LSTM is obtained by combining CNN with LSTM. We carried out experiments on the standard dataset UMN and UCSD, and achieved good results, with the accuracy reaching 99.54% and 92.36% respectively, and F1 reaching 98.41% and 87.65% respectively, it can be seen that this model has strong generalization ability. At the same time, other models were compared on the UCSD dataset. The AUC of this model reached 97.00%, far higher than that of other models. Experimental results showed that the performance of our model was better than that of other models.

Soil salinity is a major environmental constraint affecting crop growth and threatening global food security. Plants adapt to salinity by optimizing the performance of stomata. Stomata are formed by two guard cells (GCs) that are morphologically and functionally distinct from the other leaf cells. These microscopic sphincters inserted into the wax-covered epidermis of the shoot balance CO2 intake for photosynthetic carbon gain and concomitant water loss. In order to better understand the molecular mechanisms underlying stomatal function under saline conditions, we used proteomics approach to study isolated GCs from the salt-tolerant sugar beet species. Of the 2088 proteins identified in sugar beet GCs, 82 were differentially regulated by salt treatment. According to bioinformatics analysis (GO enrichment analysis and protein classification), these proteins were involved in lipid metabolism, cell wall modification, ATP biosynthesis, and signaling. Among the significant differentially abundant proteins, several proteins classified as “stress proteins” were upregulated, including non-specific lipid transfer protein, chaperone proteins, heat shock proteins, inorganic pyrophosphatase 2, responsible for energized vacuole membrane for ion transportation. Moreover, several antioxidant enzymes (peroxide, superoxidase dismutase) were highly upregulated. Furthermore, cell wall proteins detected in GCs provided some evidence that GC walls were more flexible in response to salt stress. Proteins such as L-ascorbate oxidase that were constitutively high under both control and high salinity conditions may contribute to the ability of sugar beet GCs to adapt to salinity by mitigating salinity-induced oxidative stress.

Soil salinity is a major environmental constraint affecting crop growth and threatening global food security. Plants adapt to salinity by optimizing the performance of stomata. Stomata are formed by two guard cells (GCs) that are morphologically and functionally distinct from the other leaf cells. These microscopic sphincters inserted into the wax-covered epidermis of the shoot balance CO2 intake for photosynthetic carbon gain and concomitant water loss. In order to better understand the molecular mechanisms underlying stomatal function under saline conditions, we used proteomics approach to study isolated GCs from the salt-tolerant sugar beet species. Of the 2088 proteins identified in sugar beet GCs, 82 were differentially regulated by salt treatment. According to bioinformatics analysis (GO enrichment analysis and protein classification), these proteins were involved in lipid metabolism, cell wall modification, ATP biosynthesis, and signaling. Among the significant differentially abundant proteins, several proteins classified as &ldquo, stress proteins&rdquo, were upregulated, including non-specific lipid transfer protein, chaperone proteins, heat shock proteins, inorganic pyrophosphatase 2, responsible for energized vacuole membrane for ion transportation. Moreover, several antioxidant enzymes (peroxide, superoxidase dismutase) were highly upregulated. Furthermore, cell wall proteins detected in GCs provided some evidence that GC walls were more flexible in response to salt stress. Proteins such as L-ascorbate oxidase that were constitutively high under both control and high salinity conditions may contribute to the ability of sugar beet GCs to adapt to salinity by mitigating salinity-induced oxidative stress.

As sessile organisms, plants are highly sensitive to environmental stresses. In response to stresses, globally repressed translation initiation leads to stress granule (SG) formation. Protein liquid-liquid phase separation (LLPS) contributes to SG formation, but a direct link between protein LLPS and stress resistance has not yet been found in plants. Here, we report that two RNA-binding proteins, RBGD2 and RBGD4, function redundantly to improve heat resistance in Arabidopsis. RBGD2 and RBGD4 undergo LLPS in vitro and condense into heat-induced SGs in vivo via tyrosine residue array (TRA). Importantly, disrupting LLPS by mutating TRA abolishes RBGD2/4 condensation in SGs and impairs their protective function against heat stress (HS). Further study found that upon HS, the RBGD2/4 interaction network expands with additional SG proteins and heat-responsive mRNA. Our work shows a mechanistic basis that underlies protein LLPS in HS response in plants and suggests manipulation of protein LLPS as a general strategy to improve plant stress resistance.

Summary The complexity of polyploid Saccharum genomes hindered progress of genome research and crop improvement in sugarcane. To understand their genome structure, transcriptomes of 59 F1 individuals derived from S. officinarum LA Purple and S. robustum Molokai 5829 (2n = 80, x = 10 for both) were sequenced, yielding 11 157 and 8998 SNPs and 83 and 105 linkage groups, respectively. Most markers in each linkage group aligned to single sorghum chromosome. However, 71 interchromosomal rearrangements were detected between sorghum and S. officinarum or S. robustum, and 24 (33.8%) of them were shared between S. officinarum and S. robustum, indicating their occurrence before the speciation event that separated these two species. More than 2000 gene pairs from S. spontaneum, S. officinarum and S. robustum were analysed to estimate their divergence time. Saccharum officinarum and S. robustum diverged about 385 thousand years ago, and the whole‐genome duplication events occurred after the speciation event because of shared interchromosomal rearrangements. The ancestor of these two species diverged from S. spontaneum about 769 thousand years ago, and the reduction in basic chromosome number from 10 to 8 in S. spontaneum occurred after the speciation event but before the two rounds of whole‐genome duplication. Our results proved that S. officinarum is a legitimate species in its own right and not a selection from S. robustum during the domestication process in the past 10 000 years. Our findings rejected a long‐standing hypothesis and clarified the timing of speciation and whole‐genome duplication events in Saccharum.

Background Leucine-rich repeat receptor-like protein kinase (LRR-RLK) is the largest gene family of receptor-like protein kinases (RLKs) and actively participates in regulating the growth, development, signal transduction, immunity, and stress responses of plants. However, the patterns of LRR-RLK gene family evolution in the five main Rosaceae species for which genome sequences are available have not yet been reported. In this study, we performed a comprehensive analysis of LRR-RLK genes for five Rosaceae species: Fragaria vesca (strawberry), Malus domestica (apple), Pyrus bretschneideri (Chinese white pear), Prunus mume (mei), and Prunus persica (peach), which contained 201, 244, 427, 267, and 258 LRR-RLK genes, respectively. Results All LRR-RLK genes were further grouped into 23 subfamilies based on the hidden Markov models approach. RLK-Pelle_LRR-XII-1, RLK-Pelle_LRR-XI-1, and RLK-Pelle_LRR-III were the three largest subfamilies. Synteny analysis indicated that there were 236 tandem duplicated genes in the five Rosaceae species, among which subfamilies XII-1 (82 genes) and XI-1 (80 genes) comprised 68.6%. Conclusions Our results indicate that tandem duplication made a large contribution to the expansion of the subfamilies. The gene expression, tissue-specific expression, and subcellular localization data revealed that LRR-RLK genes were differentially expressed in various organs and tissues, and the largest subfamily XI-1 was highly expressed in all five Rosaceae species, suggesting that LRR-RLKs play important roles in each stage of plant growth and development. Taken together, our results provide an overview of the LRR-RLK family in Rosaceae genomes and the basis for further functional studies. Electronic supplementary material The online version of this article (10.1186/s12864-017-4155-y) contains supplementary material, which is available to authorized users.

Hexokinase genes (HXK) encode enzymes that play an important role in fruit quality and plant development because they participate in the first step of the glycolytic pathway for sugar accumulation and metabolism. In this study, the authors conducted analysis of the phylogeny, gene structure, and expression of HXK genes in pear (Pyrus bretschneideri Rehd.), comparing them with 7 other plant species. A total of 10 HXK protein sequences from pear and 33 HXKs from 7 other species were identified and clustered into 4 groups, matching their motif classification. Duplication events and motif structure variation attributed the diversification of HXK genes, supported by differential expression pattern of PbrHXKs from pear fruit transcriptome data and qRT – PCR verification. Combining pear fruit soluble sugar content data and analyses above, PbrHXK1and PbrHXK3 were identified as important candidate genes in sorbitol, fructose, and glucose accumulation during pear fruit development. PbrHXK1and PbrHXK3 were i...

Supplementary Table S1 and Figure S1-21. (DOCX 5373 kb)

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.

Genome fractionation (also known as diploidization) frequently occurs following paleopolyploidization events. Biased fractionation between subgenomes has been found in some paleo-allopolyploids, while this phenomenon is absent in paleo-autopolyploids. Pear (Pyrus bretschneideri Rehd.) experienced a recent whole-genome duplication (WGD, ~30 million years ago); however, the evolutionary fate of the two subgenomes derived from this WGD event is not clear. In this study, we identified the two paleo-subgenomes in pear using peach (Prunus persica) as an outgroup and investigated differences in the gene loss rate, evolutionary rate, gene expression level, and DNA methylation level between these two subgenomes. Fractionation bias was not found between the two pear subgenomes, which evolved at similar evolutionary rates. The DNA methylation level of the two subgenomes showed little bias, and we found no expression dominance between the subgenomes. However, we found that singleton genes and homeologous genes within each subgenome showed divergent evolutionary patterns of selective constraints, expression and epigenetic modification. These results provide insights into subgenome evolution following paleopolyploidization in pear., Following a trail of duplicated pear DNA Genomic sleuthing reveals insights into an ancestral chromosomal duplication event from the early evolution of the pear tree. Such ‘whole-genome duplication’ is common among plants, and can arise from the combination of genetic material of plants from either the same (autopolyploid) or different (allopolyploid) species. The pear is known to have undergone a duplication event 30 million years ago, and Shaoling Zhang and colleagues at Nanjing Agricultural University in China set out to assess the nature of this genome expansion. After analysis of pear genome data and comparison against other fruit plant species, the researchers determined that the pear is an autopolyploid, and that genes found in just one ‘subgenome’ have subsequently evolved expression profiles that differ noticeably from those genes that are present as multiple homologues amongst the duplicated chromosomes.

Background Pear (Pyrus) is a globally grown fruit, with thousands of cultivars in five domesticated species and dozens of wild species. However, little is known about the evolutionary history of these pear species and what has contributed to the distinct phenotypic traits between Asian pears and European pears. Results We report the genome resequencing of 113 pear accessions from worldwide collections, representing both cultivated and wild pear species. Based on 18,302,883 identified SNPs, we conduct phylogenetics, population structure, gene flow, and selective sweep analyses. Furthermore, we propose a model for the divergence, dissemination, and independent domestication of Asian and European pears in which pear, after originating in southwest China and then being disseminated throughout central Asia, has eventually spread to western Asia, and then on to Europe. We find evidence for rapid evolution and balancing selection for S-RNase genes that have contributed to the maintenance of self-incompatibility, thus promoting outcrossing and accounting for pear genome diversity across the Eurasian continent. In addition, separate selective sweep signatures between Asian pears and European pears, combined with co-localized QTLs and differentially expressed genes, underline distinct phenotypic fruit traits, including flesh texture, sugar, acidity, aroma, and stone cells. Conclusions This study provides further clarification of the evolutionary history of pear along with independent domestication of Asian and European pears. Furthermore, it provides substantive and valuable genomic resources that will significantly advance pear improvement and molecular breeding efforts. Electronic supplementary material The online version of this article (10.1186/s13059-018-1452-y) contains supplementary material, which is available to authorized users.

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.

Quantitative real-time PCR (qRT-PCR) is the preferred method for accurate analysis of targeted gene expression. The aim of this study was to validate optimal reference genes for qRT-PCR analyses in pear pollen. Ten candidate reference genes were chosen and examined in 36 pear samples. Among these, PbrEF1α, PbrGAPDH and PbrAPT were developed to detect the purity of cDNA. Three qRT-PCR analysis programs, geNorm, NormFinder and BestKeeper, were applied for systematic determinations. The results showed that all the candidate reference genes could be used for qRT-PCR normalization in pear pollen. PbrEF1α was the most stable reference gene across all abiotic stresses of pollen samples, especially for high temperature, N-butyl alcohol and NaCl treatments, and also during pollen development. PbrTUB was the optimal reference gene for pear tissue expression. PbrGAPDH was suitable for low temperature treatment and PbrPP2A was suitable for abscisic acid and also for NaCl treatment. For CuCl2 treatment, either of reference genes PbrCYP and PbrUBI was suitable. To illustrate the suitability of these reference genes, the expression profile of PbrVHA-a was used to analyze the expression patterns in different pear tissues, pollen developing stages and NaCl treatment. The results provide a valuable resource for gene expression characterization in pear pollen.

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.

Background Cysteine-rich peptides (CRPs) are gaining recognition as regulators of cell–cell communication in plants. Results We identified 9556 CRPs in 12 plant species and analysed their evolutionary patterns. In most angiosperm plants, whole genome duplication and segmental duplication are the major factors driving the expansion of CRP family member genes, especially signal peptides. About 30% of the CRP genes were found clustered on the chromosomes, except in maize (Zea mays). Considerable collinearities between CRP genes between or within species reveal several syntenic regions on the chromosomes. Different subfamilies display diverse evolutionary rates, suggesting that these subfamilies are subjected to different selective pressures. CRPs in different duplication models also show contrasting evolutionary rates, although the underlying mechanism is unclear because of the complexity of gene evolution. The 1281 positively selected genes identified are probably generated within a certain period of time. While most of these belonged to maize and sorghum (Sorghum bicolor), new CRP functions would also be expected. Up-regulation of 10 CRPs was observed in self-pollinated pear pistils and pollen tubes under self S-RNase treatments in vitro. The expression divergence between different CRP gene duplication types suggests that different duplication mechanisms affected the fate of the duplicated CRPs. Conclusion Our analyses of the evolution of the CRP gene family provides a unique view of the evolution of this large gene family. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3948-3) contains supplementary material, which is available to authorized users.

Red fruits are popular and widely accepted by consumers because of an enhanced appearance and enriched anthocyanins. The molecular mechanism of anthocyanin regulation in red-skinned pear (Pyrus) has been studied, and the genes encoding the biosynthetic steps and several transcription factors (TFs) have been characterized. In this study, a candidate R2R3 MYB TF, PyMYB114, was identified by linkage to the quantitative trait loci (QTL) for red skin color on linkage group 5 in a population of Chinese pear (Pyrus bretschneideri). The function of PyMYB114 was verified by transient transformation in tobacco (Nicotinana tabacum) leaves and strawberry (Fragaria) and pear fruits, resulting in the biosynthesis of anthocyanin. Suppression of PyMYB114 could inhibit anthocyanin biosynthesis in red-skinned pears. The ERF/AP2 TF PyERF3 was found to interact with PyMYB114 and its partner PybHLH3 to co-regulate anthocyanin biosynthesis, as shown by a dual luciferase reporter system and a yeast two-hybrid assay. In addition, the transcript abundance of PyMYB114 and PyMYB10 were correlated, and co-transformation of these two genes into tobacco and strawberry led to enhanced anthocyanin biosynthesis. This interaction network provides insight into the coloration of fruits and the interaction of different TFs to regulate anthocyanin biosynthesis.

Comparison of Ks values for CRP genes. a) Comparison of Ks values for CRP genes from 11 angiosperm plants; b) Comparison of Ks values between positively and negatively selected genes. (PDF 191 kb)

Expression pattern of CRP genes during fruit development. S1â S6 indicate stages of fruit development at 15, 36, 80, 110, 145 and 167Â days after flowering. (PDF 586 kb)

Relationship between WGD events and the number of duplicated CRP genes. a) Brief phylogenetic tree of sequenced species; b) Duplication type of CRP genes from 12 species. WGD data were obtained from the Plant Genome Duplication Database [53]. (PDF 285 kb)

Conserved cysteine motifs in three large subfamilies: Lipid Transfer Protein, Pollen Ole I and DEFL. (PDF 20131 kb)

Comparison of evolutionary rates between four duplication models in strawberry, mei and peach. (PDF 185 kb)

Heat map of the expression patterns of 201 strawberry LRR-RLK genes in different tissues. Red and green colours correspond to up-regulation and down-regulation, respectively. Normalized gene expression values are provided in Additional file 1: Table S10. Figure S2. Heat map of the expression patterns of 244 apple LRR-RLK genes in different tissues. Red and green colours correspond to up-regulation and down-regulation, respectively. Normalized gene expression values are provided in Additional file 1: Table S11. Figure S3. Heat map of the expression patterns of 427 Chinese white pear LRR-RLK genes in different tissues. Red and green colours correspond to up-regulation and down-regulation, respectively. Normalized gene expression values are provided in Additional file 1: Table S12. Figure S4. Heat map of the expression patterns of 267 mei LRR-RLK genes in different tissues. Red and green colours correspond to up-regulation and down-regulation, respectively. Normalized gene expression values are provided in Additional file 1: Table S13. Figure S5. Heat map of the expression patterns of 258 peach LRR-RLK genes in different tissues. Red and green colours correspond to up-regulation and down-regulation, respectively. Normalized gene expression values are provided in Additional file 1: Table S14. Figure S6. Co-expression network of strawberry LRR-RLK genes. Nodes indicate genes and edges indicate significant co-expression events between genes. Figure S7. Co-expression network of apple LRR-RLK genes. Nodes indicate genes and edges indicate significant co-expression events between genes. Figure S8. Co-expression network of Chinese white pear LRR-RLK genes. Nodes indicate genes and edges indicate significant co-expression events between genes. Figure S9. Co-expression network of mei LRR-RLK genes. Nodes indicate genes and edges indicate significant co-expression events between genes. Figure S10. Co-expression network of peach LRR-RLK genes. Nodes indicate genes and edges indicate significant co-expression events between genes. (PDF 1698Â kb)

Summary A highly saturated pear genetic map was constructed using over 3000 SNP markers developed by RADseq integrated with anchored SSR markers. The reliable QTLs of several fruit traits were also identified., Pear (Pyrus spp) is an important fruit crop, grown in all temperate regions of the world, with global production ranked after grape and apples among deciduous tree crops. A high-density linkage map is a valuable tool for fine mapping quantitative trait loci (QTL) and map-based gene cloning. In this study, we firstly constructed a high-density linkage map of pear using SNPs integrated with SSRs, developed by the rapid and robust technology of restriction-associated DNA sequencing (RADseq). The linkage map consists of 3143 SNP markers and 98 SSRs, 3241 markers in total, spanning 2243.4 cM, with an average marker distance of 0.70 cM. Anchoring SSRs were able to anchor seventeen linkage groups to their corresponding chromosomes. Based on this high-density integrated pear linkage map and two years of fruit phenotyping, a total of 32 potential QTLs for 11 traits, including length of pedicel (LFP), single fruit weight (SFW), soluble solid content (SSC), transverse diameter (TD), vertical diameter (VD), calyx status (CS), flesh colour (FC), juice content (JC), number of seeds (NS), skin colour (SC), and skin smooth (SS), were identified and positioned on the genetic map. Among them, some important fruit-related traits have for the first time been identified, such as calyx status, length of pedicel, and flesh colour, and reliable localization of QTLs were verified repeatable. This high-density linkage map of pear is a worthy reference for mapping important fruit traits, QTL identification, and comparison and combination of different genetic maps.