A self-compatible pear mutant derived from γ-irradiated pollen carries an 11-Mb duplication in chromosome 17

Self-compatibility is a highly desirable trait for pear breeding programs. Our breeding program previously developed a novel self-compatible pollen-part Japanese pear mutant (Pyrus pyrifolia Nakai), '415-1', by using γ-irradiated pollen. '415-1' carries the S-genotype S4dS5S5, with "d" indicating a duplication of S 5 responsible for breakdown of self-incompatibility. Until now, the size and inheritance of the duplicated segment was undetermined, and a reliable detection method was lacking. Here, we examined genome duplications and their inheritance in 140 F1 seedlings resulting from a cross between '515-20' (S1S3) and '415-1'. Amplicon sequencing of S-RNase and SFBB18 clearly detected S-haplotype duplications in the seedlings. Intriguingly, 30 partially triploid seedlings including genotypes S1S4dS5, S3S4dS5, S1S5dS5, S3S5dS5, and S3S4dS4 were detected among the 140 seedlings. Depth-of-coverage analysis using ddRAD-seq showed that the duplications in those individuals were limited to chromosome 17. Further analysis through resequencing confirmed an 11-Mb chromosome duplication spanning the middle to the end of chromosome 17. The duplicated segment remained consistent in size across generations. The presence of an S3S4dS4 seedling provided evidence for recombination between the duplicated S5 segment and the original S4haplotype, suggesting that the duplicated segment can pair with other parts of chromosome 17. This research provides valuable insights for improving pear breeding programs using partially triploid individuals.

Self S-RNase reduces the expression of two pollen-specific COBRA genes to inhibit pollen tube growth in pear

Due to self-incompatibility (SI) prevents self-fertilization, natural or artificial cross-pollination has been conducted in many orchards to stabilize fruit yield. However, it is still puzzled which routes of self S-RNase arresting pollen tube growth. Herein, 17 COBRA genes were isolated from pear genome. Of these genes, the pollen-specifically expressed PbCOB.A.1 and PbCOB.A.2 positively mediates pollen tube growth. The promoters of PbCOB.A.1 and/or PbCOB.A.2 were bound and activated by PbABF.E.2 (an ABRE-binding factor) and PbC2H2.K16.2 (a C2H2-type zinc finger protein). Notably, the expressions of PbCOB.A.1, PbCOB.A.2, and PbC2H2.K16.2 were repressed by self S-RNase, suggesting that self S-RNase reduces the expression of PbCOB.A.1 and PbCOB.A.2 by decreasing the expression of their upstream factors, such as PbC2H2.K16.2, to arrest pollen tube growth. PbCOB.A.1 or PbCOB.A.2 accelerates the growth of pollen tubes treated by self S-RNase, but can hardly affect level of reactive oxygen species and deploymerization of actin cytoskeleton in pollen tubes and cannot physically interact with any reported proteins involved in SI. These results indicate that PbCOB.A.1 and PbCOB.A.2 may not relieve S-RNase toxicity in incompatible pollen tube. The information provides a new route to elucidate the arresting pollen tube growth during SI reaction.

Contribution of Native Plasmids of Pantoea vagans C9-1 to Epiphytic Fitness and Fire Blight Management on Apple and Pear Flowers and Fruits

Pantoea vagans C9-1 (C9-1) is a biological control bacterium that is applied to apple and pear trees during bloom for suppression of fire blight, caused by Erwinia amylovora. Strain C9-1 has three megaplasmids: pPag1, pPag2, and pPag3. Prior bioinformatic studies predicted these megaplasmids have a role in environmental fitness and/or biocontrol efficacy. Plasmid pPag3 is part of the large Pantoea plasmid (LPP-1) group that is present in all Pantoea spp. and has been hypothesized to contribute to environmental colonization and persistence, while pPag2 is less common. We assessed fitness of C9-1 derivatives cured of pPag2 and/or pPag3 on pear and apple flowers and fruit in experimental orchards. We also assessed the ability of a C9-1 derivative lacking pPag3 to reduce populations of E. amylovora on flowers and disease incidence. Previously, we determined that tolerance to stresses imposed in vitro was compromised in derivatives of C9-1 lacking pPag2 and/or pPag3; however, in this study, the loss of pPag2 and/or pPag3 did not consistently reduce the fitness of C9-1 on flowers in orchards. Over the summer, pPag3 contributed to survival of C9-1 on developing apple and pear fruit in two of five trials, whereas loss of pPag2 did not significantly affect survival of C9-1. We also found that loss of pPag3 did not affect C9-1's ability to reduce E. amylovora populations or fire blight incidence on apple flowers. Our findings partially support prior hypotheses that LPP-1 in Pantoea species contributes to persistence on plant surfaces but questions whether LPP-1 facilitates host colonization.

Specific functions of single pistil S-RNases in S-gene homozygous Pyrus germplasm

Gametophytic self-incompatibility (SI) is regulated by S-allele recognition; that is, pollen in a style with the same S-genotype will undergo programmed cell death and stop growing so that it is unable to complete double fertilization, ultimately resulting in the SI response. S-RNase is the female determinant of SI in pear (Pyrus). In the Pyrus genome, there are two different S-RNase alleles at the S-locus, which generate two different S-RNase products in the pistil. The extracted S-glycoprotein is actually a protein complex. In this study, artificial self-pollination was conducted at the bud stage to overcome SI in 'Huanghua' (S1S2) pear. Seven plants homozygous for S1-RNase and four homozygous for S2-RNase were selected from the selfed progeny of 'Huanghua' by S-gene molecular identification biotechnology. We investigated the function of single S-RNases isolated from the pistils of S-gene homozygous Pyrus germplasm. The pollen of 'Huanghua' could smoothly pass through the style of the S-gene homozygous germplasm and complete fertilization. S-RNases were extracted from flower styles of different genotypes and used to treat different types of pollen. The S-RNase from 'Huanghua' completely inhibited the growth of S1S2, S1S1, and S2S2 pollen, while the S-RNase from homozygous germplasm allowed some S1S2 pollen and different single genotypes of pollen to continue growing. These results further validate the core events of SI including cytoskeleton depolymerization and programmed cell death. By iTRAQ-based proteomic analysis of style proteins, a total of 13 S-RNase-related proteins were identified. In summary, we have created reliable S-RNase gene homozygous germplasm, which will play a crucial role in further research on SI in pear and in the development of the pear industry.

Spatial Distribution Characteristics of Suitable Planting Areas for Pyrus Species under Climate Change in China

Planting suitability determines the distribution and yield of crops in a given region which can be greatly affected by climate change. In recent years, many studies have shown that carbon dioxide fertilization effects increase the productivity of temperate deciduous fruit trees under a changing climate, but the potential risks to fruit tree planting caused by a reduction in suitable planting areas are rarely reported. In this study, Maxent was first used to investigate the spatial distribution of five Pyrus species in China, and the consistency between the actual production area and the modeled climatically suitable area under the current climatic conditions were determined. In addition, based on Coupled Model Intercomparison Project Phase 6, three climate models were used to simulate the change in suitable area and the migration trend for different species under different emission scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5). The results showed that the suitable area for pear was highly consistent with the actual main production area under current climate conditions. The potential planting areas of P. ussuriensis showed a downward trend under all emission paths from 2020 to 2100; other species showed a trend of increasing first and then decreasing or slowing down and this growth effect was the most obvious in 2020-2040. Except for P. pashia, other species showed a migration trend toward a high latitude, and the trend was more prominent under the high emission path. Our results emphasize the response difference between species to climate change, and the method of consistency analysis between suitable planting area and actual production regions cannot only evaluate the potential planting risk but also provide a reasonable idea for the accuracy test of the modeled results. This work has certain guiding and reference significance for the protection of pear germplasm resources and the prediction of yield.

A large-scale proteogenomic atlas of pear

Pear is an important fruit tree that is widely distributed around the world. The first pear genome map was reported from our laboratory approximately 10 years ago. To further study global protein expression patterns in pear, we generated pear proteome data based on 24 major tissues. The tissue-resolved profiles provided evidence of the expression of 17 953 proteins. We identified 4294 new coding events and improved the pear genome annotation via the proteogenomic strategy based on 18 090 peptide spectra with peptide spectrum matches >1. Among the eight randomly selected new short coding open reading frames that were expressed in the style, four promoted and one inhibited the growth of pear pollen tubes. Based on gene coexpression module analysis, we explored the key genes associated with important agronomic traits, such as stone cell formation in fruits. The network regulating the synthesis of lignin, a major component of stone cells, was reconstructed, and receptor-like kinases were implicated as core factors in this regulatory network. Moreover, we constructed the online database PearEXP (http://www.peardb.org.cn) to enable access to the pear proteogenomic resources. This study provides a paradigm for in-depth proteogenomic studies of woody plants.

Sodicity stress differently influences physiological traits and anti-oxidant enzymes in pear and peach cultivars

Background

The growth and physiological responses to sodicity stress of pear and peach are poorly understood. Insights into how sodicity stress alters tree physiology remain vital to developing salt tolerant scion and rootstock cultivars.

Methods

The effects of sodicity stress (soil pH_s ~8.8) on tree growth and physiological traits of field grown trees of pear cultivars Punjab Beauty and Patharnakh, and peach cultivars Partap and Shan-e-Punjab were recorded using standard procedures. Sodicity-induced changes in oxidative stressors, proline, anti-oxidant enzymes and leaf ions were measured to draw inferences.

Results

Sodicity-induced reductions in vegetative growth were particularly marked in Patharnakh pear and Partap peach compared with other cultivars. Although sodicity stress triggered a significant increase in leaf malondialdehyde (MDA) and hydrogen peroxide (H₂O₂), their levels relative to controls were much higher in peach than in pear; reflecting that peach suffered from greater oxidative stress. Interestingly, MDA and H₂O₂ levels did not seem to be deleterious enough to trigger proline-induced osmotic adjustment in pears. The activities of anti-oxidant enzymes strongly varied with the cultivar; specifically, the sodicity-induced increases in CAT and SOD activities were much higher in Punjab Beauty pear and Shan-e-Punjab peach. Principal Component Analysis revealed an explicit convergence between CAT and SOD activities in Punjab Beauty and Shan-e-Punjab cultivars in response to sodicity-induced oxidative stress. Correlation analysis revealed that leaf Na⁺ strongly inhibited tree growth in peach than in pear. Leaf K⁺ and proline were found to be the major osmolytes in sodicity-stressed pear and peach cultivars, respectively.

Conclusions

We have for the first time studied the effects of sodicity stress on important tree growth and physiological traits of commercially important pear and peach cultivars. Our findings revealed a marked suppressive effect of sodicity stress on tree growth in peach than in pear. The sodicity-induced upticks in leaf malondialdehyde, hydrogen peroxide and Na⁺ seemed to induce proline-mediated osmotic adjustment in peach but not in pear. The overall better sodicity tolerance in pear compared to peach was ascribed to increased activities of anti-oxidant enzymes catalase and superoxide dismutase enzymes together with restricted Na⁺ uptake and better leaf K⁺ levels. Further investigations are needed to elucidate the effects of sodicity stress on genetic and transcriptional changes, and on fruit yield and quality.

Evaluation of Nutritional Value, Fatty Acids and Polyphenols Profiles of Pyrus amygdaliformis L. Grown in North-East Kurdistan Regional Government, Iraq

The present research evaluated some chemical composition in tissues from the leaves, seed, and fruits P. amygdaliformis. In leaf samples, phenolics compounds are predominated, whereas, in fruit samples, sugars, organic acid predominated. Among phenolic compounds, chlorogenic acid and fumaric acid were found in the highest amount in leaf tissue. In addition, the amounts of minerals in the fruit tissue had high values. In the P. Amigdaliformis seed oil, oleic and linoleic acid were the primary fatty acids. This study shows that these wild pears have beneficial natural compounds and can also be used as a new source of natural oil.

Analysis of the Fruit Quality of Pear (Pyrus spp.) Using Widely Targeted Metabolomics

Pear is a kind of common temperate fruit, whose metabolite composition that contributes to the difference in fruit quality is unclear. This study identified and quantified the metabolites using a widely targeted LC-MS/MS approach in three pear species, including Pyrus bretschneideri (PB), Pyrus usssuriensis (PU) and Pyrus pyrifolia (PP). A total of 493 metabolites were identified, consisting of 68 carbohydrates, 47 organic acids, 50 polyphenols, 21 amino acids, 20 vitamins, etc. The results of PCA and OPLS-DA demonstrated that the metabolite compositions differed distinctly with cultivar variability. Our results also involved some metabolic pathways that may link to the fruit quality based on KEGG pathway analysis, the pathway of phenylalanine metabolism revealed significant differences between PB and PP (p < 0.05). Furthermore, the study selected D-xylose, formononetin, procyanidin A1 and β-nicotinamide mononucleotide as the major differentially expressed metabolites in the three species. The present study can open new avenues for explaining the differences in fruit quality of the major commercial pear cultivars in China.

RNA-Seq Analysis Identifies Transcription Factors Involved in Anthocyanin Biosynthesis of 'Red Zaosu' Pear Peel and Functional Study of PpPIF8

Red-skinned pears are favored by people for their attractive appearance and abundance of anthocyanins. However, the molecular basis of anthocyanin biosynthesis in red pears remains elusive. Here, a comprehensive transcriptome analysis was conducted to explore the potential regulatory mechanism of anthocyanin biosynthesis in 'Red Zaosu' pear (Pyrus pyrifolia × Pyrus communis). Gene co-expression analysis and transcription factor mining identified 263 transcription factors, which accounted for 6.59% of the total number of transcription factors in the pear genome in two gene modules that are highly correlated with anthocyanin biosynthesis. Clustering, gene network modeling with STRING-DB, and local motif enrichment analysis (CentriMo) analysis suggested that PpPIF8 may play a role in anthocyanin biosynthesis. Furthermore, eight PIFs were identified in the pear genome, of which only PpPIF8 was rapidly induced by light. Functional studies showed that PpPIF8 localizes in the nucleus and is preferentially expressed in the tissue of higher levels of anthocyanin. The overexpression of PpPIF8 in pear peel and pear calli promotes anthocyanin biosynthesis and upregulates the expression of anthocyanin biosynthesis genes. Yeast-one hybrid and transgenic analyses indicated that PpPIF8 binds to the PpCHS promoter to induce PpCHS expression. The positive effect of PpPIF8 on anthocyanin biosynthesis is different from previously identified negative regulators of PyPIF5 and MdPIF7 in pear and apple. Taken together, our data not only provide a comprehensive view of transcription events during the coloration of pear peel, but also resolved the regulatory role of PpPIF8 in the anthocyanin biosynthesis pathway.

Variations in Pedicel Structural Properties Among Four Pear Species (Pyrus): Insights Into the Relationship Between the Fruit Characteristics and the Pedicel Structure

Fruit pedicel is the bridge linking the parent tree and the fruit, which is an important channel for water and nutrients transport to the fruit. The genetic specificity determines the characteristics of the pedicel and the fruit, but the relationship between the pedicel structure and the fruit characteristics is unexplored. Combining the investigation of fruit characteristics, the statistical analysis of the pedicel structural properties, and the 2D and 3D anatomical observation of the pedicel, this study found distinctive contributions of the pedicel elements to the fruit characteristics in four pear species. The European pear (Conference) showed distinct fruit shape index and pedicel structural properties compared with the oriental pears (Akizuki, Yali, and Nanguoli). The fruit size positively correlated with pedicel length, fiber area, pedicel diameter, the area percentage of the cortex, and the area percentage of phloem; however, fruit firmness and soluble solids concentration are showed a stronger positive correlation with xylem area, pith area, the area percentage of xylem, the area percentage of sieve tube, and the area percentage of pith. Pedicel elements, including pith, fiber, and cortex, likely play a certain role in the fruit growth due to the variations of their characteristics demonstrated in the four pear species. The porosity, the ratio of the surface area to the volume, and the spatial arrangement of the vessels showed significant variations across the pear species, indicating the distinction of the hydraulic conductance of the pedicels. Our findings provided direct evidence that pedicel structural elements contributed distinctively to the fruit characteristics among pear species.

Genome-wide identification and functional analysis of U-box E3 ubiquitin ligases gene family related to drought stress response in Chinese white pear (Pyrus bretschneideri)

BACKGROUND

The plant U-box (PUB) proteins are a family of ubiquitin ligases (E3) enzymes that involved in diverse biological processes, as well as in responses to plant stress response. However, the characteristics and functional divergence of the PUB gene family have not yet been previously studied in the Chinese white pear (Pyrus bretschneideri).

RESULTS

In the present study, we identified 62 PbrPUBs in Chinese white pear genome. Based on the phylogenetic relationship, 62 PUB genes were clustered into five groups. The results of conserved motif and gene structure analysis supported the classification phylogenetic tree. The PbrPUB genes were unevenly distribution on 17 pear chromosomes, chromosome 15 housed most member of PUB family, with eight PUB genes. Cis-acting element analysis indicated that PUB genes might participate in diverse biological processes, especially in the response to abiotic stresses. Based on RNA-data from 'Dangshansuli' at seven tissues, we found that PUB genes exhibited diverse of expression level in seven tissues, and qRT-PCR experiment further supported the reliable of RNA-Seq data. To identify candidate genes associated with resistance, we conducted qRT-PCR experiment the expression level of pear seed plant under four abiotic stresses, including: ABA, dehydration, salt and cold treatment. One candidate PUB gene associated with dehydration stress was selected to conduct further functional experiment. Subcellular localization revealed PbrPUB18 protein was located on cell nucleus. Furthermore, heterologous over-expression of PbrPUB18 in Arabidopsis indicated that the over-expression of PbrPUB18 could enhance resistance in drought treatment. In conclusions, we systematically identified the PUB genes in pear, and provided useful knowledge for functional identification of PUB genes in pear.

Nitrendipine-Treatment Increases Cork Spot Disorder Incidence in Pear 'Akituki' (Pyrus pyrifolia Nakai.) by Altering Calcium Distribution Inside the Fruit

‘Akituki’ (Pyrus pyrifolia Nakai.) is a very popular and profitable pear cultivar in China. However, its high susceptibility to cork spot disorder has limited its expansion of cultivated area. The mechanisms of cork spot disorder have been discussed extensively, focusing on Ca²⁺ deficiency, yet no consensus has been made. In this study, we applied nitrendipine (NI) as a Ca²⁺ uptake inhibitor to explore the role of calcium in cork spot disorder occurrence. Results showed that NI treatment on the fruit remarkably increased the incidence of cork spot disorder; alteration of mineral contents happened at the early developmental stage of the fruit, especially on the outer flesh and the peel of the fruit; and this gap was filled gradually along with the expansion of the fruit. Significant differences in the expression levels of Ca²⁺ transport-related genes were found in the inner flesh, outer flesh and peel during the fruit growth period. The observation of free Ca²⁺ localization indicated the intracellular imbalance of Ca²⁺ in the NI-treated fruit. In conclusion, NI treatment reduced the calcium content in the fruit at an early developmental stage, altered the related expression of genes and influenced the cellular Ca²⁺ balance in the fruit, which prompted the occurrence of cork spot disorder. Measures for the prevention and control of cork spot disorder should be taken at the early stage of the fruit development in the field.

Homozygosity Mapping Reveals Population History and Trait Architecture in Self-Incompatible Pear (Pyrus spp.)

Runs of homozygosity (ROH) have been widely used to study population history and trait architecture in humans and livestock species, but their application in self-incompatible plants has not been reported. The distributions of ROH in 199 accessions representing Asian pears (45), European pears (109), and interspecific hybrids (45) were investigated using genotyping-by-sequencing in this study. Fruit phenotypes including fruit weight, firmness, Brix, titratable acidity, and flavor volatiles were measured for genotype-phenotype analyses. The average number of ROH and the average total genomic length of ROH were 6 and 11 Mb, respectively, in Asian accessions, and 13 and 30 Mb, respectively, in European accessions. Significant associations between genomic inbreeding coefficients (F_ROH) and phenotypes were observed for 23 out of 32 traits analyzed. An overlap between ROH islands and significant markers from genome-wide association analyses was observed. Previously published quantitative trait loci for fruit traits and disease resistances also overlapped with some of the ROH islands. A prominent ROH island at the bottom of linkage group 17 overlapped with a recombination-supressed genomic region harboring the self-incompatibility locus. The observed ROH patterns suggested that systematic breeding of European pears would have started earlier than of Asian pears. Our research suggest that F_ROH would serve as a novel tool for managing inbreeding in gene-banks of self-incompatible plant species. ROH mapping provides a complementary strategy to unravel the genetic architecture of complex traits, and to evaluate differential selection in outbred plants. This seminal work would provide foundation for the ROH research in self-incompatible plants.

Mining and evolution analysis of lateral organ boundaries domain (LBD) genes in Chinese white pear (Pyrus bretschneideri)

BACKGROUND

The lateral organ boundaries domain (LBD) gene is a plant-specific transcription factor that plays a critical role in diverse biological processes. However, the evolution and functional divergence of the LBD gene family has not yet been characterized for the Chinese White Pear.

RESULTS

In our study, a total of 60 PbrLBDs were identified in the pear genome. The PbrLBD gene family was divided into two classes based on gene structure and phylogenetic analysis: class I (53) and class II (7). Cis-acting element analysis results suggested that PbrLBDs may participate in various biological processes, such as flavonoid biosynthetic and stress response. Synteny analysis results indicated that segmental duplication played a key role in the expansion of the PbrLBD gene family. The mean Ks and 4DTv values showed that the PbrLBD gene family had undergone only one recent whole-genome duplication event occurring at 30-45 MYA. Purifying selection was a primary force during the PbrLBD gene family evolution process. Transcriptome data analysis revealed that 10 PbrLBDs were expressed in all six examined tissues, and 73.33% of members in the PbrLBD gene family were expressed in pear sepal. qRT-PCR was conducted to verify the expression levels of 11 PbrLBDs in these six tissues. Specifically, PbrLBD20, PbrLBD35 and PbrLBD53 genes were down-regulated when anthocyanin concentrations were high, whereas PbrLBD33 was significantly up-regulated in pear when anthocyanin concentrations were high. Furthermore, PbrLBD20, one of the candidate genes related to anthocyanins was localized in the nucleus.

CONCLUSIONS

Our analysis provides valuable information for understanding the evolution of the PbrLBD gene family, and provides new insights into the regulation of pear pigment metabolism and lays a foundation for the future disclosure of the molecular mechanism of LBD gene regulating flavonoid metabolism.

Characterization of the Auxin Efflux Transporter PIN Proteins in Pear

PIN-FORMED (PIN) encodes a key auxin polar transport family that plays a crucial role in the outward transport of auxin and several growth and development processes, including dwarfing trees. We identified a dwarfing pear rootstock 'OHF51' (Pyrus communis), which limits the growth vigor of the 'Xueqing' (Pyrus bretschneideri × Pyrus pyrifolia) scion, and isolated 14 putative PbPINs from the pear Pyrus bretschneideri. The phylogenic relationships, structure, promoter regions, and expression patterns were analyzed. PbPINs were classified into two main groups based on the protein domain structure and categorized into three major groups using the neighbor-joining algorithm. Promoter analysis demonstrated that PbPINs might be closely related to plant growth and development. Through quantitative real-time PCR (qRT-PCR) analysis, we found that the expression patterns of 14 PbPINs varied upon exposure to different organs in dwarfing and vigorous stocks, 'OHF51' and 'QN101' (Pyrus betulifolia), indicating that they might play varying roles in different tissues and participated in the regulation of growth vigor. These results provide fundamental insights into the characteristics and evolution of the PINs family, as well as the possible relationship between dwarfing ability and auxin polar transport.

High diversity of Diaporthe species associated with pear shoot canker in China

Species of Diaporthe (syn. Phomopsis) are important endophytes, saprobes and pathogens, infecting a wide range of plants and resulting in important crop diseases. However, the species occurring on pear remain largely unresolved. In this study, a total of 453 Diaporthe isolates were obtained from branches of Pyrus plants (including P. bretschneideri, P. communis, P. pyrifolia and P. ussuriensis collected from 12 provinces in China) showing shoot canker symptoms. Phylogenetic analyses based on five loci (ITS, TEF, CAL, HIS, and TUB) coupled with morphology of 113 representative isolates revealed that 19 Diaporthe species were isolated, representing 13 known species (including D. caryae, D. cercidis, D. citrichinensis, D. eres, D. fusicola, D. ganjae, D. hongkongensis, D. padina, D.pescicola, D. sojae, D. taoicola, D. unshiuensis and D. velutina) and six new species described here as D. acuta, D. chongqingensis, D. fulvicolor, D. parvae, D. spinosa and D. zaobaisu. Although Koch’s postulates confirmed all species to be pathogenic, a high degree of variation in aggressiveness was observed. Moreover, these species have a high diversity, plasticity, and prevalence related to the geographical location and pear species involved.

The complete plastome sequence of a subtropical tree Pyrus betulaefolia (Rosaceae)

The genus Pyrus, comprising several popular fruit crops worldwide, includes over 30 tree species. Here we determined the complete plastid genome sequence of Pyrus betulaefolia. The plastome consists of 160,184 bp, including a pair of inverted repeats (IRs) with a length of 26,384 bp separated by a large single-copy region (LSC) and a small single-copy region (SSC) of 88,121 bp and 19,295 bp, respectively. Further phylogenetic analyze was conducted using 11 complete plastid genomes of Rosaceae with KVM + F + I model, which supports Pyrus betulaefolia as a sister to all other eight Pyrus taxa with published plastomes.

Evolutionary Rate Heterogeneity and Functional Divergence of Orthologous Genes in Pyrus

Negatively selected genes (NSGs) and positively selected genes (PSGs) are the two types of most nuclear protein-coding genes in organisms. However, the evolutionary rates and characteristics of different types of genes have been rarely understood. In the present study, we investigate the rates of synonymous substitution (Ks) and the rates of non-synonymous substitution (Ka) by comparing the orthologous genes of two sequenced Pyrus species, Pyrus bretschneideri and Pyrus communis. Subsequently, we compared the evolutionary rates, gene structures, and expression profiles during different fruit development between PSGs and NSGs. Compared with the NSGs, the PSGs have fewer exons, shorter gene length, lower synonymous substitution rates and have higher evolutionary rates. Remarkably, gene expression patterns between two Pyrus species fruit indicated functional divergence for most of the orthologous genes derived from a common ancestor, and subfunctionalization for some of them. Overall, the present study shows that PSGs differs from NSGs not only under environmental selective pressure (Ka/Ks), but also in their structural, functional, and evolutionary properties. Additionally, our resulting data provides important insights for the evolution and highlights the diversification of orthologous genes in two Pyrus species.

Development of a highly efficient Axiom™ 70 K SNP array for Pyrus and evaluation for high-density mapping and germplasm characterization

Background

Both a source of diversity and the development of genomic tools, such as reference genomes and molecular markers, are equally important to enable faster progress in plant breeding. Pear (Pyrus spp.) lags far behind other fruit and nut crops in terms of employment of available genetic resources for new cultivar development. To address this gap, we designed a high-density, high-efficiency and robust single nucleotide polymorphism (SNP) array for pear, with the main objectives of conducting genetic diversity and genome-wide association studies.

Results

By applying a two-step design process, which consisted of the construction of a first ‘draft’ array for the screening of a small subset of samples, we were able to identify the most robust and informative SNPs to include in the Applied Biosystems™ Axiom™ Pear 70 K Genotyping Array, currently the densest SNP array for pear. Preliminary evaluation of this 70 K array in 1416 diverse pear accessions from the USDA National Clonal Germplasm Repository (NCGR) in Corvallis, OR identified 66,616 SNPs (93% of all the tiled SNPs) as high quality and polymorphic (PolyHighResolution). We further used the Axiom Pear 70 K Genotyping Array to construct high-density linkage maps in a bi-parental population, and to make a direct comparison with available genotyping-by-sequencing (GBS) data, which suggested that the SNP array is a more robust method of screening for SNPs than restriction enzyme reduced representation sequence-based genotyping.

Conclusions

The Axiom Pear 70 K Genotyping Array, with its high efficiency in a widely diverse panel of Pyrus species and cultivars, represents a valuable resource for a multitude of molecular studies in pear. The characterization of the USDA-NCGR collection with this array will provide important information for pear geneticists and breeders, as well as for the optimization of conservation strategies for Pyrus.

Electronic supplementary material

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

Characterization of genome-wide microRNAs and their roles in development and biotic stress in pear

Using a genome-wide analysis of miRNAs in 'Yali' pear (Pyrus bretschneideri) via the next-generation high-throughput sequencing of small RNAs with a bioinformatics analysis, we found that pbr-miR156, pbr-miR164, pbr-miR399, and pbr-miR482 and their target genes function in viral defense in 'Duli' and 'Hongbaoshi'. pbr-miR160, pbr-miR168, pbr-miR171, and pbr-miR319 and their targets function in auxin signaling pathways in 'Zhongai 4' and 'Zhongai 5'. Successful fruit production in pear (Pyrus spp.) depends on the use of optimal combinations of rootstocks and scions. Deciphering plant-pathogen defense mechanisms and hormone signaling pathways is an important step towards developing pear rootstocks and varieties with improved qualities. In the current study, we combined next-generation sequencing of small RNAs with a bioinformatics analysis to systematically identify and characterize 298 miRNAs in the pear scion cultivar 'Yali' (Pyrus bretschneideri). We also analyzed miRNAs in three rootstock varieties ('Duli', 'Zhongai 4', and 'Zhongai 5') and one scion cultivar ('Hongbaoshi'). We found that pbr-miR156, pbr-miR164, pbr-miR399, and pbr-miR482 are induced following infection with the pear virus Apple stem pitting virus (ASPV), and identified their target genes (pbRPS6, pbNAC, pbTLR, and pbRX-CC, respectively), which participate in viral defense pathways in 'Duli' and 'Hongbaoshi'. Furthermore, we identified pbr-miR160, pbr-miR168, pbr-miR171, and pbr-miR319, and found that the production of these miRNAs was suppressed under low levels of synthetic auxin. The targets of these miRNAs (pbARF, pbAEC, pbSCL, and pbTCP4) respond to auxin signaling pathways in 'Zhongai 4' and 'Zhongai 5'. Our results lay the foundation for breeding improved pear cultivars.

Genome-wide identification, evolution, and expression analysis of the KT/HAK/KUP family in pear

The K⁺ transporter/high-affinity K⁺/K⁺ uptake (KT/HAK/KUP) family, as one of the largest K⁺ transporter families in higher plants, plays an essential role in plant growth, mineral element absorption, salt stress tolerance, and other physiological processes. However, little is known about this family in pear (Pyrus). Here, we identified 20 K⁺ transporter genes in pear (P. bretschneideri) using genome-wide analysis. Their gene structure, chromosomal distribution, conserved motifs, phylogenetics, duplication events, and expression patterns were also examined. The results of phylogenetic analysis showed that PbrKT/HAK/KUP genes were clustered into three major groups (Groups I-III). Among the 20 PbrKT/HAK/KUP genes, 18 were mapped to nine chromosomes and two to scaffolds. Four WGD/segmental gene pairs were identified, indicating that WGD/segmental duplication may have contributed to the expansion of the KT/HAK/KUP family in pear. Among the four pairs of WGD/segmentally duplicated genes, both members of three pairs had been subjected to purifying selection, whereas the fourth pair had been subjected to positive selection. Furthermore, phenotypic experiments showed that the growth of pear seedlings was affected by potassium deficiency treatment. Expression patterns of 20 PbrKT/HAK/KUP genes in roots were further assayed with qRT-PCR. PbrHAK1 and PbrHAK12/16 were significantly expressed in response to K⁺ deficiency, suggesting that these genes are crucial for K⁺ uptake in pear, especially under the condition of K⁺ starvation. Our results provide a foundation for further study on the function of KT/HAK/KUP genes in pear.

Comparative Analysis of Transcriptomes to Identify Genes Associated with Fruit Size in the Early Stage of Fruit Development in Pyrus pyrifolia

Pear (Pyrus L.) is an important commercial fruit in the world. The fruit size is one of the important characters in fruit quality. The previous research reported that the fruit size of pear was mainly caused by the number of cell in about 40 days after blossom (DAB) in nature. However, studies about the mechanisms underlying cell division in young fruit development are very limited in pear. Two pear accessions codenamed ‘GH59B’ with big fruit and ‘GH81S’ with small fruit in three stages were sampled and the RNA-seq high-throughput sequencing was used to evaluate changes of gene transcription levels in the early stage of fruit development. The difference of cell size among two samples was little in 40 DAB, implying that the difference of the fruit size was caused by the number of the cell. More than 274,517,982 high quality reads from six libraries of fruit development were sequenced. A total of 797 differentially expressed genes (DEGs) were identified. Three cytokinin dehydrogenase genes and two gibberellin 2-beta-dioxygenase gene were identified in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways related to zeatin and gibberellin. Their expression was upregulated at 20 DAB in ‘GH81S’ and at 30 DAB in ‘GH59B’, suggesting that the small fruit size might be related to the early degradation of cytokinin and gibberellin inducing a short period of cell division. A total of 38 DEGs of transcription factors were found and 23 DEGs including NAC, ERF and bHLH transcription factors were highly related with cytokinin dehydrogenase and gibberellin dioxygenase genes. Altogether, the results of the present study provide information from a comprehensive gene expression analysis and insight into the molecular mechanism underlying the difference of fruit size in Pyrus pyrifolia.

Colletotrichum species associated with anthracnose of Pyrus spp. in China

Colletotrichum species are plant pathogens, saprobes, and endophytes on a range of economically important hosts. However, the species occurring on pear remain largely unresolved. To determine the morphology, phylogeny and biology of Colletotrichum species associated with Pyrus plants, a total of 295 samples were collected from cultivated pear species (including P. pyrifolia, P. bretschneideri, and P. communis) from seven major pear-cultivation provinces in China. The pear leaves and fruits affected by anthracnose were sampled and subjected to fungus isolation, resulting in a total of 488 Colletotrichum isolates. Phylogenetic analyses based on six loci (ACT, TUB2, CAL, CHS-1, GAPDH, and ITS) coupled with morphology of 90 representative isolates revealed that they belong to 10 known Colletotrichum species, including C. aenigma, C. citricola, C. conoides, C. fioriniae, C. fructicola, C. gloeosporioides, C. karstii, C. plurivorum, C. siamense, C. wuxiense, and two novel species, described here as C. jinshuiense and C. pyrifoliae. Of these, C. fructicola was the most dominant, occurring on P. pyrifolia and P. bretschneideri in all surveyed provinces except in Shandong, where C. siamense was dominant. In contrast, only C. siamense and C. fioriniae were isolated from P. communis, with the former being dominant. In order to prove Koch’s postulates, pathogenicity tests on pear leaves and fruits revealed a broad diversity in pathogenicity and aggressiveness among the species and isolates, of which C. citricola, C. jinshuiense, C. pyrifoliae, and C. conoides appeared to be organ-specific on either leaves or fruits. This study also represents the first reports of C. citricola, C. conoides, C. karstii, C. plurivorum, C. siamense, and C. wuxiense causing anthracnose on pear.

Comparative and Expression Analysis of Ubiquitin Conjugating Domain-Containing Genes in Two Pyrus Species

Ripening affects the nutritional contents and quality of fleshy fruits, and it plays an important role during the process of fruit development. Studies have demonstrated that ubiquitin-conjugating (UBC or E2) genes can regulate fruit ripening, but the characterization of UBCs in pear is not well documented. The recently published genome-wide sequences of Pyrus bretschneideri and Pyrus communis have allowed a comprehensive analysis of this important gene family in pear. Using bioinformatics approaches, we identified 83 (PbrUBCs) and 84 (PcpUBCs) genes from P. bretschneideri and P. communis, respectively, which were divided into 13 subfamilies. In total, 198 PbrUBC paralogous, 215 PcpUBC paralogous, and 129 orthologous gene pairs were detected. Some paralogous gene pairs were found to be distributed on the same chromosome, suggesting that these paralogs may be caused by tandem duplications. The expression patterns of most UBC genes were divergent between Pyrus bretschneideri and Pyrus communis during pear fruit development. Remarkably, the transcriptome data showed that UBC genes might play a more important role in fruit ripening for further study. This is the first report on the systematic analysis of two Pyrus UBC gene families, and these data will help further study the role of UBC genes in fruit development and ripening, as well as contribute to the functional verification of UBC genes in pear.

Pest categorisation of Botryosphaeria kuwatsukai

The Panel on Plant Health performed a pest categorisation of Botryosphaeria kuwatsukai, the causal agent of fruit rot and wart bark on apple and pear, for the EU. The pathogen, which was recently characterised, is a well-defined fungal species affecting mainly Pyrus pyrifolia (Japanese pear), although Pyrus communis (European pear) and apples (Malus domestica) can also be affected. The host status of other plant species reported in the literature, i.e. Cydonia oblonga, Chaenomeles japonica, Malus micromalus, Vitis vinifera and Prunus spp., is unclear. B. kuwatsukai is currently present in Japan, China, Korea, Taiwan and the USA, and uncertainty exists about its presence in other areas, where the disease has been associated with other Botryosphaeria spp. The pathogen is not known to occur in the EU and is listed in Annex IIAI of Directive 2000/29/EC. It could potentially enter the EU on host plants for planting and fruit originated in infested countries. Climatic conditions in the EU are suitable for the establishment and spread of the pathogen, as its epidemiology is similar to that of other Botryosphaeria spp. present in the EU. Pears and apples are widely distributed in the EU. In the infested areas, B. kuwatsukai causes branch dieback and fruit rot resulting in yield/quality losses. Its introduction and spread in the EU could impact pear and apple production, although the magnitude is unknown. Cultural practices and chemical measures may reduce the inoculum sources but cannot eliminate the pathogen. Phytosanitary measures are available to mitigate the risk of introduction and spread of the pathogen in the EU. B. kuwatsukai meets all criteria assessed by EFSA for consideration as a potential Union quarantine pest. As B. kuwatsukai is not known to occur in the EU, this criterion to consider it as a Union regulated non-quarantine pest is not met.

The role of polyamines during exocarp formation in a russet mutant of 'Dangshansuli' pear (Pyrus bretschneideri Rehd.)

Differential genes of suberin, polyamine and transcription factors in transcriptome sequences and the contents of H 2 O 2 , spermidine, spermine, and putrescine changed significantly after treating with MGBG. Russeting is a commercially important process that restores the control of water loss through the skin via the formation of a waterproofing periderm just beneath the microcracked skin of pear primary fruit. A spontaneous russet skin mutant, the yellow-green 'Dangshansuli' pear, has been identified. To understand the role of polyamines in the formation of the russet skin of the mutant-type (MT) pear, it was treated with methylglyoxal-bis-(guanylhydrazone) (MGBG) for 4 weeks after full bloom. One week later, differentially expressed genes among the wild-type (WT), MT, and MGBG-treated MT pears were screened, hydrogen peroxide (H2O2) was localized using CeCl3, and the contents of H2O2 and polyamine were measured. A total of 57,086,772, 61,240,014, and 67,919,420 successful reads were generated from the transcriptomes of WT, MT, and MGBG-treated MT, with average unigene lengths of 701, 720, and 735 bp, respectively. Differentially expressed genes involved in polyamine metabolism and suberin synthesis were screened in 'Dangshansuli' and in the mutant libraries, and their relative expression was found to be significantly altered after treatment with MGBG, which was confirmed by real-time PCR. The expression patterns of differentially expressed transcription factors were identified and were found to be similar to those of the polyamine- and suberin-related genes. The results indicated that the H2O2 generated during polyamine metabolism might contribute to russet formation on the exocarp of the mutant pear. Furthermore, the contents of H2O2, spermidine, spermine, and putrescine and H2O2 localization provided a comprehensive transcriptomic view of russet formation in the mutant pear.

Isolation and characterization of putative functional long terminal repeat retrotransposons in the Pyrus genome

Background

Long terminal repeat (LTR)-retrotransposons constitute 42.4 % of the genome of the ‘Suli’ pear (Pyrus pyrifolia white pear group), implying that retrotransposons have played important roles in Pyrus evolution. Therefore, further analysis of retrotransposons will enhance our understanding of the evolutionary history of Pyrus.

Results

We identified 1836 LTR-retrotransposons in the ‘Suli’ pear genome, of which 440 LTR-retrotransposons were predicted to contain at least two of three gene models (gag, integrase and reverse transcriptase). Because these were most likely to be functional transposons, we focused our analyses on this set of 440. Most of the LTR-retrotransposons were estimated to have inserted into the genome less than 2.5 million years ago. Sequence analysis showed that the reverse transcriptase component of the identified LTR-retrotransposons was highly heterogeneous. Analyses of transcripts assembled from RNA-Seq databases of two cultivars of Pyrus species showed that LTR-retrotransposons were expressed in the buds and fruit of Pyrus. A total of 734 coding sequences in the ‘Suli’ genome were disrupted by the identified LTR-retrotransposons. Five high-copy-number LTR-retrotransposon families were identified in Pyrus. These families were rarely found in the genomes of Malus and Prunus, but were distributed extensively in Pyrus and abundance varied between species.

Conclusions

We identified potentially functional, full-length LTR-retrotransposons with three gene models in the ‘Suli’ genome. The analysis of RNA-seq data demonstrated that these retrotransposons are expressed in the organs of pears. The differential copy number of LTR-retrotransposon families between Pyrus species suggests that the transposition of retrotransposons is an important evolutionary force driving the genetic divergence of species within the genus.

Electronic supplementary material

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

Phylogeny and evolutionary histories of Pyrus L. revealed by phylogenetic trees and networks based on data from multiple DNA sequences

Reconstructing the phylogeny of Pyrus has been difficult due to the wide distribution of the genus and lack of informative data. In this study, we collected 110 accessions representing 25 Pyrus species and constructed both phylogenetic trees and phylogenetic networks based on multiple DNA sequence datasets. Phylogenetic trees based on both cpDNA and nuclear LFY2int2-N (LN) data resulted in poor resolution, especially, only five primary species were monophyletic in the LN tree. A phylogenetic network of LN suggested that reticulation caused by hybridization is one of the major evolutionary processes for Pyrus species. Polytomies of the gene trees and star-like structure of cpDNA networks suggested rapid radiation is another major evolutionary process, especially for the occidental species. Pyrus calleryana and P. regelii were the earliest diverged Pyrus species. Two North African species, P. cordata, P. spinosa and P. betulaefolia were descendent of primitive stock Pyrus species and still share some common molecular characters. Southwestern China, where a large number of P. pashia populations are found, is probably the most important diversification center of Pyrus. More accessions and nuclear genes are needed for further understanding the evolutionary histories of Pyrus.

Diversity of odor-active compounds from local cultivars and wild accessions of Iwateyamanashi (Pyrus ussuriensis var. aromatica) revealed by Aroma Extract Dilution Analysis (AEDA)

Some local cultivars and wilds of Iwateyamanashi (Pyrus ussuriensis var. aromatica) that grows wild in Northern Tohoku, Japan have good aromatic fruit. Iwateyamanashi may be valuable germplasms as a donor of odor compounds in breeding of Japanese pear (Pyrus pyrifolia), because almost all Japanese pear cultivars have faint odor. Fruits odors from a local cultivar 'Sanenashi', a wild accession (i0830) in Iwateyamanashi, cultivars of 'Kosui' and 'La France' were characterized at first with comparative Aroma Extract Dilution Analysis (AEDA). Application of AEDA, based on Gas chromatography/Olfactometry analysis (GC/O), on the odor concentration prepared from 'Sanenashi' indicated the presence of 33 odor-active compounds including methyl and ethyl esters, aldehydes and alcohol. The eleven odor compounds from 16 accessions of Iwateyamanashi showed various combinations and wide range of odor concentrations by Principal Component Analysis (PCA). Especially 2 accessions of local cultivar 'Natsunashi' plotted in the highly ethyl ester group might be useful for Japanese pear breeding.

Development of cultivar-specific DNA markers based on retrotransposon-based insertional polymorphism in Japanese pear

We developed retrotransposon-based insertional polymorphism (RBIP) markers based on the long terminal repeat (LTR) sequences of copia-like retrotransposon Ppcrt4 and flanking genome sequences, which were derived from 454 sequencing data from Japanese pear (Pyrus pyrifolia) 'Hosui'. Out of 40 sequences including both LTR and flanking genome regions, we developed 22 RBIP markers and used them for DNA profiling of 80 pear cultivars: 64 Japanese, 10 Chinese (Pyrus ussuriensis) and 6 European (Pyrus communis). Three RBIP markers were enough to differentiate 'Hosui' from the other Japanese pear cultivars. The 22 RBIP markers could also distinguish 61 of the 64 Japanese pear cultivars. European pears showed almost no amplification of the 22 RBIP markers, which might suggest that retrotransposons had transposed during Asian pear evolution or reflect the genetic relationship between Asian and European pears. Sixteen of the RBIP markers could be positioned on a genetic linkage map of 'Hosui'. The RBIP loci were distributed in 10 linkage groups, and some loci were very closely located within the same linkage group. The information obtained will be applicable to developing cultivar-specific RBIP marker sets in plants.

Host Suitability of Eight Prunus spp. and One Pyrus communis Rootstocks to Pratylenchus vulnus, P. neglectus, and P. thornei

The effects of Pratylenchus vulnus on rootstocks of eight commonly used Prunus spp. and one Pyrus communis were evaluated under greenhouse conditions during a 15-month period. In a first experiment, two almonds (Moncayo and Garrigues), one peach (GF-305), and two peach-almond hybrids (GF-677 and Adafuel) inoculated with 2,000 nematodes per plant proved to be good hosts of P. vulnus. Highest (P < 0.05) numbers of nematodes per gram of fresh root weight were recovered from Adafuel and GF-677. Root weights were higher in uninoculated compared to inoculated plants of all rootstocks, whereas top weights of uninoculated Garrigues, GF-305, and GF-677 differed (P < 0.05) from those of inoculated plants. In a second experiment, three plum (Marianna 2624, Myrobalan 605, and San Julian 655-2) and one pear (OHF-333) rootstocks were also found to be good hosts of P. vulnus, although significantly fewer nematodes were recovered from Myrohalan 605 roots than from the other three materials. Inoculated OHF-333 and San Julian 655-2 differed (P < 0.05) in root weights over uninoculated plants. Only inoculated San Julian 655-2 showed differences in top weights over uninoculated treatments. Rootstocks were poor or non-hosts for P. neglectus and P. thornei.